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Introduction Générale 


Introduction : 

Dans la plupart des applications industrielles, le développement d’un modèle 
mathématique exacte décrivant le mieux le comportement réel d’un système réel a toujours 
été un défi à cause de la diversité et la complexité des systèmes physiques à commander. En 
revanche, le modèle de synthèse est souvent entaché par des incertitudes de modélisation [1] 
ce qui peut restreindre le domaine d’application des contrôleurs standards développés par des 


méthodes conventionnelles de synthèse. 


Au fil des années, grâce à ses propriétés et à sa simplicité de mise en œuvre, la commande 
H à renforcé sa présence en termes d’applications dans le milieu industriel. C’est une 
méthode très efficace pour résoudre le problème de la commande des systèmes physiques 
basés sur la modélisation incertaine. Ceci par la satisfaction de plus des spécifications de 
performance et de robustesse, non seulement pour le modèle nominal de synthèse, mais 
également en présence de diversités des effets indésirables qui sont causés par l’existence des 
incertitudes de modélisation qui sont inévitables lors de la modélisation du comportement réel 
du système physique, l’existence des modes inconnus qui doivent les négliger lors de la 
linéarisation de la dynamique non linéaire du système autour d’un ou plusieurs points de 
fonctionnement, l’existence des bruits de mesure causés par la fatigue des instruments de 
mesures installés dans la boucle de commande,...etc. Un autre avantage de cette méthode de 
commande c’est qu’on peut aussi l’appliquer sur plusieurs types des systèmes. Parmi les types 
les plus saillants, on peut citer les processus simples, complexes, mono-entrée et mono-sortie 
SISO (Single Input Single output) et les systèmes multi-entrées, multi-sorties MIMO (Multi 
Input Multi Output). 


En général, le concept de cette stratégie de commande s’appuie préalablement sur la 
représentation du modèle incertain sous une forme telle que les incertitudes soient extraites du 
modèle nominal. Ces incertitudes peuvent être classées en deux catégories on l’occurrence les 
incertitudes structurées et les incertitudes non structurées [1,2,3,4]. Pour les incertitudes 
structurées, elles sont dues souvent par la variation des certains paramètres intervenant dans la 
structure du modèle de synthèse. Sachant que leur effet est important notamment en basse 
fréquence. C’est pour cette raison qu’on doit imposer des bornes minimales et maximales 
limitant chaque paramètre incertain du modèle. En revanche, les incertitudes non structurées 


quand elles agissent sur le modèle de synthèse d’une manière globale sans pouvoir quantifier 
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précisément le phénomène. Elles sont dues souvent par des approximations effectuées sur le 
retard pur du modèle ainsi que la simplification d’un certain nombre de modes inconnus lors 
de la linéarisation de la dynamique réel du système physique. Sachant que leur effet est 
important en haute fréquence. C’est pourquoi qu’on doit imposer un majorant limitant toutes 


les incertitudes susceptibles qui affectent le modèle de synthèse [5]. 


En présence des incertitudes de modélisation, la première question concerne la stabilité du 
système bouclé. A cet effet, l’analyse de la robustesse en stabilité consiste donc à établir si le 
système corrigé demeure stable ou non. Dans ce cas on cherche à satisfaire une condition de 
robustesse appelée « condition sur la stabilité robuste RS » [4,6, 7, 8]. De plus, on peut aussi 
désirer que le modèle nominal maintienne certaines performances (comme la bande passante, 
précision, etc.). En absence des incertitudes, l’analyse de la robustesse en performance 
consiste donc à satisfaire une condition de robustesse appelée «condition sur les 
performances nominales NP » [9]. Par contre, si les incertitudes de modélisation sont prises 
en considération [10], l’analyse de robustesse cherche donc à établir si le système bouclé 
maintient les performances prévues pour tous les éventuels régimes perturbés ou non. Dans ce 


cas, on cherche à satisfaire une condition de robustesse appelée «condition sur les 


performances robuste RP » [11]. 


Habituellement, les incertitudes de modélisation ont la capacité de déstabiliser le système 
bouclé et diminuer ses performances. C’est pour cette raison que l’analyse de robustesse 
cherche donc à assurer un bon compromis entre le suivi de la trajectoire de référence et la 
robustesse vis-à-vis le type d’incertitudes affectant le modèle de synthèse. Dans ce cas, on 
cherche à satisfaire une condition de robustesse appelée « compromis RS/NP on en anglais 


Trade-off » [12,13]. 


Dans ce travail, il est possible d’aller encore plus loin lors de la synthèse du contrôleur 
robuste. Nous allons chercher à satisfaire la condition sur RP qui englobe les trois conditions 
de robustesse sur RS, NP et sur RS/NP. Le problème de la commande H,, est de comment 
assurer la condition sur RP du système bouclé avec une marge de robustesse très élevée? 
Nous verrons que la solution fournie par la commande H., basée sur des améliorations 
proposées présentera l’objet de ce travail. Le paragraphe suivant représente un bref aperçu 
historique sur l’évolution de la commande H,, depuis ses principes de base « bien connus » 
jusqu’aux développements les plus récents en terme de la recherche des pondérations 


adéquates utilisées dans leur problème de synthèse. 
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La commande H,, représente un moyen relativement simple d’aborder une analyse et 
synthèse du contrôleur dans le domaine fréquentiel. Elle a été introduite, en premier lieu, par 
le mathématicien Doyle en 1982 [2] pour commander des systèmes mono-variables SISO 
dont les principales règles de réglage qui permettent l’analyse du système bouclé dans le plan 
fréquentiel ont été suggérées par d’autres chercheurs. Ensuite, en 1983 [14], les deux 
automaticiens Zames et Francis ont résolu le problème de la commande des systèmes 
SISO basés sur des modèles nominaux simplifiés qui ne prennent pas en compte toutes les 
complexités des systèmes physiques. Enfin, en 1984, Francis , Helton et Zames [15] ont 
abordé à l’extension de l’application de la commande H., sur des systèmes MIMO. En 1986, 
les trois chercheurs Vidyasagar et Kimura et Glover [16] ont synthétisé un contrôleur 
robuste MIMO avec lequel des spécifications fréquentielles ont été satisfaites par l’application 
des algorithmes très performants, qui sont aujourd’hui disponibles dans la bibliothèque 


du Matlab. 


En général, la synthèse du contrôleur robuste H,, nécessite, a priori, la minimisation de la 
norme H, d’un critère présenté par transfert complexe. Ce dernier reliant plusieurs entrées 
exogènes, à savoir : les perturbations, bruits, références,.…, etc., avec des sorties exogènes on 
l’occurrence: les erreurs pondérées, commandes pondérées,..., etc. Ce critère peut être 
formulé via l’application de la transformation linéaire fractionnaire LFT [17,18] (ou en 
anglais Linear Fractional Transformation) inférieure sur un système d’interconnexion 
composé par deux transferts complexes à savoir: celui du modèle généralisé (modèle 
augmenté) et celui du contrôleur robuste à synthétiser. Habituellement, le modèle généralisé 
est composé du modèle nominal et certain nombre de spécifications fréquentielles qui sont 
souvent imposées par un tel cahier de charges. Parmi les spécifications les plus considérées 
lors de la synthèse par la commande H.,., on peut citer : spécification sur RS, spécification 
sur NP, spécification sur la limite des signaux de commandes....etc. Le pas décisif, pour 
minimiser ce critère, réside dans le choix d’un algorithme adéquat avec lequel la 


commande H, à la capacité de satisfaire toutes les performances du cahier de charges. 


Habituellement, La commande H,, fait un appel à l’un des deux algorithmes 
d'optimisation, qui sont disponibles dans la bibliothèque du Matlab. On aura donc la 
commande H., utilisant les deux équations algébriques de Riccati (ARESs) [19, 20, 21,22,23] 
et la commande H,, utilisant les Inégalités matricielles linéaires (LMISs) [24, 25, 26, 27, 28]. 
Une fois le contrôleur robuste est développé, une analyse de robustesse peut être effectuée via 


le tracé d’un certain nombre des lieux des valeurs singulières maximales et minimales [29, 30, 
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31] de quelques sensibilités fournies par le système bouclé. Parmi ces sensibilités les plus 


souvent testées, on peut citer: sensibilité directe, sensibilité complémentaire, ...etc [32, 33]. 


Par conséquent, Ces sensibilités permettant d’analyser l’impact de différentes entrées 
exogènes sur le système bouclé. Par exemple, le tracé des lieux des valeurs singulières 
maximales de la sensibilité complémentaire permet de découvrir à quelle marge de robustesse 
de la stabilité robuste est assurée en présence des erreurs de modélisation de type 
mulüplicative. En effet, une marge de robustesse faible de cette condition peut se traduire par 
la dégradation de la stabilité robuste ce qui est traduit [33], dans le plan temporel, par l’impact 
important de l’effet des bruits de mesure ainsi que l’effet des modes négligée en hautes 
fréquences sur la boucle de commande. De plus, le tracé des valeurs singulières de la 
sensibilité directe permet de découvrir les informations qui concernent les performances 
nominales [33]. D’où, une faible marge de robustesse peut être se traduit, dans le plan 
temporel, par la dégradation de la dynamique de suivi de référence ainsi que l’abaissement 


lent de l’entrée des perturbations. 


Toutefois, pour certaines applications industrielles, la commande H, basée sur les 
algorithmes précédents présente l’avantage de fournir un bon suivi de la trajectoire de 
référence dont la réponse temporelle est souvent caractérisée par un dépassement faible, un 
temps de montée rapide et une erreur statique presque nulle en régime permanant. Elle 
présente aussi l’avantage de fournir également des bonnes dynamiques de rejection des 
perturbations et atténuation de l’effet des bruits de mesure. La question posée réside dans 
l'assuré d'une bonne marge de sécurité du compromis ” Poursuite/Régulation" en présence 


des incertitudes de modélisation?) 


Selon des travaux antérieurs, la réalisation d'un bon compromis entre NP et RS dépend 
largement du choix appropriée des pondérations fixes d’ordre entier [34] avec lesquelles la 
formulation du problème de sensibilité pondérée peut prendre en compte la satisfaction de 
toutes les spécifications imposées par le cahier de charges. Les paramètres de chaque 
pondération doivent être choisis préalablement par l’utilisateur où des règles de réglages 
existant dans la littérature doivent être respectées. Le rôle de ces pondérations réside dans la 
limitation de l’évolution de chaque sensibilité en termes des valeurs singulières maximales. 
Cette limitation permet de ressortir un autre problème de commande appelé « problème 
sensibilité pondérée », ou en anglais «weighted sensitivity problem». La solution 
optimale de ce problème par la méthode H,, peut offrir la représentation d’état du contrôleur 


robuste, appelé par la suite « contrôleur H,, primaire ». 
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Toutefois, pour certaines exigences industrielles, la satisfaction de la condition sur 
la RP est indispensable [11, 35] lors de la synthèse du contrôleur robuste. Dans ce cas, le 
critère fourni par cette condition de robustesse est souvent très difficile à minimiser. C’est 
pour cette raison, qu’on doit chercher à d’autres types des pondérations, car la commande H,, 
peut fournir une solution vide ce qui produit plusieurs anomalies, telles que: L’obtention 
d’une commande fluctuante sur un intervalle du temps très courtes, ce qui pouvant 
déstabiliser le système bouclé, l’obtention d’une marge réduite du compromis de robustesse 
réduisant la sécurité contre les incertitudes de modélisation, l'effet indésirable des bruits de 


mesure et les dynamiques négligées en hautes fréquences. 


Par conséquent, le choix manuel des paramètres de ces pondérations est souvent difficile. 
C’est pour cette raison qu’on doit introduire un algorithme d’optimisation sous contrainte 
assurant le choix automatique de chaque pondération intervenant dans la construction du 
modèle généralisé. À cet effet, la plupart de recherches avancées ont été orientées vers 
l’utilisation des pondérations irrationnelle ajustable dans le problème de sensibilité pondérée 
où les paramètres de chaque pondération sont d'abord optimisés par l’algorithme 
d'optimisation  évolutionnaire, appelé «algorithme génétique» GA (où en 
anglais Genetic Algorithm ). Ensuite, la représentation d'état du contrôleur robuste est 
assurée par l’application de la commande H,, sur le système d’interconnexion composé par le 
modèle nominal et les pondérations irrationnelles ajustables fournies par GA. Ces deux étapes 


précédentes sont réitérées jusqu'à ce que la condition RP soit satisfaite. 


Habituellement, la sélection automatique des pondérations irrationnelles ajustables est 
assurée via la satisfaction de quelques règles existant dans la littérature en utilisant 
l’algorithme génétique. A cet effet, plusieurs chercheurs ont suggéré des formes idéales pour 
chaque pondération irrationnelle [36, 37] ajustable intervenant dans la formulation du modèle 
généralisé. Toutefois, la commande H., basée sur ce type des pondérations à la capacité de 
satisfaire la condition sur RP. Néanmoins, pour un nombre élevé des spécifications 
fréquentielles, ce type des pondérations peut restreindre l’application de la commande H,, sur 
certains processus industriels. C’est pour cette raison, qu’on doit chercher à d’autres types des 


pondérations plus générales que celles citées précédemment. 


Notre objectif étant plus précis pour pallier à ces anomalies. Nous allons introduire, dans le 
problème de sensibilité pondérée précédent, des pondérations fractionnaires ajustables dont 
ses paramètres sont optimisés systématiquement par un algorithme d’optimisation appelé 


« algorithme d’essaim particulaire PSO (ou en anglais Particle Swarm Optimization » [38, 
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391. le processus d’optimisation doit satisfaire les mêmes réglages proposées pour les 
pondérations irrationnelles ajustables citées ci-dessus. Notons ici que l'avantage d’imposer 
des pondérations rationnelles ajustable dans la synthèse du contrôleur robuste réside dans le 
fait qu'elles peuvent assurées des atténuations fractales dans les lieux des valeurs singulières 
maximales des sensibilité du système bouclé. Ces dernières jouent un rôle important dans la 


plage de fréquence ou le compromis de robustesse est difficile à atteindre. 


En général, un transfert fractionnaire est souvent rapproché par une série des transferts 
usuels d'ordre entier ou chacun est présenté par un nombre infini de pôles et de zéros. A cet 
effet, l’algorithme PSO peut bénéficier de cette propriété importante, car la dimension du 
problème d’optimisation assurant la sélection automatique de ces pondérations est réduite de 
façon très remarquable. Ceci représente, plus tard, une importance capitale pour l’accélération 
de la convergence du PSO vers des solutions optimales. Il est à noter que la méthode H,, 
nécessite l’utilisation des pondérations d’ordre entier. Cela peut poser des problèmes sérieux, 
car l’utilisation des pondérations irrationnelles est inévitable pour améliorer la marge de 
robustesse RP. À cet effet, nous allons rapprocher [40] ces pondérations par des transferts 


irrationnels ce qui permet l’application de la méthode H... 


Récemment, plusieurs méthodes d’approximation ont été proposées pour approximer les 
transferts d’ordre non entier. Parmi ces méthodes d’approximation les plus utilisées dans 
l’automatique avancée, on peut citer : les approximations d’Oustaloup [41], Matsuda [42] 
et Charef [43] ainsi que les méthodes utilisant l’identification fréquentielle ou temporelle 


[44]. 


Notons ici que l’amélioration requise de la marge de robustesse de RP permettant de 
robustifier toutes les classes des contrôleurs H,, primaires dont leurs paramètres sont obtenus 
par la solution du problème de sensibilité pondérée basée sur les pondérations rationnelles 
ajustables. La méthode de robustification proposée, dans ce travail, sera basée sur les quatre 
étapes suivantes : dans la première étape, un contrôleur H, primaire est synthétisé via la 
résolution du problème de sensibilité pondérée formulée via des pondérations irrationnelles 


fixes. 


Dans la deuxième étape, le problème précédent est reformulé par d’autres pondérations 
rationnelles ajustables inconnues. Ensuite, à partir d’un ensemble des règles de réglage, les 
paramètres de ces pondérations seront optimisés par l’algorithme du PSO. Dans la troisième 
étape, la réponse fréquentielle de chaque pondération rationnelle PSO sera calculée dans une 


bande de fréquence où l’approximation est nécessaire. Cette réponse permettant de fournir le 
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filtre d’ordre entier rapprochant le module et la phase du transfert rationnel. Les paramètres de 


ce filtre seront déterminés en utilisant l’approche basée sur l’identification fréquentielle. 


Dans la dernière étape, le modèle généralisé sera développé à partir du modèle nominal 
ainsi que les filtres d’ordre entier fournis par l’étape précédente. Le problème de sensibilité 
pondérée correspondant est ensuite reformulé et la méthode H,, sera réutilisée pour fournir la 
représentation d’état du nouveau contrôleur. La mise à jour de la marge de robustesse 
de RP est assurée par ce dernier contrôleur et les quatre étapes ci-dessus seront répétées en 
fonction de la décroissance de cette marge de robustesse. Ce processus de robustification 
proposé sera terminé par la synthèse d’un contrôleur H, robustifié assurant la condition RP 


avec une marge de sécurité très élevée. 


Le contrôleur H,, primaire et celui de H,, robustifié seront appliqués sur deux procédés 
réels dont ses dynamiques ont été modélisées par des modèles linéaires incertains et ses 
incertitudes ont été évaluées par la forme multiplicative en sortie. Les résultats de simulation 
obtenus seront comparés dans les plans fréquentiel et temporel. Ces résultats obtenus sont très 
innovants chose qui nous a permis de les publié dans des revues scientifiques de renommée 


internationale. 
Ce travail est organisé selon le modèle suivant : 


Le premier chapitre sera subdivisé en trois parties. Dans la première partie nous allons 
présenter les outils mathématiques fondamentaux de la commande robuste d’ordre entier tels 
que: la représentation des systèmes physiques modélisés par des modèles incertains, les deux 
transformations linéaires fractionnaires basse et haute, définition de différentes normes ainsi 
que les valeurs singulières d'un transfert complexe et finalement la classification des 
incertitudes de modélisation. Dans la deuxième partie, nous allons présenter le théorème du 
faible gain avec lequel on pourra tirer la condition sur la stabilité robuste pour chaque type 
d’incertitudes de modélisation, la condition sur les performances nominales et les 
performances robuste, l’analyse de robustesse du système bouclé à travers le tracé du lieux 
des valeurs singulières maximales de chaque sensibilité fournie via le bouclage du modèle de 
synthèse avec le contrôleur robuste à synthétiser et finalement la formulation du problème de 
sensibilité pondérée ainsi que les deux algorithmes utilisés dans la commande H, pour 
résoudre le problème de la commande. La troisième partie sera consacrée à l’analyse et 
synthèse du contrôleur H, primaire pour la commande d’un procédé hydraulique mal 
conditionné (ill — conditionned system) simulant le fonctionnement d’un canal à surface 


libre. Le comportement réel de ce système est modélisé par un modèle fractionnaire possédant 
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un retard pur. De plus, des incertitudes de modélisation de type multiplicatives agissant sur le 
modèle de synthèse seront d’abords analysées puis les majorées par un transfert complexe qui 
permettra le développement de la condition sur la stabilité robuste du système bouclé et ce par 
l’application du théorème du faible gain. Ensuite, le modèle généralisé sera développé à partir 
du modèle nominal et quelques pondérations irrationnelles fixes fournies par les spécifications 
du cahier de charges. Finalement le contrôleur robuste sera synthétisé par la commande H,, et 
l’analyse du système corrigé sera abordée dans le plan fréquentiel et temporel. Ce chapitre 
sera terminé par une conclusion dans laquelle le problème de la commande H,, basée sur les 
pondérations irrationnelles fixes sera présenté ainsi que les éventuelles solutions proposées 


faisant l’objet du deuxième chapitre. 


Le deuxième chapitre sera consacré à la synthèse du contrôleur H,, via la résolution du 
problème de sensibilité pondérée basé sur les pondérations irrationnelles ajustables. Il se 
compose par les quatre parties. Dans la première partie, nous allons présenter les nouvelles 
structures fractales proposées pour la conception des pondérations irrationnelles ajustables. 
Ces dernières à pour but de palier le problème du choix manuel des pondérations adéquates 
satisfaisant la condition sur RP. Sachant que les paramètres de ces pondérations peuvent être 
optimisés par l’algorithme génétique tout en satisfaisant un ensemble des règles de réglage 
existant dans la littérature. Pour cela, la deuxième partie sera consacrée à la présentation de 
ces règles ainsi que la formulation du problème de la sélection automatique des pondérations 
irrationnelles sous forme d’un problème d’optimisation sous contraintes. Sachant que la 
solution optimale de ce problème fait l’objet de la troisième partie de ce chapitre. Pour cela 
nous allons donner un bref aperçu sur l’algorithme génétique qui est destiné à la recherche des 
pondérations dont la puissance de ses transferts est d’ordre entier. Finalement, dans la 
dernière partie de ce chapitre, nous allons contrôler le même système hydraulique cité en 
chapitre précédent par la commande Æ,, utilisant les pondérations irrationnelles ajustables. La 
sélection automatique de ces pondérations sera effectuée par l’algorithme génétique et les 
résultats de simulation seront comparés à ceux fournis par le contrôleur H, primaire 
synthétisé en premier chapitre. Ce chapitre sera terminé par une conclusion dans laquelle 
l’anomalie apparue par l’utilisation de ce type de pondération sera résolue dans le chapitre 


suivant. 


Le troisième chapitre à pour but d'améliorer la robustesse du système bouclé en utilisant 
les pondérations rationnelles ajustable lors de la synthèse du contrôleur robustifié. Pour cela, 


ce chapitre sera dédié à exposer les définitions de base utilisées dans le calcul fractionnaire, 
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en particulier les définitions de l’opérateur intégro-différentiel d’ordre non entier, les 
différentes méthodes d’approximation de l’opérateur intégro-différentiel d’ordre non entier à 
savoir les méthodes de Oustaloup, Charef, ..….. etc., et la méthode utilisant l’identification 
fréquentielle. Sachant que la structure fractale proposée pour chaque pondération comportant 
moins de paramètres inconnus. Ceci peut réduire la dimension du problème d’optimisation 
garantissant la sélection automatique des pondérations rationnelles désirées. Cette réduction 
peut améliorer la convergence de l’algorithme du PSO à des solutions optimales avec 
lesquelles on assurera le choix des pondérations rationnelles souhaités. La troisième partie de 
ce chapitre sera consacrée à la validation de cette robustification proposée sur le même 
système hydraulique précédent et les résultats obtenus par le contrôleur H,, robustifié seront 
comparés avec ceux fournies le contrôleur H,, synthétisé à partir des pondérations 


irrationnelles ajustables présentées en deuxième chapitre. 


Le quatrième chapitre est dédié à la commande robuste d’un moteur synchrone à aimant 
permanant PMSM (Permanant Magnet Synchronous Moto) où sa dynamique réelle a été 
linéairisée autour d’un point de fonctionnement au moyen du logiciel Simulink de Matlab. 
Ce dernier produit un modèle MIMO de synthèse caractérisé par le découplage des signaux 
d’entrées avec ceux de sorties. Pour la synthèse du contrôleur H, robuste (primaire et 
robustifié), nous allons appliquer la stratégie basée sur la commande du flux orienté FOC 
(Oriented Flux Control) pour découpler le système à commander en deux sous 
systèmes SISO. La vitesse mécanique du PMSM est donc contrôlée par les deux contrôleurs 
robustes cités ci-dessus est les résultats obtenus seront comparés, dans le plan fréquentiel et 
temporel, non seulement pour le fonctionnement nominal du moteur, mais également pour un 
ensemble des régimes perturbés. Finalement, ce travail sera achevé avec une conclusion 
générale et quelques perspectives proposées pour les futurs candidats intéressés par la 


robustification d’une classe des contrôleurs H,, primaires. 
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CHAPITRE Î : 
COMMANDE ROBUSTE UTILISANT 


LES PONDÉRATIONS IRRATIONNELLES 


FIXES 


CH.1: COMMANDE ROBUSTE UTILISANT LES PONDÉRATIONS IRRATIONNELLES FIXES 


Chapitre 1 


Commande robuste utilisant les pondérations 
irrationnelles fixes 


1. Introduction : 

Ce chapitre porte les notions de base de la synthèse des contrôleurs robustes pour 
commander des processus industriels présentés par des modèles incertains. A cet effet, la 
méthode de synthèse qui permet la synthèse du contrôleur stabilisant la boucle de commande est 
celle de H, avec laquelle des conditions de robustesses doivent être assurées et cela non 
seulement pour le régime de fonctionnement nominal du système physique à commander, mais 
également pour tous les régimes perturbés générés par les incertitudes de modélisation, les 
dynamiques non linéaires négligées en haute fréquence ainsi que la présence les effets des bruits 


de mesures provoqués par la fatigue les capteurs installés dans la boucle de commande. 


En général, la synthèse du contrôleur robuste par la méthode H,, s’articulera sur les étapes 
de base suivantes : Tout d’abord, les incertitudes de modélisation affectant le modèle de synthèse 
doivent les classer sous forme des incertitudes structurées et/ou non structurées. Ensuite, elles 
sont évaluées à partir de tracer de ses lieux des valeurs singulières maximales (1.e., ses modules) 
dans le plan fréquentiel où un transfert stable majorant toutes ses incertitudes doit être choisi et à 
partir duquel le théorème des fiables gains est appliqué pour fournir la condition de robustesse 
sur la stabilité robuste. Cette dernière doit être accompagnée par d’autres conditions de 
robustesse à savoir celle sur les performances nominales et celle sur la limitation des signaux de 


commande qui assurent le compromis entre la poursuite et la régulation du système bouclé. 


La satisfaction de ces conditions peut être assurée par un choix a priori des fonctions de 
pondérations décrivant toutes les spécifications imposées par le cahier de charge. Ces 
pondérations sont utilisées pour augmenter le modèle nominal, ce qui fournit le modèle 
généralisé. L’application de la transformation linéaire fractionnaire inférieure sur ce modèle 
permettant d’élaborer le critère de sensibilité mixte généralisée où sa solution optimale peut être 
assurée par la méthode H,, utilisant les équations algébriques de Riccati ou la méthode H,, 


utilisant les inégalités matricielles linéaires (LMIS). 


Ce chapitre sera achevé par une partie de simulation dans laquelle la dynamique réelle du 
processus à commander est décrite par un modèle incertain. Le but principal de cette partie de 


simulation est de découvrir les anomalies ainsi que les difficultés d’application de cette stratégie 
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de commande où la solution de ces anomalies présentera l’objet de la partie suivante de ce 


travail. 
1.1. Représentation d’un système linéaire invariant dans le temps: 


La dynamique d’un processus industriel comportant p grandeurs d’entrée et q grandeurs de 


sortie est souvent décrite par la représentation d’état ci-dessous [45,46]: 


ie = A:x(t)+B'u(t) (1.1) 
y(E) = C-x(t) + D -u(t) | 


x(t) € R” : représente le vecteur d’état; 


u(t) € R? : représente le vecteur d’entrées du système à commander (grandeurs réglant); 


y(t) € R‘ : représente le vecteur de sorties du système (grandeurs à régler) ; 


e AE R"*”": la matrice d'évolution du système ; 
e BE R"”“?: la matrice de commande du système ; 
e CE R?*": la matrice de sortie (observation) du système ; 


e DE RP*1: la matrice de transmission directe du système. 


La matrice de transfert du système G(s) € CP*4 est obtenue par l’application de l’expression 


(1.2) comme suit: 
G(S) = C(s.lyxn — A) !B + D (1.2) 


Notant ici que l’existence de G(s) nécessite la non singularité de la matrice de 
transfert (5. lixn — À) € C**. De ce fait, la réalisation minimale de la matrice d’évolution A 
est souvent nécessaire ou le couple (4, B) doit être contrôlable ainsi que le couple (4, C) doit 


être observable. 


1.1.1. Valeurs singulières et normes d’un transfert: 


Les valeurs singulières d’une matrice de transfert, notée ici par oi(G (jw)), sont données par 


[471]: 


G(GGa)) = JA(GGw) : GC-jo)T) (13) 


Où 1; sont les valeurs propres du transfert [G(jw) : G(—-jæ)"]. Sachant que oi(G (jw)) sont 


toujours positives, il est possible de les classer par ordre croissant telle que : Ghin(G) est la plus 
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petite valeur singulière et ohn4x(G) est la plus grande valeur singulière. On obtient donc les 


inégalités ci-dessous : 
Omax(G) = 01(G) Z 02(G) Z 03(G) 2: Z Omin(G) (1.4) 


1.1.2. Normes d’un transfert : 


La norme est un opérateur qui permet de quantifier une fonction où une matrice de transfert. 
Une norme du transfert G(s) est généralement notée par |[G(jw)||; où € d’esigne la nature de la 
norme considérée. De ce fait, plusieurs normes différentes peuvent donc être utilisées dans la 


théorie de la commande robuste. Parmi lesquelles, on peut citer : La norme H, et La norme 


H>. La norme H, du transfert G(s) est, en générale, définie par [47]: 


L 


; Y D6 Ce) Gps ; 2 
1GGa)le = Supucoen: (ce) = Gris eracel6 (a) Gw)du) 


(1.5) 


Elle peut également s’exprimer, en fonction des valeurs singulières, par la relation suivante : 


1 
. IIY(S)Ilco 1 p+o : 2 
IGGw)l2 = Supuces, (QE) = GIE 5 (GGw)dw 


(1.6) 


La norme H,, du transfert G(s) représente le gain, noté y, qui correspond le pic de résonance 


de la réponse fréquentielle G(jw). Ce gain est souvent donné par l’expression ci-dessous: 


IYCS)1l2 
IU(s)Il2 


7 = IIGGo)Ile = Supuçyen, (RTE) = Max 0maxlG G&)] 


(LD 


1.2. Interconnexion d’un système bouclé : 


Dans la synthèse du contrôleur robuste, on désigne sous le nom de représentation LFT ( 
Linear Fractional Transformation) un système pouvant s’écrire sous la forme d’une matrice 
d’interconnexion, notée P(s). Cette dernière peut être connectée en haut par un transfert, 
noté A,(s) modélisant toutes les incertitudes qui affectent le modèle de synthèse. Elle peut être 
également connectée en bas par un transfert, noté K(s), décrivant le contrôleur à synthétiser. En 
général, suivant le bouclage de P(s), on peut distinguer deux représentation de LFT [18, 46], 


[48] à savoir: LEFT basse et LFT haute. 


1.2.1. Transformation linéaire fractionnaire basse LFT , : 
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Soit P(s) le modèle généralisé comportant toutes les spécifications du cahier charge et le 
modèle nominal de synthèse [47, 49]. Le système d’interconnexion, illustré par la Fig.1.1, est 
souvent adopté pour formuler le problème de la commande robuste. Cette dernière est aboutie 
via l’application de LFT, sur le système d’interconnexion reliant P(s) par le transfert à 
déterminer K(s). Par conséquent, on obtient le transfert T,,,,,(s) qui représente la relation entre 
le vecteur d’entrée exogène w(s), qui peut contenir plusieurs signaux à savoir : la trajectoire de 
référence, la perturbation amenant en entrée et/ou en sortie du modèle de synthèse, les bruits de 
mesures, .. etc., et le vecteur de sortie exogène z(s), qui peut contenir plusieurs signaux à 
savoir : l’erreur modérée de poursuite, la commande modérée, .… etc. En général, le bouclage en 


bas de P(s) par K(s) est souvent présenté par la Fig.1.1 ci-dessous : 


Fig.I1.1. Configuration standard de LFT,. 


Pi Pi2 
Sachant que la matrice complexe P(s) est partitionnée comme suit P(s) = | MSN } 
Po } P22 


Z w 
telle que (-) = P(s): (-) d’où, la relation reliant la sortie z à l’entrée w peut s’écrire comme 
e u 


la transformation linéaire fractionnaire basse suivante : 


z(s) = {LFT/(P(s), K(S))}.w(s) (1.8) 
Twz(s) 
Avec : 
LETy(P,K) = Py2.(1 — K.P32) 2. K.P31 + Pu (1.9) 


Remarque : L’existence du transfert T,,,,(s) est conditionné par la non singularité da la 


matrice (1 — K.P,2) 1. 
1.2.2. Transformation linéaire fractionnaire haute LEFT, : 


En général, le modèle de synthèse est souvent entaché par des erreurs de modélisation. La 
représentation LFT,, peut servir à représenter un modèle incertain sous la forme du modèle 


nominal bouclé sur une matrice d’incertitudes représentant ces erreurs [47, 49]. Cette 


Page | 13 


CH.1: COMMANDE ROBUSTE UTILISANT LES PONDÉRATIONS IRRATIONNELLES FIXES 


transformation est appliquée sur le système d’interconnexion reliant la matrice de transfert 
complexe P(s) par le transfert d’incertitudes A,(s). En général, le bouclage en haut de P(s) par 


A,(s) est souvent présenté par la Fig. 1.2 ci-dessous : 


Fig.1.2. Configuration standard de LFT, 


De même il est possible de définir le transfert T,,,,,(s) comme suite : 


z(s) = {LFT, (P,A,)}.w(s) (1.10) 


Avec : 
LFT,(P, À,) = Pos. ( — A pe RAP; + P5 (1.11) 


A noter ici que la synthèse du contrôleur robuste basée sur un modèle incertain consiste à 
déterminer le transfert du contrôleur K(s) minimisant la norme H,, de T,,,,(s). Ce dernier est 
obtenu via l’application de LFT, sur le système d’interconnexion reliant la matrice de transfert 
complexe M(s) par le transfert K(s), sachant que la matrice M(s) est obtenue via l’application a 
priori de LFT,, sur le système d’interconnexion reliant la matrice de transfert complexe P(s) par 
le transfert A,(s). En général, la configuration standard complète de (P — K — A,) est souvent 


présentée par la Fig.1.3 ci-dessous : 


Fig.1.3. Configuration standard (P — K — À.) 


Afin de pouvoir aisément séparer la conception de l’analyse de robustesse, il est habituel de 


représenter le modèle incertain sous une forme telle que les incertitudes soient extraites du 
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modèle nominal. Le paragraphe suivant donne les différentes structures possibles pour ces 


incertitudes et leurs natures. 


1.3. Les incertitudes de modélisation: 


ES 


Il existe deux types d’incertitudes à savoir : incertitudes non structurées et incertitudes 


structurées. 
1.3.1. Incertitudes non-structurées : 


Ce type d'incertitude se produit souvent à haute fréquence. Ceci est dû aux deux facteurs 
principaux [46]. Le premier est celui causé par les signaux de perturbations agissant sur la 
commande et la sortie du système bouclé. Le deuxième facteur est celui causé par les modes 
négligés lors de la modélisation du système réel, en particulier l’étape de la linéarisation autour 
d’un point de fonctionnement, la négligence des retards purs dans la boucle de commande. Il 
convient de noter que l'extraction des incertitudes à partir des modèles unidimensionnels 
incertains devient assez facile par rapport à celle à partir des modèles multidimensionnels 
incertains, puisque toutes les opérations mathématiques appliquées sur les matrices de transfert 
deviennent sous forme des produits, sommes et divisions appliquées sur des fonctions de 
transfert. Pour cela, le paragraphe suivant s’articulera sur les différents types d’incertitudes 
agissant sur les systèmes multivariables. A cet effet, les nomenclatures ci-dessous seront utilisées 


dans la suite de ce chapitre. On aura donc : 


e Gn(s) : est la matrice de transfert du modèle nominal ; 
e_G,(s) : est la matrice de transfert du modèle perturbé ; 


e_A,(s) : est une matrice d’incertitudes agissant sur le modèle de synthèse. Elle représente 
l’écart absolu ou relatif entre le fonctionnement nominal et perturbé du système réel. 


1.3.1.1.Incertitude additive : 


Ce type d’incertitudes est bien adapter pour la modélisation des dynamiques non linéaires [49, 


50]. En général, il existe deux représentations de ce type d’incertitudes à savoir : 
> Incertitude additive directe : 


Pour cette forme d’incertitudes la famille de modèles G, se définit comme la LFT, de la 


matrice complexe M bouclée par la matrice d’incertitude additive A, (voir la Fig.1.4). Sachant 


: VA : 0 I Aw . 0 ll  … 
que Aw = À, (s):Az et ( = = È | ( : ) , avec M = Ë cl) d’où, le transfert du 
modèle perturbé est donné par l’expression ci-dessous : 
Gps) = LFT,(M, Aa) (1.12) 
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(a) : transfert perturbé (b) : G, donné sous forme LTF, 
Fig. 1.4. Incertitude additive directe 


De l’équation (1.12), on obtient la relation suivante : 


Gps) = Gn(s) + AQ(s) (1.13) 
> Incertitudes additive inverse : 


Pour la forme additive inverse, la famille de modèles est présentée par la Fig. 1.5 ci-dessous : 


Fig. 1.5. Incertitude additive inverse 


A partir de la Fig.1.5, l’application de LFT,, permet de fournir l’expression ci-dessous : 


—Gn Gn 
GS) = LFTi(M, Aa), avec M = E G, 2 (1.14) 


D'où, la matrice de transfert du modèle perturbé, en présence des incertitudes additive 


inverse, est obtenue par la relation ci-dessous : 
-1 
G(S) = [Gy*(s) + Aa,(s)] (1.15) 


1.3.1.2. Incertitude multiplicative : 


Ces formes sont utilisées pour prendre en compte les erreurs de modélisation des actionneurs 
et des capteurs ainsi que les dynamiques négligées [50]. Il existe quatre classes de formes 


mulüplicatives: 
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> Incertitudes multiplicative en sortie : 


Pour la forme multiplicative en sortie, la famille de modèles est présentée par la Fig.1.6 ci- 


dessous : 


Fig. 1.6. Incertitudes multiplicative en sortie 


A partir de la Fig.1.6, l’application de LFT,, permet de fournir l’expression ci-dessous : 


Ü G 
G(s) = LFT, (M, A), avec M = | ï (1.16) 
p F CG 


D'où, la matrice de transfert du modèle perturbé, en présence des incertitudes multiplicative 


en sortie, est donnée par l’expression ci-dessous: 
G(S) = [1 + A;(s)]. Gy(s) (1.17) 
> Incertitudes multiplicative inverse en sortie : 


Ce type d’incertitudes est présenté par la Fig.1.7 ci-dessous : 


Fig. 1.7. Incertitudes multiplicative inverse en sortie 


A partir de la Fig.1.7, l’application de LFT,, permet de fournir l’expression ci-dessous : 


-I G 
G(S) = LFT,(M,A,,), avec M = ce | (1.18) 


D'où, la matrice de transfert du modèle perturbé, en présence des incertitudes multiplicative 


inverse en sortie, est donnée par l’expression ci-dessous: 


GS) = [1+A,()] .Gwts) (1.19) 
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> Incertitudes multiplicative en entrée : 


De la même façon il est aussi possible de définir la forme multiplicative en entrée. Ceci peut 


être illustré par la Fig.1.8 ci-dessous: 


Fig. 1.8. Incertitudes multiplicative en entrée 


A partir de la Fig.1.8, l’application de LFT,, permet de fournir l’expression ci-dessous : 


G,(s) = LFT,(M, AS), avec M = | . | (1.20) 
F Gv  Gn 


D'où, la matrice de transfert du modèle perturbé, en présence des incertitudes multiplicative 


en entrée, est donnée comme suit: 
Gps) = Gn(s). [I + A,(s)] (1.21) 


> Incertitudes multiplicative inverse en entrée : 


De la même façon il est aussi possible de définir la forme multiplicative inverse en entrée. 


Ceci peut être illustré par la Fig.1.9 ci-dessous: 


Fig. 1.9. Incertitudes multiplicative inverse en entrée 


A partir de la Fig.1.9, l’application de LFT,, permet de fournir l’expression ci-dessous : 


1 1 
G(s) = LFT (M, A), avec M = É Gy cl (1.22) 


D'où, la matrice de transfert du modèle perturbé, en présence des incertitudes multiplicative 


inverse en entrée, est obtenue comme suit : 
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—1 
G(S) = Gn(s). [1 + A, (s)] (1.23) 
1.3.2. Incertitudes structurées: 


Dans l’industrie, l’incertitude d’un élément physique peut être apparue à tous moment 
pendant le fonctionnement du système réel à commander [2]. Cette dernière est due au 
changement du point de fonctionnement et/ou la présence des phénomènes physiques affectant la 
boucle de régulation. Tous ces cas peuvent être regroupés de façon à ce que la matrice 
d’incertitudes soit structurée de telle sorte qu’elle soit diagonale par blocs, chaque bloc 
correspondant à une formulation différente des incertitudes. Dans notre travail, nous ne nous 
intéresserons que par le type d’incertitudes non structurées, la détermination de la matrice de 
transfert décrivant les incertitudes structurées pourra faire l’objet à d’autres recherches 


ultérieures. 


1.4. Analyse de robustesse des systèmes dynamique : 

1.4.1. Stabilité Nominale (NS) : 

La stabilité nominale est indispensable pour chaque système de commande. Elle permet de 
façonner la réponse du système à une référence donnée et garder le fonctionnement nominal du 
système physique stable en boucle fermée. En général, on distingue deux définitions pour la 
stabilité nominale : 

e La stabilité externe exige que pour toute entrée bornée injectée au système bouclé, on 
obtient une sortie bornée. Cette stabilité est encore appelée stabilité au sens EBSB (en 
anglais Bounded Input Bounded Output B1BO). 

e La stabilité interne est plus large que celle mentionnée précédemment. Elle exige que 
tout signal injecté en n’importe quel point de la boucle génère une réponse bornée en 
tout autre point. 


1.4.2. Stabilité Robuste (SR) : 


La stabilité est dite robuste si, en plus du modèle nominal, la stabilité doit être garantie pour 
toutes les perturbations sur le modèle nominal du système, jusqu'au pire cas des incertitudes du 
modèle. La condition sur la stabilité robuste pour un type spécifique de perturbations affectant le 
modèle nominal de synthèse est obtenue par l’application d’un théorème appelé « théorème du 


faible gain ou de petit gain » aux diverses représentations des incertitudes citées précédemment. 


1.4.2.1.Théorème du faible gain : 
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Considérons À, un transfert décrivant les incertitudes de modélisation, et M\(s) la matrice de 
transfert du système entre l’entrée et la sortie de À, qui est supposé stable. Le schéma général, 


qui permet l’étude de la stabilité robuste, peut être illustré par la Fig.1.10 ci-dessous [511] : 


Fig.1.10. Schéma bloc permettant l'analyse de la stabilité robuste d’un système bouclé 


MA(s) est le transfert qui apparait lorsqu’on isole dans un bloc A, toutes les incertitudes de 
modèle. Alors, le système présenté par la Fig. 1.10 demeure stable en présence de la perturbation 
A,si seulement si |A, ||, < 1, d’où la condition sur la stabilité robuste est assurée si et 


seulement si : 


I Mas).4, (5) < 1 & GnaxlMaÜw).A,(jw)] < 1, Vw € R** 
(1.24) 


Une autre condition de la stabilité robuste également suffisante mais plus conservative est la 


suivante : 


I Ma(S)Ilco. 1] 4x (5) < 1, Vo € R* 
(1.25) 


1.4.2.2. Conditions sur la stabilité robuste du système bouclé : 


Pour déterminer la condition sur la stabilité robuste en présence d’incertitudes de 
modélisation, nous devons appliquer le théorème du faible gain sur les différentes configurations 
des incertitudes citées précédemment. Dans ces condition, 1l faut calculer la sortie Aw de A, en 
fonction de sa sortie Az pour obtenir le transfert MA de chaque configuration. Finalement, la 


condition (1.25) permet de fournir les conditions suivantes: 


e Incertitude additive directe | K. Far + Gn.K l” | < ACTE 
L 


e _Incertitude additive inverse | G\:. Hé D K.Gx] | < Ta oL 
Le] ai ” 


1 


RE 1 
e _Incertitude multiplicative directe en sortie | Gn.K. (xs + Gn.K ] | oO 
[ee] S ©O 
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1 


. . . . Li 1 We 
e Incertitude multiplicative inverse en sortie | Le + Gn.K ] | = ] A )| ’ 
co s;\S a 


LE —_— —; 
co I Ae(s)llo ” 
Se 
TO 


: nt : : ; —1 
e Incertitude multiplicative directe en entrée | K.G\. fs a +K.G x] | 


. EE : -1 
e _Incertitude multiplicative inverse en entrée | Ho + K.G »| | 


1.4.3. Performance Nominale (WP) : 


Soit le système bouclé présenté par la Fig.1.11 comme suit : 


Fig.1.11. Structure d'un système bouclé pour l'analyse de robustesse 


Les différents signaux utilisés dans le système bouclé sont les suivants : 
eu : représente le signal de commande fournie par le contrôleur; 
e y: représente le signal de sortie ; 


+ d,, : représente le signal de bruit sur la sortie y; 


d,, : représente le signal de bruit sur la sortie u ; 

en, : représente le bruit de mesure ; 

e r: représente le signal de référence ; 

e e: représente le vecteur des erreurs de poursuite entre les sorties mesurées y et les 
références r. 


Notons ici que ces signaux sont associés entre eux par le système d’équation ci-dessous : 


e = (Lxp + GneK) (r—d,—n,)— Gr (lux + K.Gn) du 
u=K.(Exp+GnK) .(r—d,-n,) —K.Gy.(lyxa + KG) du (1.26) 
y = GnK.(Lxp + GK) (r=n,)+ (Ex + GnK) dd, + Gn(luxa + K.Gne) du 


A partir du système d’équations (1.26), il est claire de remarquer qu’il existe quatre transferts 
qui se répètent dans les expressions de l’erreur, la commande et la sortie du système. Ces 


transferts appelées par la suite « sensibilités ». Elles sont définies par [12] : 


(beton) (1.27) 
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—1 
Su = (loxa + K.Gn) (1.28) 
D RC GR) (1.29) 
L=erGCtRo) (1.30) 


Pour comprendre le rôle des ces sensibilités dans l’analyse de robustesse du système bouclé, 


en essayant de les introduire dans le système d’équations précédent. On aura: 


e = Sal —d, —n,) — Gun. Su dy 
i=KS, (r-d,-7)-L:d, (1.31) 
y = La — ny) +S,.d, + Gn.Su' du 


Lors de l’analyse de robustesse du système bouclé, il est important d’assurer un niveau de 
performances nominales. Ces dernières sont caractérisées souvent par la bonne dynamique de 
poursuite de la trajectoire de référence, la rejection rapide de la perturbation, l’atténuation de 
l’effet des bruits de mesure, la désensibilisation du système bouclé aux effets des modes négligés 
lors de la linéarisation de la dynamique réelle du système physique autour de son point de 
fonctionnement. 

Pour assurer les bonnes performances nominales du système bouclé, il est claire que les 
valeurs singulières maximales de S, , S,,, T, et T,, doivent être minimisées dans toute la plage de 


fréquences. Ces sensibilités sont définies comme suit: 
1.4.3.1.Sensibilité directe en sortie S.,(s): 


Cette sensibilité fournie des informations sur deux dynamiques différentes, à savoir: le rejet 
de la perturbation d., sur la sortie du processus (1.e., S,(s) —= Ta, -y) ainsi que la poursuite de la 
trajectoire de référence r (i.e., S,(s) = T;..). Alors, la première performance qu'on doit 
demander à la régulation est d'assurer un bon rejet de cette perturbation. Autrement dit, on 
souhaite que d,, ait une faible influence sur la sortie y. De plus, la deuxième performance qu’on 
doit demander à l’asservissement est d’assurer une bonne minimisation de l’erreur de poursuite. 
Ces deux performances peuvent être vérifiées à condition que les valeurs singulières de la 
sensibilité directe soient suffisamment faibles, aux fréquences où la perturbation et l’erreur de 


poursuite sont importantes. 
1.4.3.2. Sensibilité directe en entrée S,,(s): 


Cette sensibilité présente l’effet de la perturbation en entrée d, sur la sortie du processus 


(ie, Gn°Suts) = Ta, y). À cet effet, la performance qu'on doit exiger à la régulation est 
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d'atténuer l’effet de cette perturbation. Cela est vérifié à condition que le transfert Gy * S,(s) soit 


suffisament faible où la perturbation d,, est notable. 
1.4.3.3. Sensibilité complémentaire en sortie T.,(s): 


Cette sensibilité serve à nous renseigner sur la stabilité robuste de la boucle fermée vis-à-vis 
les incertitudes multiplicatives en sortie qui affectent le modèle nominal. Elle serve également à 
nous donner des informations sur l’effet de bruits et les modes négligés lors de la linéarisation du 
comportement non linéaire du système. La performance qu'on doit demander à la régulation est 
d'assurer une bonne minimisation de 7,,. Cela est vérifié à condition que les valeurs singulières 
de la sensibilité complémentaire soient suffisamment faibles, aux fréquences où les effets 
susmentionnés sont importants. D’après (1.31), on constate qu’un système bouclé ayant d’une 
bonne poursuite de r assurée avec une atténuation acceptable de d., est tel que S,, = 0. De plus, 
la bonne minimisation de l’effet des bruits de mesure impose que T,, = 0. Ceci n’est pas 


réalisable dans la même plage de fréquence à cause de l’identité ci-dessous : 
T,(s) + S,(s) = 1 (1.30) 


Par conséquent, un compromis de robustesse entre les objectifs ci-dessus peut être réalisé 
suivant la plage de fréquence ou la présence des signaux exogènes précédents est importante. En 
général, l’effet d’incertitudes de modélisation est souvent important en basse fréquence. De plus, 
l’erreur de poursuite doit être assez faible que possible dans la même cette plage de fréquence. 
En revanche, les effets des bruits de mesure ainsi que la dynamique non linéaire négligée sont 
souvent important en haute fréquence. C’est pour cette raison que les valeurs singulières 
maximales de la sensibilité directe S:, doivent être faibles en basse fréquence et tendent vers à 
l’unité (048) en haute fréquence. Par contre, les valeurs singulières de la sensibilité 
complémentaire T,, doivent être faibles en haute fréquence et tendent vers à l’unité en basse 


fréquence. 
1.4.3.4. Sensibilité complémentaire en entrée T,,(s) : 


Cette sensibilité représente le transfert entre la perturbation d,, et la commande u fournie par 


le contrôleur à synthétiser. Elle est liée avec S,,(s) par la relation suivante : 
T(s) =1—S,(S) (1.31) 


Ceci peut imposer un compromis entre les deux sensibilités sur l’entrée T,,(s) et S,,(s). 


1.4.4. Performance robuste (RP): 
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Le système en boucle fermée satisfait le critère de robustesse en performance si et seulement 
si le système répond aux spécifications de performance pour une classe d'incertitudes structurées 
ou non structurées du système appliquée sur le modèle nominal, jusqu'au pire cas d'incertitude du 
modèle [11,35]. Nous verrons, dans la suite de ce travail, que la vérification de la condition sur 
RP assure de façon systématique la satisfaction du compromis de robustesse entre la stabilité 
robuste et les performances nominales. La réciproque de cette implication n’est pas toujours 


vérifiée. 


1.5. Problème de sensibilité pondérée : 


Considérons le système bouclé présenté par la Fig.1.11. L’objectif de la suit de ce travail, est 


de synthétiser un contrôleur robuste satisfaisant [52] : 


e la condition sur la stabilité robuste face aux incertitudes do modélisation de type 


multplicatives amenant en sortie du modèle de synthèse ; 
e La condition sur les performances nominales ; 
e La condition sur la limitation de l’énergie de commande ; 
e __ Compromis de robustesse entre la stabilité robuste et les performances nominales. 


Il convient de noter ici que le système augmenté assurant la vérification des trois conditions 


précédentes peut être présenté par le schéma de la Fig.1.12. 


Fig.1.12. Formulation du problème de sensibilité pondérée 


AvecWs(s), W,(s) et Wz(s) sont les pondérations irrationnelles supposées fixes. Ces 
pondérations sont a priori choisies par l’utilisateur afin de limiter omax [S, jw)|, Tax EX jo)| 


Las [K (jo). S,, ( jw)|. Soit le transfert stable £,(w) qui représente le gabarit désiré assurant 
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la limitation du signal de l’erreur de poursuite (1.e., la limitation de l’effet de la perturbation). La 


condition assurant les bonnes performances nominales est définie par [13] : 
max [Sy (J@)] < IC) (1.32) 
A partir de l’expression (1.32), on peut tirer la condition ci-dessous : 
l£5* (wo) : Gmax[Sy(jo)] < 1 (1.33) 


Remplaçons, maintenant, le transfert £51(w) par celui de W;(w). L'expression (1.33) peut 


être réécrite comme suit : 
Omaxl Ws(&)] : Gmax[Sy(jo)] 1 |Ws.S,f, <1 (1.34) 


Soit le transfert stable £7(w) qui représente le gabarit désiré majorant toutes les incertitudes 


multiplicatives. La condition sur la stabilité robuste du système bouclé est donc définie par: 
Omax[Ty(jw)] < I/r(w)| (1.35) 
A partir de l’expression (1.35), on peut tirer la condition ci-dessous : 
87° (@)| : omax[7yCjo)] < 1 (1.36) 


Remplaçons, maintenant, le transfert £71(w) par celui de Wz(w), d’où l’expression (1.36) 


peut être réécrite comme suit : 
OmaxL Wr(@)] : Gmax[Ty (ja) 1 [Wr.T,f, <1 (1.37) 


Dans le but d’éviter la saturation des actionneurs installés dans la boucle de commande, la 
sensibilité décrivant l’énergie de commande doit être limitée par un gabarit, noté ici par £,(w). 


Dans ce cas, la condition sur la limitant l’énergie de commande est donnée par : 
Omax|K(j@).S,(jw)] < |£u(w)| (1.38) 


A partir de l’expression (1.38), on obtient : 


Ha! (w)| - Gmax[K(j&@).S,(jw)] < 1 (1.39) 


Remplaçons le transfert stable £,1(w) par celui de W,,(w), d’où l’expression (1.39) est réécrite 


comme suit : 
Omaxl Wa(@)] * Gmax[K(j&).S,(jo)] <1S  [M.K.S,], <1 (1.40) 


Les trois conditions (1.34), (1.37) et (1.40) peuvent être réarrangées en seul critère appelé 


« critère de sensibilité pondérée ». D’où le contrôleur robuste à synthétiser à la capacité d’assurer 
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le bon compromis de robustesse RS/NP [32] avec énergie de commande réduite si et seulement 
si la condition ci-dessous est vérifiée dans toute la plage fréquentielle. On obtient donc : 
Ws.sS,, 
WaK.S,|| = max (||Ws.5,||., [Wa K.5,1f..11Wr. 7,11.) <» (1.41) 
7 Ty 
Avec y est une valeur réelle qui doit préalablement choisir entre 0 et 1 par l’utilisateur. De 
plus, l’inégalité présentée par (1.41) est souvent appelée « problème H,, sous optimal ». Il 
peut également présenter sous forme d’un problème d’optimisation dont la solution optimale 
fournie la représentation d’état du contrôleur robuste souhaité. On obtient donc: 
Ws(jow).S, (jo) 
minx(s)ÜLFTe(P, Ko} = Minx(s) | MAX | Omax Wa (jo). K(jo).S, (jo) (1.42) 
Wr(jw).T, (jo) 
Dans certaines applications industrielles, une condition plus rigoureuse que celle de RS/NP 
est imposée lors de la synthèse du contrôleur robuste à synthétiser. C’est la condition sur RP [35] 


qui est définie par la condition suffisante ci-dessous : 
Wr.T, fl. +1Ws.S, 1 < 1 (1.43) 

En tenant compte la relation suivante : 
max (IWs.s, We. K.S, ||, wr. 7,1.) < ]|Ws.s, 1, + [Ma KS, fl + WT. TL 44) 


D'où, la satisfaction de la condition sur RP conduit automatiquement à satisfaire le bon 
compromis de robustesse RS / NP, quelles que soient les incertitudes du modèle affectant le 


modèle nominal. 


1.5.1. Solution du problème H,, : 
Considérons le schéma de la Fig.1.13 qu’à partir duquel le transfert du système augmenté est 


représenté par le bloc ci-dessous : 


Fig.1.13. Transfert du système augmenté (modèle généralisé) 


D'après la Fig.1.13, les entrées exogènes sont liées par les sorties exogènes comme suit : 
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VA Pin © Puy a) 
(-) 2 ( ‘ *]() (1.45) 
y Poi © Pr u 


La représentation d’état du modèle augmenté est définie par : 


S * I — À : (Bi B:2) 
(&) 
C2 


D'où, la représentation d’état décrivant de ce modèle est donnée par: 


h 1.46 
( Du) (98) 
D;: D;2 


X(t) = Ax(t) + Biw(t) + Bu(t) 
Z(t) = Cix(t) + Disw(t) + Duu(t) (1.47) 
VC) = Cx(t) + D1w(E) + D,zu(t) 
Considérons, maintenant, le schéma de la Fig.1.14 qu’à partir duquel le transfert du 


contrôleur à synthétiser est représenté par le bloc ci-dessous: 


Fig.1.14. Représentation d'état du contrôleur à synthétiser 


En tenant compte la relation y = K * u4. D’où la représentation d’état du contrôleur est 


donnée par : 


ee = Azxx(t) + Brux(t) (1.48) 


YkE) = Cexe(t) + Drux(E) 


Le système d’interconnexion reliant le transfert du système augmenté avec celui du contrôleur 


est présenté par la Fig.1.15 comme suit [33] : 


Fig.1.15. Système d'’interconnexion décrivant le système en boucle fermée 
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D'où, la représentation d’état du système en boucle fermée est obtenue à partir des équations 


(1.47) et (1.48) comme suit [33, 53, 54]: 
Ge) { G 4) : Go) : (6 a) | /G) 
D | G 0: Du @e O |: &(®) (1.49) 
Cey \ 0) : Go) ? où 0) Gt) 


En tenant compte les notations ci-dessous : 


| é 1) (:) Ce :) 
A : B : B 0 A4 0 0 Bk 
CE, : Du : Du l= (C1 0) : Dis : (Di 0) (1.50) 
Bi Pa 52e (6): (oi) : ( 0) 
C2 0 D; Dy> 0 
On aura donc: 
A=(( :) Dis = Dis 
AU Ak Ci = (Ci 0) Di2 = (Di 0) 
B 0 
B=(0) le(s Guy Da=(},) (1.51) 
5, = ( 4 Ga 0 ; a D) 
| 0 Bk 27 \D; 0 
A partir de l’équation (1.49), on peut tirer les relations ci-dessous: 
Yr “ Yk 
(y)= CC Da) (ue) + De (7) (1.52) 
x 
Da) (7) = (C2 Du(sx) (1.53) 
a) 


De (1.52) et (1.53), la matrice (1 — D;2) doit être inversible. Cela permet de définir les deux 


expressions suivantes : 


() = (1—D;>) 1(C> D) (x) (1.54) 
(eo) 
(x) = CE . )+ o }a- 2297 D2)| (x) (1.55) 
Ci Di D:2 mn) 
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D'où, l’existence du contrôleur robuste est liée directement par la vérification des deux 


conditions nécessaires ci-dessous : 


ll —D 
e La matrice (J = Da2) = D . doit être inversible. 
— D22 
A 0 B; 0 1  —D;\ "f0 C 
_ ie — 2 k k\ 4 
e La matrice À + B,(1 — D;:) C) a)+( a) es ñ ) É o°) doit 


être stable. 

Dans ce qui est suit, nous allons présenter deux algorithmes d’optimisation avec lesquels la 
représentation d’état du contrôleur robuste est assurée par la commande H,,. Ces algorithmes 
sont les suivants : 


1.5.1.1.Méthode X,, basée sur la résolution des deux équations algébriques de Riccati : 
Supposons ici qu’il existe un contrôleur robuste K(s) stabilisant le système augmenté P(s). Cela 
implique que : 

e Les couples (4, Z.) et (A, #)) soient stabilisables 

e Les couples (4, C;) et (4, C;) soient détectables 


Supposons également que : 


e _rang(D:2) égal le nombre d’entrée u 
+ rang(D;:) égal le nombre de sortie y 
© Di = Da = 0 


e P,,(s)et P;.(S) n’ont pas des zéros sur l’axe imaginaire. 
12 21 


Théorème : 


La détermination de la représentation d’état du contrôleur robuste K(s) dépend de l’existence 


des deux matrices définies positives X,k > 0 et Y, > 0 qui représentent la solution des deux 


équations de Riccati ci-dessous : 
ATX + XA + X(y72B,BT — B,B1 )X + CIC, = 0 (LS 
AY + YAT +Y(y72CTC, — CTC,)Y + B,BT = 0 (1.56) 


Tenant compte le rayon spectral p(X.., Y.) qui doit vérifier la condition p(X%, Yo) < Y?°, 


d’où, les solutions X,, et Y, peuvent être données par : 


A — B}D1C1 7 ?B1 B1 — B2Bz 


Xo — Ric T T T TPNT 
—(( —D:2D12)C1) (A — Di2Di2)C1) —(A — BD12C1) 


(1.57) 
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= Ri A — B,D;1C2 p Ci Ci = CC 


ic (1.58) 
—(Bi( — D21D71)B1Q — D21D21)) ((U — Di2D12)C1)  — (A — B1D21C2)" 


La valeur optimale de y € [Ymin Ymax], (notée ici par yon) est obtenue de façon itérative par 
l’algorithme de dichotomie dans lequel la distance dy, = as : il est décroit dans chaque 
itération k. La valeur finale y,,+, qui est fournie par cet algorithme à l’itération #, représente la 


médiane de l’intervalle finale, on aura donc: 


1 
Vopt — 5 Wminy na Ymaxy) (1.59) 


De plus, la représentation d’état de K(s) est obtenue par l’application de LFT inferieure sur le 


système d’interconnexion {K,, b} présentée par la Fig.1.16, sachant la condition ci-dessous : 


Fig.1.16. Ensemble des régulateurs qui résous le problème de synthèse H, 


En tenant compte la condition ci-dessous : 
0 < IMOIE + Vopt (1.60) 


On obtient donc une infinité de contrôleurs stabilisant dans l’intervalle [0 Pc Sachant 
que le contrôleur qui correspond le cas particulier D(s) = 0 est appelé « le contrôleur central », 
noté par K,(s). D’où sa représentation minimale est donc déterminée en fonction des trois 
matrices ci dessous : 


e La matrice de retour d’état F, qui est donnée par : 

F = -(DT, : C1 + BI * Xe) (1.61) 
e La matrice d’observation H, qui est donnée par : 

H=-(B,: 04 +70) (1.62) 
e La matrice de mise à l’échelle Z,,,, qui est donnée par : 


PART CES OCR 2 1 Du (1.63) 
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Sachant que la représentation d’état souhaitée du contrôleur robuste est définie par : 


(s°1=24,,) 5 Bya 


KA (s) = | (1.64) 


Cra : Dya 
D'où, Aya» Bras Cka ©t Dxa Sont données comme suit : 


Axa = A+B'F+y 2: B :B-Xe +270" H:(C +7 2: D: BT'XS) (1.65) 


Bya —= = H Zoo(B2 + | Re DL (E : D:2)] (1.66) 
F 
re | : 1.67 
de —(C, +7 2: D" BT' Xe) ( ) 
RE 
Den = L : (1.68) 


1.5.1.2. Méthode X.,, utilisant la formulation LMIs: 


La technique LMIS permet de résoudre le problème de sensibilité pondérée par des 
algorithmes d’optimisation sous contraintes. La détermination de la représentation d’état du 


contrôleur robuste est donc basée sur le théorème ci-dessous [24, 26] : 
Théorème: 


La détermination de la représentation d’état du contrôleur robuste K(s) dépend de l’existence 
des deux matrices définies positives X,. > 0 et Y, > 0 qui représentent la solution du problème 


LMIS ci-dessous : 


on ee . jo 


ke u BXS —y.1 x. Ce ‘)<0 (1.69) 
[C1 Dai] 2 
Tr [Mo AT + AY Yo: si] B; 
co ‘) | Ye —7y. 1 D Ce )<0 (1.70) 
[C1 Du] [—y.1] 
CR 
( 0 (1.71) 


Où N, et N,, sont les bases des noyaux de (85 Dl,)et (C2 D21) respectivement. 
1.6. Simulations et interprétations des résultats : 


Dans le but d’évaluer les performances fournies par un système bouclé par un contrôleur 


robuste synthétisé par la commande H, utilisant l’approche LMISs, nous allons commander un 
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processus hydraulique avec lequel le niveau d’eau d’un canal d’irrigation à surface libre est 
commandé par une tension alimentant une électrovanne, qui est installée au début de ce canal. A 
cet effet, le but de la commande robuste est de transmettre de l’eau de sa source amont (le 
barrage, la rivière) à ses utilisateurs finaux aval (des agriculteurs). Il est à noter ici que le 
comportement réel de ce système physique est fortement non linéaire. La linéarisation est donc 
effectuée en plusieurs points de fonctionnement et cela par l’utilisation d’un modèle incertain 
présenté par le transfert rationnel suivant: 


K-e-0's 
1+T:54 


G(s) = 


(1.72) 

Lors de la modélisation du comportement réel de ce système, le gain statique K du modèle 
incertain est changé dans l’intervalle [K,,51, Kmax | = [0-2648, 0.4145 ]. La constante du temps 
T est variée également dans l’intervalle [Tin , Tnax | = [2-5002,3.3620 |. De plus, le retard pur 
0 peut prendre n’importe quelle valeur dans l’intervalle [din , max | = [43262 ,4.9746] 
ainsi que l’ordre fractal a, qui est porté sur l’opérateur de Laplace, peut être changé dans 
l'intervalle [tin , Amax | = [0.-6408 ,0.8740 ]. Sachant que la synthèse basée sur la commande 
H nécessite a priori le développement d’un modèle linéaire d’ordre entier. Pour cela, le terme 
T : s*® est préalablement approximé par un transfert irrationnel d’ordre n où ses pôles et ses zéros 


sont obtenus par l’approximation de Oustaloup. On aura donc la structure générale ci-dessous: 


a _ Hn1(S-Zi) 
T's — TE, G-P) (1.73) 
Le détail de cette approximation sera détaillé en chapitre 3. De plus, le terme eŸ" est 


préalablement approximé par un transfert irrationnel stable d’ordre m où ses paramètres sont 
obtenus par la méthode de Padé, qui est disponible dans la bibliothèque du Matlab. On aura 
donc la structure générale suivante : 
gs \X 
+5 |( 7) à 


k 


12 [(5) à) 


Supposons ici que le modèle nominal possède la même structure que celui du modèle 


e 7 — (1.74) 


incertain. D'où le transfert décrivant ce modèle est donné comme suit : 


Gn(s) un 


Où les paramètres K, , 0, , T, et &, intervenant dans la construction de ce modèle sont 


Knre”0ns 
1+Tn°s4n 


(1.75) 


obtenus comme suit : 
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K _ Kmax+K min T. _—_ Tmax+Tmin 
n n— 


,et - (1.76) 


__ 0max+0min _ Amax*Amin 


n— 2 An 2 
De (1.75) et (1.76), la fonction de transfert du modèle nominal est donnée par : 


0.33969-e 74-6504 


GS) = L9812507877 (70 
En tenant compte les deux approximations : 
. <0.75747 _ 30.491-(5+1.136)(s+0.04121)(5+0.001495) 
2.9312:5 = (s+14.01)(5+0.5082)(5+0.01844) 1.78 
_4.65045 _ -(5+0.9987)(52-1.582:5+1.195) (1.78) 
É 7 (s+0.9987)(52+1.582-5+1.195) 
On obtient donc un modèle d’ordre 6 présenté par le transfert équivalent ci dessous: 
—0.0107-(5+14.01)(s—0.998)(s+0.508)(s+0.018)(s2—1.582-5+1.195 
ON te Rte) (1.79) 


(s+1.404)(5+0.998)(s+0.181)(5+0.016)(52+1.582:5+1.195) 


Sachant que la synthèse du contrôleur robuste nécessite a priori l’évaluation de toutes les 
incertitudes de modélisation agissant sur le modèle (1.79). A cet effet, La forme multiplicative 
est choisie pour prendre en considération la variation susceptible de chaque paramètre du modèle 
et les approximations nécessaires du retard pur et l’ordre fractal survenant dans le transfert du 


modèle de synthèse. 


Il est à noter que ces incertitudes peuvent fournir une classe des régimes perturbés où chaque 
régime perturbé peut produire une erreur relative par rapport au modèle nominal de synthèse. 
L’analyse fréquentielle est donc effectuée par le tracé des lieux des valeurs singulières 
maximales de ces erreurs dans le plan fréquentiel. Ces dernières sont ensuite majorées par le 
transfert stable W(s), qui est choisie par l’utilisateur. Ce transfert permettant de fournir la 


condition sur RS du système en boucle fermée. 


La Fig.1.17 présente les valeurs singulières maximales de 20 formes d’incertitudes 


multiplicatives directes en sortie ainsi que les lieux de GnaxlLWr(jæw)] majorant celles-ci. 
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Erreurs relatives 
20 — ï TT 
Spécification sur 
= 19 I tabilité bust 
CA a stabilité ro met || 
> 
& Oo 
SDS -3 
Æ 
ES 
-10 Un — Æ 
= = = —= 
A = = 
LL = = = 
“© 20Ë sÉ 
— 
TD 
li == É — 
e = = = 
ae -# +" 
a) =. 
(as) 
7 40 RE NN 1—< / A 
\ / W % 
Kg v A 
-50 i i i Fr it g L 2 in 
1072 1071 10° 107 
Pulsation (rad/s) 


Fig. 1.17. Lieux des valeurs singulières maximales de Omax [a Gw)] et OmaxLWr(jw)] 


D’après la Fig.1.17, les lieux des valeurs singulières maximales de toutes les incertitudes 
multiplicatives sont assez faibles dans les basses fréquences, en particulier dans la plage de 
fréquence w € [1072, 0.91586 ] radians/secondes, là où les performances nominales soient 
exigées. En revanche, Ces incertitudes deviennent importantes dans les hautes fréquences, là où 
les dynamiques non linéaires sont négligées. Sachant que la majoration des incertitudes de 


modélisation est aboutie comme suit : 


Omax As, (&)] — Omax [(Go Ga) s GnG&)) ‘ GuGa)| < GmaxlWr(jw)] (1.80) 


Avec k = 1,:,20 et w € [1072, 10] radians/secondes. La spécification qui assure la 


condition (1.80) est donc donnée par le transfert stable du 1° ordre suivant : 


W(s) = 4.2482-(s+0.4218) 
E s+3.891 


(1.82) 


Afin de prendre en compte les spécifications nominales dans la synthèse du contrôleur H,, 
primaire, on augmente le modèle nominal (1.79) par les gabarits de pondération fréquentielle 
suivants : 

e La contrainte de rejet de perturbation est prise en compte en pondérant le transfert de 
sensibilité S,(s) par le gabarit de pondération fréquentielle W,(s) choisie par le transfert 


stable du 1°7 ordre suivant : 


0.4-(5+0.15) 
S+6:10—° 
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e Cette contrainte s’exprime en imposant au gain de W,(s) d’avoir un comportement de type 
intégrateur en basse fréquence. Le lieu des valeurs singulières maximales de cette pondération 


est représenté sur la Fig.1.18 comme suit : 


Spécification sur les Performances Nominales 
A0 T T T T 


Valeurs Singulières Maximales (dB) 
a 


-10 1 1 1 1 
107$ 107? 1071 10° 10! 
Pulsation (rad/s) 


Fig. 1.18. Lieux des valeurs singulières maximales de Onax|Ws(jw)] 


e Les contraintes de robustesse la limitation de l’énergie de commande se traduisent par le 
gabarit de pondération W,,(s), portant sur le transfert de sensibilité KS,, et dont le lieu des 


valeurs singulières maximales représenté sur la Fig. 1.19 comme suit : 


Spécification sur la limitation de l'énergie de commande 
20 TT T T T TT 


Valeurs Singulières Maximales (dB) 
[e] 


_-20 it 1 : 1 it 
107$ 1072 1071 109 101 
Pulsation (rad/s) 


Fig. 1.19. Lieux des valeurs singulières maximales de Omax|Wi(jw)] 


Sachant que la pondération W{y(s) est choisie comme suit : 


15.152:(s+0.5883) 
s+50.4 


Wi(s) = (1.84) 


Le système augmenté est obtenu à partir des transferts irrationnels de ces pondérations ainsi 


que celui du modèle nominal. Le problème de sensibilité pondérée est ensuite formulé dans le 
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but de satisfaire un bon compromis de robustesse NP/RS avec une limitation de l’énergie de 
commande. Ce problème est résolu par la méthode Æ,, basée sur l’approche LMISs où sa solution 


fournie un contrôleur H., primaire d’ordre 9 et un niveau de performance égal à y = 0.945. 


La représentation d’état de ce contrôleur est finalement convertie sous forme d’un transfert 


d’ordre entier présenté comme suit : 


ro 70.35-(5+50.48)(s+3.89)(5+1.4)(s+0.99)(s+0.18)(s+0.016)(52+1.582-5+1.195) 
7 (s+1252)(5+3.85)(s+0.51)(5+0.02)(5+5.9-10-5)(s2+0.51:5+0.91)(52+3.98-5+5.58) 


(1.85) 


Par conséquent, La condition sur SR du système bouclé par le contrôleur H,, primaire est 
vérifiée dans toute la plage de fréquence «w € [1072, 10] radians/secondes. La Fig.1.20 
représente le lieu des valeurs singulières maximales du transfert de la sensibilité complémentaire 
T,, et celui de l’inverse de la spécification sur SR. A partir de la Fig.1.20, il est très claire de 
confirmer que la condition Gmax [T, ( jo)| < (Omax|Wr(jw)l)"?. Ceci est traduit, dans le plan 
temporel, que la stabilité du système en boucle fermée est garantie pour toutes les perturbations 
multiplicative sur le fonctionnement nominal du système réel, jusqu'au pire cas des incertitudes 


du modèle. On obtient donc la Fig.1.20 comme suit : 


Condition sur la Stabilité Robuste 
10 T T T T 


omax (T_y) 
omax(1/W_T) 


-30 


Valeurs Singulières Maximales (dB) 


-50 


-60 1 1 1 1 
10 1072 1071 10° 10! 
Pulsation rad/sec (rad/s) 


Fig. 1.20. Condition sur RS du système bouclé par le contrôleur H, primaire 


De plus, La condition sur NP du système en boucle fermée est vérifiée par le 
contrôleur H, primaire et cela en toutes les fréquences qui se situent dans la plage w € 
[1072, 10] radians/secondes. La Fig.1.21 représente le lieu des valeurs singulières maximales de 
la sensibilité directe S,. Ce dernier est comparé par l’inverse de la spécification sur NP. On 


obtient : 
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Condition sur les Performances Nominales 
10 F T TTTTTITI T h TOFTTTTTI T T FTTTTTT T T TFTTTITI 


omax(S _d) 
omax(1/W_S) 


Valeurs Singulières Maximales (dB) 


-A0 i ii ii ii : ii ii iii : CORRE A ER ES ET : PA er 
107$ 102 10°! 100 10! 
Pulsation rad/sec (rad/s) 


Fig. 1.21. Condition sur NP du système bouclé par le contrôleur H,, primaire 


D’après la Fig.1.21, il est très claire de confirmer que la condition oyyx [S, jo) | < 
(OnaxlWs(jw)])"{. Ceci est traduit, dans le plan temporel, par l’assuré d’une bonne dynamique 
de poursuite caractérisée par un bon suivi de la trajectoire de référence dans lequel l’erreur de 
poursuite est très réduite en régime permanant. De plus, le système bouclé devient insensible aux 
effets des bruits de mesure et les modes inconnus qui sont négligés lors de la phase de 
linéarisation du système physique. 

Il convient de noter ici que ces performances sont assurées seulement en absence des 
incertitudes de modélisation. Concernant, maintenant, la limitation de l’énergie de commande, le 


contrôleur H., primaire présente l’avantage d’assurer tous les objectifs précédents avec énergie 


de commande très raisonnable en termes du coût de commande et saturation des actionneurs. 


La Fig.1.22 représente le lieu des valeurs singulières maximales de la sensibilité K + S:, qui 


est comparé celui de l’inverse de la spécification sur la limitation de l’énergie de commande. On 


obtient : 


Page | 37 


CH.1: COMMANDE ROBUSTE UTILISANT LES PONDÉRATIONS IRRATIONNELLES FIXES 


Condition sur le signal de commande 
20 : LR ELLE ï TT TC | EP E r TT TT] 


omax(KS_d) 
omax(1/W_u!| - 


15 


Valeurs Singulières Maximales (dB) 
O 


5+ 

-10 + 

-15 + 

_-20 à nl de joints ii iii iiul in iii 
10 102 10°! 10° 10! 


Pulsation rad/sec (rad/s) 


Fig. 1.22. Condition sur la limitation de commande du système bouclé par le contrôleur H,, primaire 
D’après la Fig.1.22, il est très claire de confirmer que l’objectif de commande est atteint avec 
une commande moins coûteuse. Ceci à partir de la vérification de la condition 
suivante GouxlR SC jo)| < (GmaxlWa(jw)])"*. La partie suivante à pour but d’examiner est 
ce que le contrôleur H,, primaire est capable d’assurer un bon compromis de robustesse NP\RS 
ou non? A cet effet, nous devons tracer le lieu des valeurs singulières maximales du critère 
suivant max (IWs.S, IL I Wz. T1.) et cela dans toute la plage de fréquence w € [1072, 10] 


radians/secondes. On aura donc la Fig.1.23 : 


Compromis de robustesse 
-0.48 T T T 
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o 
01 
N 
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107$ 1072 1071 10° 10! 
Pulsation rad/sec (rad/s) 


Fig. 1.23. Vérification du compromis de robustesse NP\RS du système bouclé par le contrôleur H,, primaire 
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D’après la Fig.1.23, on peut constater que le lieu de omux [max (IlWs.s, 1 | 


Wr.T,ll,.)| se 
situe au dessous de l’unité (ie., O3) dans toute la plage de fréquence [1072, 10] 
radians/secondes. Ceci peut nous confirmer que la synthèse par la méthode H,,basée sur 
l’approche LMISs s’avère très efficace pour résoudre ce genre des problèmes. L’objectif de la 
suite de ce travail est de vérifier est-ce que ce contrôleur est capable de conserver les 
performances observées précédemment en présence des incertitudes de modélisation? Pour cela, 
on doit tracer le lieu de la condition la plus rigoureuse que celle du compromis de 


robustesse NP\RS. C’est la condition RP qu’on doit la vérifier à travers le tracé du lieu des 


valeurs singulières maximales du critère suivant (IlsS, IL, + wr.T,1.). On aura donc la 


Fig.1.24 : 
Condition sur les Performances Robuste 
2 r T T 
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Fig. 1.24. Vérification de la condition RP du système bouclé par le contrôleur H,, primaire 


D’après la Fig.1.24, on peut voire nettement la violence de la condition RP dans les basses 
fréquences, en particulier dans la plage [1072, 0.0761] radians/secondes. Ceci se traduit, dans le 
plan temporel, par la détérioration des performances du système corrigé à cause les effets des 
incertitudes de modélisation. 

Dans la suite de ce travail, nous allons présenter les réponses temporelles du système corrigé 
en présence 192 régimes perturbés avec lesquels on a pu tracer toutes les erreurs relatives 
présentées dans la Fig.1.17. A cet effet, la boucle de commande présentée par le schéma bloc de 
la Fig.1.11 est excitée par les trois entrées exogènes suivantes : la première entrée représente la 
référence (à suivre) imposée sous forme d’un échelon unitaire. Cette entrée est appliquée au long 
de l’intervalle du temps t € [0, 450 +xT,] secondes avec le pas de discrétisation T, = 0.15 
secondes. La deuxième entrée représente la perturbation amenant en sortie de modèle de 


synthèse. Elle est imposée sous forme d’un échelon d'amplitude 0.25 (1.e., 25% de l’entrée de 
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_ Le . : 250 " : : 
référence) appliqué à partir de l’instant du temps t = a secondes. La dernière entrée représente 
e 


le signal des bruits de mesures caractérisé par la valeur moyenne nulle et la variance Var — 


0.01. Ce signal est appliqué à l’instant du temps t = — secondes. La Fig.1.25 représente les 


sorties globales fournies par 192 régimes perturbés. En revanche, la Fig.1.26 représente la partie 


transitoire des réponses indicielles fournies par le contrôleur H,, primaire. 


Fig. 1.25. Sorties fournies par le contrôleur H,, primaire pour 192 régimes perturbés. 


D’après la Fig.1.25, on peut remarquer clairement que la dynamique de poursuite de la 
trajectoire de référence est assurée par la présence des erreurs statiques avant l’excitation par la 
perturbation d.,, et après l’injection du signal des bruits de mesure. De plus, la réponse fournie 
par le pire des régimes perturbés est caractérisée par un temps de montée égal à th — 
T5% 


26.33 : — 150 4 Ô 
Es secondes, un temps de réponse égal à t,.., = = secondes, et un dépassement maximal égal 
e e 


à Dax = 12% (voir la Fig. 1.26). 
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FF 


Fig. 1.26. Agrandissement de la partie transitoire de la dynamique de poursuite fournie par le 
contrôleur H,, primaire pour 192 régimes perturbés. 


La Fig.1.27 représente les commandes assurant les sorties globales citées précédemment. On 


obtient : 


IL! 


Fig. 1.27. Commandes fournies par le contrôleur H., primaire pour 192 régimes perturbés. 


D’après la Fig.1.27, on peut remarquer que toutes les commandes fournies par le 
contrôleur H, primaire sont insensibles aux effets des bruits de mesure en régime stationnaire. 
En revanche, le pire cas de ces commandes peut atteint un pic Umax — 4.92 en régime 


transitoire. 


1.7. Conclusion 
Dans ce chapitre nous avons présenté les notions théoriques de base utilisées lors de la 


synthèse du contrôleur H,, primaire. Ce dernier est appliqué sur les systèmes physiques 
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modélisés par des modèles incertains. Nous avons également classé les différents types 
d’incertitudes agissant sur le modèle de synthèse. L’analyse fréquentielle de ces incertitudes à 
donner la possibilité de fournir la spécification sur la stabilité robuste du système en boucle 
fermée. Cela est assuré au moyen du théorème du faible gain. Sachant que l’inverse de cette 
spécification permet de limiter l’évolution des valeurs singulières maximales de la sensibilité 
complémentaire dans le plan fréquentiel. À cet effet, la condition sur la stabilité robuste est 
assurée à partir d’une satisfaction d’une inégalité fournie par cette pondération. Parallèlement, 
une autre condition sur les performances nominales est nécessaire lors de la synthèse du 
contrôleur robuste. Cette condition dépend du choix d’une pondération où son inverse permet de 
limiter les valeurs singulières maximales de la sensibilité directe du système corrigé. De plus, 
une autre condition de robustesse doit être prise en considération lors de la synthèse du 
contrôleur robuste. C’est la condition assurant un coût réduit de la commande ainsi que la non- 
saturation des actionneurs de la boucle. La vérification de ces trois conditions de robustesse est 
impossible dans la même plage de fréquence. C’est pour cette raison qu’on doit chercher à 
réaliser un compromis de robustesse dans des fréquences bien déterminées. Ce compromis est 
ensuit formulé sous forme d’un problème appelé « problème de sensibilité pondérée » et sa 
solution optimale est assurée par deux algorithmes avec lesquels la commande H,, peut fournir la 
représentation d’état du contrôleur H, primaire. Ce dernier a été utilisé pour la commande d’un 
système hydraulique simulant le fonctionnement d’un canal à surface libre. Le comportement 
réel de ce système a été modélisé par un modèle incertain comportant un terme d’ordre non 
entier et un retard pur. L’analyse des incertitudes est effectuée via le choix de la forme 
multiplicative et l’objectif de la commande réside dans la satisfaction de plusieurs condition de 
robustesse, en particulier celle de RP . Les résultats de simulation montrent que le choix de la 
méthode H,, utilisant l'approche LMIsest très utile quand il s’agit de satisfaire un bon 
compromis de robustesse entre RS/NP. Par contre, la robustesse sur les performances robustes 
est très difficile à aboutir par le choix des pondérations irrationnelles fixes. C’est pour cette 
raison qu’on doit les remplacer par d’autres pondérations plus flexibles et ajustables. La 
sélection des ces dernières doit être assurée via un algorithme d’optimisation et la combinaison 
de ce dernier avec la commande H, peut améliorer la qualité de commande fournie par le 


contrôleur à synthétiser. Cette amélioration faisant donc l’objet du chapitre suivant. 
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Chapitre 2 


Commande robuste utilisant les pondérations 
irrationnelles ajustables 


1. Introduction : 


Dans ce chapitre, nous allons présenter une solution alternative qui peut être améliorée les 
performances fournies par la commande A, basée sur les pondérations irrationnelles fixes. Cette 
solution sera basée sur l’algorithme génétique lors de la synthèse du contrôleur 4, robustifié. 
Cet algorithme permettra de fournir la sélection automatique des pondérations irrationnelles 
possédant des paramètres ajustables optimisés préalablement par l’algorithme génétique. Cette 
optimisation sera fondée sur la satisfaction de quelques règles de réglage qui existent déjà dans la 
littérature. Ces dernières seront après utilisées pour augmenté le modèle nominal de synthèse 
afin de produire le système augmenté qui avec lequel on pourra ressortir le problème de 
sensibilité pondérée en utilisant la transformation linéaire fractionnaire. Ce problème sera résolu 
par la méthode H,, basée sur l’un des deux algorithmes cités en premier chapitre dont la solution 
optimale permettra de fournir la représentation d’état du contrôleur robuste assurant la mise à 
jour du contrôleur H,, primaire. Ces étapes seront se réitèrent en tant que la condition sur RP 


n’est pas satisfaite. 


Il est à noter que la satisfaction de la condition sur RP dépend fortement du respect des 
règles directrices existantes. Ceci permet de transformer le problème de la sélection automatique 
des pondérationns irrationnelles ajustables en un problème d'optimisation sous contraintes dans 
lequel les paramètres optimaux de ces pondérations peuvent être situés dans un espace de 
recherche, ce qui est imposé par des limites inférieure et supérieure contraignant la solution 
optimale souhaitée. Sachant que la fonction de fitness du problème d'optimisation sous 
contrainte représente la condition RP. Le processus de minimisation garantissant par l'algorithme 
génétique maximisera donc les marges de robustesse sur RS, NP et RS\NP, ce qui améliorera 


ultérieurement la marge de robustesse sur RP. Ceci faisant l’objet de ce chapitre. 


Il est nécessaire de mentionner que les ordres de poids ajustables irrationnels sont des 
inconnues entières. Sachant que l'algorithme génétique a l'avantage d'optimiser ce type de 
variables. Pour cette raison, la robustification demandée utilisera cet algorithme pour 
sélectionner automatiquement les pondérations requises pour satisfaire la condition sur RP. 


Finalement, ce chapitre sera achevé par une partie de simulation dans laquelle les performances 
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fournies par le contrôleur H,, robustifié seront comparées, dans les plans fréquentiel et temporel, 


avec celles fournies par le contrôleur H, primaire cité en premier chapitre. 


2.1. Règles de réglage des pondérations irrationnelles ajustables 


En général, la synthèse du contrôleur H,, robustifié à partir des pondérations irrationnelles 
ajustables nécessite a priori la connaissance de l’impact des signaux d’entrées exogènes sur la 
boucle de commande. Pour cela, nous devons définir les différentes sensibilités fournies par le 
système bouclé. A cet effet, le schéma de la Fig.2.1 permet de définir les expressions de ces 
sensibilités à partir des entrées exogènes agissant sur l’entrée et la sortie du modèle de synthèse 


[34, 55]. On obtient donc : 


Fig.2.1. Les différentes fonctions de sensibilités fournies par le système bouclé 


D'après la Fig.2.1, 1l est claire de remarquer que l’entrée et la sortie du modèle de synthèse 
sont liées par les entrées exogènes via six sensibilités qu’à partir desquelles on distingue deux 


classes des performances, à savoir : les performances d’entrée et les performances de sortie. 
a. Performances d’entrée : 


Elles représentent l’impact des signaux d’entrées exogènes sur le signal de commande et c’est 
pour cette raison que la satisfaction de ces performances est conditionnée par la vérification des 


conditions de robustesse ci-dessous : 


e Le lieu de Gmax [K (jw) : S,Gw)| doit être minimisé dans toute la plage de fréquence. Ceci 
permet d’éviter la saturation des actionneurs et éviter la détérioration des instruments de 
mesure installés dans la commande. 

e Le lieu de OmaxlK (jw) : S,(jw)] doit être minimisé afin de réduire l’effet des bruits de 


mesure 7, sur le signal de commande, notamment en haute fréquence. Ceci permet d’assurer 
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une bonne atténuation de l’effet des incertitudes sur la commande, en particulier celles de type 
multiplicatives non structurées. 
e Le lieu de omaxlS, (Üjw)] doit être minimisé afin de rejeter l’effet de la perturbation d,, sur le 


signal de commande. 
b. Performances de sortie: 


Elles représentent l’impact des signaux d’entrées exogènes sur le signal de sortie du système 
en boucle fermée. À cet effet, Ces performances sont assurées si le contrôleur robuste à 


synthétiser à la capacité de fournir les propriétés ci-dessous : 


e Un bon suivi de la trajectoire de référence r qui peut être assuré quand le lieu de 
Pas [É, (jw)| est faible dans tout la plage de fréquence ; 

e Un bon rejet de la perturbation d,, qui peut être assuré si le lieu de oax [S,Gow)| est assez 
faible dans tout la plage de fréquence ; 

e Une bonne suppression de l’effet de la perturbation d,, sur la sortie y, qui peut être aboutie si 
le lieu de Omax [S, Go). G(jw) ] est faible dans toute la plage de fréquence ; 

+ Une bonne atténuation de l’effet des bruits de mesures 7, et l’effet des modes négligés en 
haute fréquence, qui sont garanties si lieu de oux [ÎT, (Gw)| est faible dans toute la plage de 
fréquence ; 

A partir de ces conditions, le choix des gabarits limitant ces sensibilités peuvent être effectué 

comme suit : 


2.1.1. Sélection du gabarit limitant le lieu de oax [Ty Gw)|: 


En général, la limitation du lieu de Omax [Ty Gw)] est indispensable pour sécuriser le 
système bouclé contre les effets des bruits de mesure et les dynamiques négligées en haute 
fréquence [34, 55] . L'amélioration de la marge de robustesse sur RS est donc liée directement 
par le choix adéquat de la pondération irrationnelle ajustable W,(s). En effet, ses paramètres 
optimaux doivent être optimisés de tel sort que le lieu des valeurs singulières maximales de son 
inverse est toujours se situe au dessous de la sensibilité complémentaire souhaité T.,(s). Par 
conséquent, la forme idéale de cette sensibilité peut être aboutie si le lieu de omux [T, Ga) | est 
s’abaissé autant que possible en haute fréquence [56-58], puis il s’acroit progressivement jusqu'à 
l'unité (1.e., 0 dB) en basse fréquence. Cette dernière doit être assurée en fonction de la 
décroissance de la pulsation w. La Fig.2.2 représente l’allure désirée de Onax [ÎT, Gw)| anisi que 


le gabarit fréquentiel correspondant (Gmax[Wr(jæ)])"!. On obtient : 
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1 
Omax [Wr] 


Fig.2.2. Forme désirée de Omax [T, | comparée avec celle de (Gmax|WrD71 


En général, le transfert complexe définissant la pondération W,(s) est donné par [34] : 
+, 7 
Wr(s) = [ST (2.1) 


Sachant que : 
e Mr E R‘* : représente le gain maximal de 6x Wr(jæ)] obtenu en basse fréquence ; 
e Er E R‘* : représente le gain minimal de GmyxlWr(jw)] obtenu en haute fréquence ; 
e pr € R'* : représente la pulsation où le gain de OaxlWr(jw)] = —3dB ; 
e ne N'*: représente l’ordre entier de la pondération W(s). 

En général, la pondération fixe W,(s) est systématiquement obtenue via la majoration de 
toutes les incertitudes multiplicatives non structurées. En conséquence, les paramètres de cette 
pondération sont ainsi obtenus pour définir les limites supérieure et inférieure contraignant les 
paramètres de la pondération 1irrationnelle ajustable souhaité. Nous verrons plus tard comment 


choisir chaque paramètre de cette pondération souhaitée de telle sorte que la marge RS puisse 


être bien améliorée. 


e Choix du gain maximal M,: 


ES 


En général, le gain maximal Mr représente souvent le pire cas à atteindre par le lieu 


dde [ÎT, (ja) | en basse fréquence. Il est exprimé par l’équation (2.2) ci-dessous: 


Mr = |A — MAX (Gmax [Ty G ‘ ©min)|) (2.2) 
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De l’équation (2.2), Il est nécessaire de décroitre le paramètre M4 autant que possible. Ceci 
permet d’apparaitre des pics néanmoins faible dans le lieu de omax[T (jw)]. En général, le choix 
de ce paramètre est souvent effectué de telle sorte que la marge du gain MG et celle de la phase 
MP du système en boucle ouverte doivent satisfaire les deux conditions ci-dessous [34]: 

d 
MG >1+ ER 
à (2.3) 
MP > 2 arcsin(——) 
2Mr 
Dans la théorie de la commande robuste, le paramètre mentionné précédemment est souvent 


choisi dans l'intervalle M} € [1.25 ,2.75]. 
e Choix du gain minimal £; : 


En général, 1l est préférable de mettre le gain minimal £r égal à zéro. Ceci permet de réduire 
l’influence des bruits de mesure sur la sortie du système bouclé notamment en haute fréquence. 
En revanche, le choix de £7 = 0 est pratiquement impossible à cause des conditionnements 
numériques imposées dans la méthode Æ,, et c’est pour cette raison qu’on doit le choisir comme 


suit [59]: 
0 = €T KC@ My (2.4) 
e Sélection de la pulsation wgr : 


En général, 1l est préférable de translater cette pulsation w47 graduellement vers les basses 
fréquences et cela pour élargir, au maximum possible, la plage de fréquence dans laquelle le lieu 
déve [ÎT, (jo) | soit faible. Cette translation permet d'améliorer le filtrage des bruits de mesure 
et désensibiliser également le système corrigé dans le cas de présence des dynamiques non 
linéaires négligées lors de la phase de linéarisation du comportement réel du système physique à 


commander. 
e Choix de l’ordre entier n : 


L’augmentation de l’ordre entier n survenant dans la puissance du transfert complexe W,(s) 
permet d’assurer l’amélioration de la marge de robustesse sur RS [34]. Ce choix conduit à 
accentuer la pente du lieu de Ball 0) de l’ordre de 20 x n dB/décades en basse 
fréquence. En revanche, si ce paramètre s'accroit fortement et plus que nécessaire, la marge de 
sécurité du NP peut être considérablement dégradée, ce qui se détériore la condition de 
robustesse sur les performances nominales et s’affaiblit, par conséquent, les deux marges de 


robustesse RS/NP et RP du système bouclé. 
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2.1.2. Sélection du gabarit limitant le lieu de 6,4x[5,(j&)] : 


Comme nous l’avions signalé auparavant que, pour synthétiser un correcteur robuste, il 
faut tout d’abord choisir un gabarit fréquentiel dans lequel l’évolution du lieu de 6ax [S, (jo) | 
doit être limitée dans toute la plage de fréquence [34,55]. Ce choix est en général effectué sur les 
paramètres de la pondération irrationnelle ajustable W (s). Sachant que l’inverse de cette 
pondération fourni l’allure souhaitée du lieu de omux [S, Ga)|. La forme idéale de ce gabarit est 
donc choisie telle que la pente du lieu de ox [Se (jo) | soit inclinée autant que possible, et 
surtout en basse fréquence, puis elle s’acroite progressivement jusqu'à l'unité (1.e., 0OdB) en haute 
fréquence et cela en fonction de la croissance de la pulsation w. A cet effet, la Fig.2.3 représente 
l’allure désirée de ôux [S, (jo) | anisi que le gabarit imposée par le lieu de (Gnax[Ws(jæw)])"{. 


On obtient : 


Fig.2.3. Forme désirée de Omax [s,] comparée avec celle de (Gmax|Ws1) 71 


En général, le transfert complexe définissant la pondération irrationnelle ajustable W.(s) peut 


être s’exprimé par [34]: 


Ss m 
my +wB 
Ws(s) = (IE ) (2.5) 


Sachant que : 


e Ms E R'* : représente le gain maximal de OmyxlW,(jw)] obtenu en haute fréquence ; 
e €, E R'* : représente l’erreur de poursuite (1.e., gain minimal) fournie en basse fréquence ; 


e wp E R'* : représente la pulsation dans laquelle le gain de GmaxLWs(Gjw)] = —3dB ; 
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e me N'*: représente l’ordre entier de la pondération W/(s). 
Les paramètres méconnus de cette pondération peuvent être choisis comme suit : 
e Choix du gain maximal M,: 


En général, le gain maximal M, représente le pire cas à atteindre par le lieu de ox [S,Go)|. 


Il est s’exprimé par l’expression suivante: 


Ms — IlS,1L — max CALE ‘ @max) |) (2.6) 


@ER*+ 


De l’équation (2.6), Il est nécessaire de diminuer le paramètre M7 autant que possible, ce qui 
permet de tolérer l’apparition d’un pic de résonance admissible dans le lieu de omax [S,Go)]. 
En général, le choix de ce paramètre est souvent effectué de telle sorte que la marge du gain MG 
et celle de la phase MP du système en boucle ouverte doivent satisfaire les deux inégalités ci- 


dessous [34]: 


MG > LM. 
Ms-1 


1 (2.7) 
MP >2 arcsin() 


Dans la théorie de la commande robuste, le paramètre susmentionné est souvent choisi dans 


l'intervalle Ms € [1.5,2.75]. Ce choix permet d’apparaitre des pics plus ou moins admissibles 


dans le lieu de Onax Es Go). 
e Choix du gain minimal £,: 


Dans le cas idéal, du gain minimal &, est souvent choisi nul. Ceci permet d’avoir une erreur de 
poursuite assez faible notamment en régime permanant. En revanche, dans le but d’éviter le 
problème de division sur zéro lors du calcul du système augmenté, il est indispensable de choisi 


€, très faible devant le paramètre M. Ceci se traduit par l’inégalité suivante [59]: 
O&Ee «M, (2.8) 
e Choix de la pulsation w : 


En général, le choix de la pulsation w peut influer directement sur la bande passante du 
système corrigé. À cet effet, 1l est plus pratique de translater cette pulsation toujours vers les 
hautes fréquences de telle sorte qu’elle devrait être inférieure à la pulsation wr. Il convient de 
noter que l’augmentation progressive de la pulsation wZ peut assurer les avantages ci-dessous: 

Un rejet rapide de l’effet des incertitudes multiplicatives non structurées; 


Ÿ__ Une bonne dynamique de poursuite de la trajectoire de référence. 
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e Choix de l’ordre entier m: 


L’augmentation de l’ordre entier m survenant dans la puissance du transfert complexe de la 
pondération W,(s) permet d’assurer une bonne amélioration de la dynamique de rejection des 
perturbations [34]. Elle assure également d’une atténuation de l’ordre de 20 x m dB/décades 


dans la pente du lieu de Gux [S, (ja) | en basse fréquence. 


En revanche, plus m est s’accroit plus que nécessaire, la robustesse de SR se dégrade en haute 
fréquence ce qui peut affaiblir d’un côté, la marge du compromis de robustesse RS/NP et d’un 


autre côté, il peut fournir une violence de la condition sur RP du système bouclé. 
2.1.3. Sélection du gabarit limitant le lieu de ax [K (w) : S,(jw)] : 


D’après la Fig.2.1, il est clair que la commande fournie par le contrôleur H,, robustifié est 


obtenue comme suit [55]: 
u=KS:T-n-d)-T A4, (2.9) 


Dans la plupart des méthodes de synthèse, l’objectif est d’assurer une commande très 
raisonnable de point de vue coût de réalisation et effort de commande. L’objectif de la 
modération de la commande permet, d’un côté, d’éviter la saturation des actionneurs et 
provoquer leur destruction. D’un autre côté, il permet de limiter la bande passante du contrôleur 


et atténuer le gain d’énergie du signal de commande notamment en haute fréquence. 


Sachant que la réalisation de cet objectif est directement liée par l’amélioration de la marge de 
robustesse de la condition imposée sur la limitation de l’énergie de commande. C’est pour cette 


raison qu’on doit effectuer un bon choix sur la pondération irrationnelle ajustable Wy(s). 


Il convient à noter que l’amélioration requise de cette marge de robustesse nécessitant que le 
gabarit fréquentiel imposé par lieu de GmyxlWy(jw)]"? devient pratiquement similaire de celui 
imposé par le lieu de Gnax[Wr(jw)]"*. A cet effet, la forme idéale de cette pondération exige 
que le lieu de omax [K (Go) : S,Gw)| soit très faible en haute fréquence et stationnaire (plate ) 


autant que possible en basse fréquence. 


Par conséquent, la Fig.2.4 représente l’allure souhaitée de cette sensibilité qui est comparée 


avec celle imposée par le lieu de (Gmax[Wy(jw)])"!. On obtient : 
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Fig.2.4. Forme désirée de Omax [K | sl comparée avec celle de (omaxlWyl)"* 
En général, le transfert complexe définissant la pondération W,,(s) est donné par : 


V 


T7 
WC) | 2e (2.10) 


Sachant que : 

e My E R‘* : représente le gain maximal de Omaxl Wy(jw)] fourni à la pulsation © = wnin ; 
e £y € R'* : représente le gain minimal de Gmgxl Wy(jw)] fourni à la pulsation © = max ; 
e pc € R'* : représente la pulsation dans laquelle le gain de OmaxlWy(w)] = —3dB ; 
e ve N'* : représente l’ordre entier de la pondération Wy(s). 

Les paramètres méconnus de cette pondération peuvent être choisis comme suit : 

e Choix du gain maximal M,,: 
En tenant compte la relation suivante: 


Mu = Î|K.S, ff, = max (onax[K(&).S,(w)]) (2.11) 


@ER*+ 


Sachant que M, donne une information pertinente sur le coût admissible du signal de 
commande, il est donc aisé de le réduire autant que possible afin que le contrôleur robustifié 


puisse produire des commandes économiques d’un point de vue pratique. 
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e Choix du gain minimal £,: 


Le choix de £y = 0 permet de palier l’influence des bruits de mesure sur le système bouclé en 
haute fréquence. En revanche, pour des raisons pratiques, ce paramètre est choisi comme étant 
négligeable devant le gain maximal de W,,(s). Ce choix est traduit par la vérification de la 


contrainte ci-dessous : 
0 = Eu LC My (2.12) 
e Choix de la pulsation wc : 


En général, il est préférable de translater la pulsation wc Vers les basses fréquences tout en 
satisfaisant la contraintewgc < wc. Ceci qui permet d’assurer une stationnarité prolongée du 
lieu de Oynax [K (w).S, (w)] dans un large bande comportant les basses fréquences. Cela permet, 
d’un côté, de fournir une meilleure limitation de l’effet des bruits de mesure avant la 
pulsation w3.. Il permet, d’un autre côté de réduire l’effet des incertitudes de modélisation après 


cette pulsation. Par conséquent, on obtient un compromis de robustesse appelé « Roll — Off » 
e Choix de l’ordre entier v : 


L'augmentation de l’ordre entier v survenant dans le transfert de la pondération Wy(s) à 
pour but de synthétiser un contrôleur H, robustifié assurant des commandes plus réduites en 
comparaison avec celles fournies par le contrôleur H,, primaire. En revanche, l’augmentation de 
ce paramètre plus que nécessaire peut facilement détériorer le compromis de robustesse NP/RS 


ce qui implique une violence de la condition de robustesse sur RP. 
2.1.4. Sélection du gabarit limitant le lieu de onax[Sy(s)G(s)] : 


A partir de la Fig.2.1, il est clair de remarquer que la sensibilité S,(s)G (s) représente le 
transfert entre la perturbation d’entrée et la sortie du système bouclé [55]. II convient à noter que 
l’influence de cette perturbation devienne considérable notamment en basse fréquence. Dans le 
cas ou l’énergie de cette entrée est importante, il est indispensable de l’atténuer dans cette plage 
de fréquence afin de réduire autant que possible l’effet des bruits agissant sur la sortie du 
contrôleur robustifié. En général, la forme la plus simple utilisée pour cette pondération est 
souvent définie par un gain d’atténuation présélectionné par l’utilisateur. Par contre, dans 
certaines d’autres méthodes de synthèse, l’allure du gabarit fréquentiel décrivant l’inverse de 
cette pondération peut être schématisé suivant la Fig.2.5 dans laquelle l’allure de la sensibilité 
souhaitée est imposée de telles sortes que ses valeurs singulières maximales doivent être toujours 


se situes au-dessous de celles fournies par le gabarit suscité. On obtient donc : 
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es, bd 


Fig.2.5. Forme désirée de Omax [5 : G| comparée avec celle de (omax|Wscl)7* 


En général, le transfert irrationnel décrivant la pondération W,£(s) est donné comme suit : 


Wie(s) = | — (2.13) 


S+OSsG'VEsG 


Sachant que : 


Mse € R°* : est le gain maximal de Omax[Wsc(jw)] fourni en haute fréquence ; 


se € R'* : est le gain minimal de opux|Wsc(jw)] fourni en basse fréquence ; 


wse € R°* : est la pulsation fournie lorsque Onax[Wy(jw)] = —-3dB ; 


e LE N'* : est l’ordre entier de la pondération Ws(s). 
Les paramètres de cette pondération doivent être choisis de manière convenable sinon, la 
commande du système bouclé peut devenir très fluctuante. Le choix de cette pondération est 


généralement similaire que celui proposé pour la pondération S,,(s). 


2.2. Sélection automatique des pondérations irrationnelles ajustables: 


En général, le choix manuel des paramètres intervenant dans les transferts des quatre 
pondérations précédentes peut poser un problème sérieux [36, 60], en particulier lors du 
développement d’un modèle généralisé approprié qui assurant la synthèse du contrôleur 
robustifié avec lequel toutes les conditions à savoir NP, RS, NP/RS et RP sont satisfaites. Dans 


ce qui est suit, nous allons proposer une méthode assurant la sélection automatique des ces 
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pondération en utilisant l’algorithme génétique GA [61, 62, 63, 64,65]. Cet algorithme à pour 
but d’optimiser les paramètres de chaque pondération cité précédemment tout en assurant les 


formes souhaitées qui sont présentées par la Fig.2.6. On obtient : 


Omax [Ws 1] 


sasesceseresesese--ss 


l L 
Os V Esc os Msc 


Fig.2.6. Allures souhaités assurant la formulation du problème de la sélection automatique des pondérations 


2.2.1. Règles de réglage assurant la sélection automatique des pondérations souhaitées : 
Dans cette partie, nous allons citer quelques règles de réglages existant dans la littérature. Ces 
dernières sont inspirées à partir des quatre formes idéales présentées dans la Fig.2.6. D’où, on 
peut envisager cinq règles qu’à partir desquelles on peut définir l’espace de recherche du 
problème d’optimisation qui assure cette sélection automatique. On aura donc les règles 


suivantes : 
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e La désensibilisation du système bouclé aux effets des entrées exogènes nécessitant 
l’atténuation des maximas M , Mr, My et Msc. Ce qui produit une minimisation de pire cas 


(worst case en anglais) de chaque lieu de sensibilité. 


e L’amélioration de la marge de robustesse nécessitant l’atténuation des minimas &, , £r, Ey 
et &sc. En effet, la réduction de €, et £sG doit être assurée en basse fréquence, tendis que la 


réduction de £r et £y doit être assurée en haute fréquence. 


e L’amélioration souhaitée impose la translation des deux pulsations wZ et w$ç vers les 
hautes fréquences dans lesquelles l’effet des incertitudes de modélisation et celui des bruits sur la 
commande sont importants. La translation des deux pulsations précédentes à donner la possibilité 
d’accentuer autant que possible les deux pentes de Oux [S,] toux [S,6| respectivement et 


cela dans une large zone situe en basse fréquence. 


e L’amélioration souhaitée impose également la translation des deux pulsations wzr, ©gc 
vers les basses fréquences dans lesquelles l’influence des bruits de mesure sur le système bouclé 


ainsi que la réduction de l’énergie de commande sont importantes en régime permanant. 


e L’augmentation des ordres entiers à savoir: m , n , v et | est indispensable dans les 
transferts irrationnels WSs(s), Wr(s) , Wy(s) et WsG(s) respectivement. Cette étape doit être 
effectuée avec précaution, car l’augmentation de l’un des ordres dans une bande de fréquence 


peut influer sur l’autre ordre dans l’autre bande de fréquence. 


La Fig.2.7 et la Fig.2.8 représentent graphiquement les conditions de robustesse du système 
bouclé avant et après la robustification du contrôleur H,, primaire respectivement. On aura 


donc : 
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ll : 
D Dei CES DQL EE par 
CRERRPREN ET AE EEE 


CCE SEL REL RC ECL SC RCL PEACE SEES : 


rer Hs. : L à ; … Wukosn 
| À \NLyN geo pome em me 
#4: e Le 
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ut . UT Onax [Ws,Sy1] + Cinax [Wr Ty1] + Gimax [Wu K1Sy1 ] + Gmax [Wsc, GS] 
RE Et 
nr rinunl 
CAM ET 
4 titi) MARNE. 
LÉ Led y CSKA = 

| A Omax [Wu, K15,, | 

= 0mmmmni( Fr NLINl one mm ———————_— — 

£: N Omax [We GS, ] 


EE SE RES i Omax [Wr, Ty] 


Fréquence (rad/s) 


Fig.2.8 représentation des conditions de robustesse après la phase de la robustification proposée 


2.2.2. Formulation du problème d’optimisation sous contraintes : 


Dans la suit de ce chapitre, nous allons exploiter les cinq règles de réglage, citées 
précédemment, afin de déterminer l’espace de recherche dans lequel on cherchera à l’ensemble 


des paramètres optimaux construisant les pondérations irrationnelles souhaitées. Ce dernier est 
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présenté par des limites inférieures et supérieures limitant les solutions possibles du problème. 
Parmi ces solutions, on doit satisfaire la condition RP via la recherche à des pondérations 
adéquates fournissant, avec le modèle nominal, le développement du système augmenté 
approprié. Sachant que la formulation du problème de robustification s’articulera a priori sur la 
définition de la fonction objective ainsi que l’espace de faisabilité des solutions possibles 
assurant la sélection automatique des pondérations. Le paragraphe ci-dessous présentera les 


différentes contraintes imposées dans ce problème d’optimisation. 


Supposons que les pondérations irrationnelles fixes (1.e., pondérations initiales) citée du 


premier chapitre sont notées comme suit : 


e W,(s), en tenant compte ses paramètres : Ms, Es, ©, €t Mo (2.14) 
e W,,(s), en tenant compte ses paramètres : Mr,, Er,, ©pr, Et no (2:13) 
e _W;,(s), en tenant compte ses paramètres : My, Eu,, &gc, Et Vo (2.16) 
e Wic,(s), en tenant compte ses paramètres : Msc,, Est» Wsc, €t Lo (2.17) 


Supposons maintenant que les pondérations irrationnelles ajustables sont notées comme suit : 


e W;,(s), en tenant compte ses paramètres : Ms, Es,, ©p, et Mi (2.18) 
e Wr,(s), en tenant compte ses paramètres : Mr, Er,, @pr, Et n1 (2.19) 
e Wy,(s), en tenant compte ses paramètres : My, Eu,, ©Bc, et V1 (2.20) 
e Wsc,(s), en tenant compte ses paramètres : MsG,, &sc,; @sc, et li (2.21) 


D'où, l’espace de recherche du problème d’optimisation peut être construit, pour chaque 


pondération présentée par la Fig.2.6, par : 


> Pour la pondération Ws. (s), ses paramètres doivent être limités comme suit : 


Ms in — Ps” Ms, 
° ME [Ms Ms, : Sachant que Me MS (2.22) 


e &, € [ ] , Sachant que ESmax — So (2.23) 


0 <ps K1 


ESmin’ €Smax 
OBmin — Bo 

OB — Pos ‘ VB 
max op 0 (2.24) 


OBT) 


° WB, € [a in: Bar] , Sachant que Le » 
Pop PA) 
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Mmin — Mo 
eo ME [Min Mmax] ; Sachant que frs = Mo + Pm (2.25) 
Pm € N° 


> Pour la pondération W7. (s), ses paramètres doivent être limités comme suit : 


Morin — PMr° Mr, 
° Mr E [Mrs Mrpux] : sachant que Mr max = Mr, (2.26) 


ËTmin — Per‘ ÊTo 
° Er € [ër in: ET max] , Sachant que ÉTmax — ÊTo (2.27) 
0 < Per K1 


OBTmin — Por  VBETO 
OBr, € [war in OBT mx] , Sachant que OBTmax — ®BTO (2.28) 
0 < Poyr € 1 


Nmin — No 
M € [Mmin; max], Sachant que {max = No + Pn (2.29) 
Pn € N° 


> Pour la pondération Wy, (s), ses paramètres doivent être limités comme suit : 


Migin — PMy My, 
My, € [Mumnin: Mu max) , Sachant que Mois — My, (2.30) 


EUmin — Peu EUo 
e Ey, € [Euin: EU max) , Sachant que EUmax — EUo (2.31) 
0 < Er K 1 


OBCmin — Powgc  VBCO 


WBc, € [a zcin OBC max , Sachant que OBCmax — BCo (2.32) 
0 < Pope € 1 


Vmin — Vo 
V1 € [Unin; Vmax] ; Sachant que [me = Vo + Py (2.33) 
Po EN 


> Pour la pondération WG, (s), ses paramètres doivent être limités comme suit : 


MsG min — Ps ‘ MsG 


° Misc, € [Mai MSG max] , Sachant que Lee = Msc (2.34) 
0 < PMsc < 1 
ESGmin — Pesc ‘ ESGo 
e Esc, € [Escin EsGmuxl , Sachant que ÉSGmax — €SGo (2.35) 
0 < Pee 1 
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DSGmin — VSGo 


mn = Poce * © 
o We, € [oscmin Ssemal » Sachant que 4 *Fmax (2.36) 
1 < (RS re 
be = Lo 
e LE [lnin: max]; Sachant que {max = lo + Pi (2.37) 
1 E N* 


Nous rappelons ici que l’objectif de cette partie s’articulé sur l’amélioration de la robustesse 
RP tout en assurant la sélection automatique appropriée de toutes les pondérations irrationnelles 
ajustables citées précédemment. Pour cela, on doit commencer d’abord par la synthèse du 
contrôleur H, primaire à partir des pondérations irrationnelles fixes (initiales). On fait ressortir 
la marge initiale RP, Ensuite, on fait un appel à l’algorithme génétique pour optimiser les 
paramètres des pondérations irrationnelles ajustables souhaitées. Ces dernières sont ensuite 
utilisées avec le modèle nominal pour calculer le système augmenté correspondant avec duquel 
le problème de sensibilité pondérée est formulé par LFT,. Après, on fait un appel à la méthode 
H> pour résoudre ce problème, on obtient par conséquent, la représentation d’état du nouveau 
contrôleur H,. La nouvelle marge de robustesse, notée RP,, est ensuite calculer à partir de ce 
contrôleur puis la comparée avec RP,. Ces étapes sont réitérées à chaque itération k et 
l’algorithme s’arrête si la condition RP,,: < RP, n’est pas satisfaite. Il fournit, par conséquent, 
la représentation d’état du contrôleur H,, robustifié, noté ici par K,(s). Nous rappelons ici que 
notre objectif est de robustifier le contrôleur H,primaire afin d’améliorer leur marge de 
robustesse RP,. A cet effet, d’après le processus de robustification cité ci-dessus, le problème 
d'optimisation sous contraintes est définit comme suit: 

Omax[Ws, (Xe: ©). Sy, (j&)] + max [Wr, (Xe ©). Ty, (jw)] 


min,, [max + (2.38) 
œ : : : 
Omax | Wu, ( Xe) &).K ° S,(jw)| Go Omax [Wse, ( Xe) &).S, (jo) ° Gn(Gw)] 


En tenant compte la contrainte de saturation ci-dessous : 
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Pms' Ms, Ms, Ms, 
Pes ES Es, ES, 
OB OB, Pog VB 
Mo mi Mo ca Pm 
PMr' Mr, Mr, Mr, 
Per‘ ÊT ET, ÊTo 
Pogr ‘ VBTO &BT, OBT, 
de el a NET Pa (2.39) 
Puy Mu My, M | 
Peu ‘ EU EU, EU 
Pogc  VBCo BC; OBCo 
Vo V1 Vo + Py 
Pmsc * Ms Msc, M5Gs 
Pesc ‘ SG ËSG: ESGo 
OsGo OsG: Pose SCO 
Lo L Lo + 1 


ES 


Il convient à noter que le nombre de variable à optimiser du problème cité ci-dessus est 
considéré élevé. A cet effet, sa solution peut être attrapée (ou localisée) dans d'éventuels minima 
locaux existant dans l’espace de recherche. C’est pour cette raison qu’on doit résoudre le 
problème précédent par des algorithmes évolutionnaires on l’occurrence les algorithmes 


génétiques, l’algorithme d’essaim particulaire,.…., etc. 
2.2.3. Solution du problème par les algorithmes génétiques: 


Nous rappelons que le but de ce chapitre réside dans la recherche sur des pondérations 
ajustables adéquates dont ses structures sont décrites par des transferts d’ordre entier. A cet effet, 
seulement les algorithmes génétiques sont capables optimiser tous les paramètres de ces 


pondérations y compris ceux donnés avec ordre entier. 
2.2.3.1. Algorithmes génétiques : 


En général, les algorithmes génétiques consistent à faire évoluer un ensemble de points se 
trouvant dans l’espace de recherche, vers l’optimum du problème tout en commençant par un 
tirage initial aléatoire. L’idée des algorithmes génétiques consistent à générer une population de 
chromosomes représentant les paramètres à optimiser de chaque pondération irrationnelle 
souhaitée. La fonction objective (2.38), appelée Fitness, est ensuit évaluée dans chaque individu 
constituant cette population afin d’examiner le critère d’arrêt de cet algorithme. Le paragraphe 
suivant décrit les étapes à suivre de la robustification du contrôleur H,, primaire dans laquelle on 
introduit les algorithmes génétiques pour assurer la sélection automatique des pondérations 


irrationnelles adéquates assurant l’amélioration de la marge de robustesse RP,. Ces étapes sont 


résumées comme suit : 
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* Initialisation de l’algorithme de robustification : 


Dans cette étape, les transferts du modèle nominal G,(s) et les pondérations irrationnelles 
fixes Ws,(s), Wr,(s) , Wu, (s) et Wsc, (s) sont supposés connus. Sachant que, le transfert du 
contrôleur H, primaire K,(s) est assuré par la méthode H,, utilisant l’un des deux algorithmes 
cités en premier chapitre. Ce dernier permet de fournir la marge de robustesse initiale RP;. 
Finalement, les limites inférieure et supérieure de l’espace de recherche sont données à partir des 
règles de réglages précédente ainsi que le nombre des paramètres à optimiser est supposé égal 
à N. 

* Initialisation de l’algorithme génétique 


A l’itération k = 0, la population initiale X, est construite via la génération de N,, vecteurs 


 …. . 
aléatoires limités par l’intervalle [re indé us | 


+ Etape 1 : Codage de la population }, : 


En général, il existe plusieurs types de codage assurant la codification de chaque paramètre 
constituant la population X,. Parmi lesquels, on peut citer: le codage binaire, le codage réel 
(symbolique), codage sous forme d'arbre,..…., etc. Dans cette étape, le codage binaire est adopté 
pour coder chaque paramètre de X, en N, bits. D’où, N, chromosomes codés peuvent être 
générés dont la longueur de chacun est égale à N, X N. La concaténation de ces chromosomes 


binaires peut ressortir la structure générale de la population codée Y, ci-dessous: 


0,k 0,k 0,k 0,k 
Ve Ye, 2 Vo 
1,k 1,k 1,k 1,k 
Y, = Yes Ve Vcs Vca (2.40) 
(Np-1)k (Np-1)k (Np—1)k (Np-1)k 
Ci C2 C3 C4 


Sachant que le £°"€ chromosome codé fourni à l’itération k, est donné par: 
q P 


{chromosome}sx = {y5* ee 1 se") (2.41) 


Où £ = 0,1, , Lmax AVEC Pmax = Ny — 1. Sachant que chaque sous vecteur 7 du 49e 
chromosome fournissant les quatre paramètres codés d’une telle pondération irrationnelle 


sb : à _ | 
codée W,7", peut être donné, pour j = 1,:::,4, par : 
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Np bits Np bits Np bits Np bits 
 i — | [100...00] [100...01] 1[111...10] [100...11| | : Pour 120) 
UE np 
[110...00] 1[110...01] [110...10] [110...11| ” 
M — | 2k k êk Fk : Pour W* 
Mr To BTo ño (2.42) 
101...00] [101...01] [101...11] [101...11] e 
a = ( Ex FR Îk Pk : Pour We 
MU EUo ®BCo Vo’ À 
011...00 011...01 011...11 011...11 ” 
VE" — ( 1 Ek 1 Êk l 1 Êk , : Pour es 
MSGo EsGo WG Lo i 


* Etape 2 : Reproduction de la nouvelle population codée : 


Dans cette étape, une nouvelle population est assurée par deux opérations de base à savoir : le 


croisement génétique et la mutation génétique. 


a. Croisement : 


Le croissement permet de produire les chromosomes enfants 20 K héritant partiellement des 


caractéristiques des chromosomes parents existant dans la population Y,. Cela dépend d’un choix 
de la probabilité P. choisie a priori par l’utilisateur dans l’intervalle [0,1]. En général, il existe 
plusieurs méthodes de croissement. Parmi les méthodes les plus utilisées dans le croisement des 
chromosomes, on peut citer : la méthode de « croissement uniforme ». Nous verrons, dans la 
suite de ce travail, que le croisement d’un certain nombre des c individus est assuré à partir d’un 
test probabilistique supérieur à la valeur de choisie de P., ce qui permet de produire des 
nouveaux individus améliorant la convergence de l’algorithme génétique. La Fig.2.9 représente 


(-p),(k-1) (—-q),(K-1) ; 


un exemple de croisement des deux chromosomes parents y et Ye à l’itération 


£-p}k £-q)k : a 
k — 1 pour générer les deux chromosomes enfants 2 TPE et 2£ o respectivement à l’itération 


k. On obtient donc : 


Chromosome parent : 0 LD =[110...00] [110...01] 
ne De 
Chromosome parent : on p}(k-1 _ [110...00] [110...01] 
MED UD 
“à 1 
Chromosome enfant : Du = [110...00] .00] ] [110.01] 01] 
CUS q)k VOS q)k 
Chromosome enfant : 20 = [110.00] 00] ] [110...01] 01] 


PO p},(k—1) 


Fig.2.9. Exemples de croisement des deux chromosomes parents 


MD p),(k—1) 
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b. Mutation : 


Cette étape est indispensable lors de la génération de la nouvelle population Y,, ce qui permet 
de générer les nouveaux chromosomes he dont certains ses éléments sont totalement différents 


à ceux existant pendant tous le processus d’optimisation. En général, la mutation des 

chromosomes dépende de la probabilité P,, qui est souvent choisie préalablement par l’utilisateur 

dans l’intervalle [0,1]. Cette étape permettant: 

e d'atteindre une propriété appelée érgodicité avec laquelle la solution du problème peut être 
convergée vers un minimum global. 

e de garder une diversité dans l’évolution des individus, ce qui est primordial pour les 
algorithmes génétiques. 


e de limiter le risque de la convergence prématurée des solutions. 


Dans notre travail, la mutation de certains chromosomes est assurée par l’introduction de 


certains nouveaux paramètres qui ne figurent plus dans toutes les populations existantes. Ces 


paramètres doivent être choisis aléatoirement dans l’espace de recherche | | La 


Xcmin Xcmax 
: - : £,k — 
Fig.2.10 représente un exemple de mutation du chromosome y; pour générer le nouveau 


chromosome LA On obtient : 


8,k [o11...00]'[011...01]:[011...11] [011...11] 


Yes = 
ÊK Ù BK | mn LR 


MiGo L sGo_ L d SGo 


Sélectionné à partir 
de l’espace de recherche 


[re min ,XCmax ] 


h£k = [011.00] 01) [011...11] [011...11] 
mt » , 


SGo 


Fig.2.10. Exemples de mutation du chromosome yÉ “4 


+ 


+ Etape 3 : Evaluation de la fonction objective 
A partir des deux opérateurs cités ci-dessus, les nouveaux chromosomes codés peuvent 
| , ; dk : 0 £k ; , 
contenir les sous ensembles codés, on l’occurrence Ve: » he; etzc; . Pour évaluer la fonction 


objective, on doit décodés tous les sous-ensembles précédents en N, individus décimales. D’où, 


la structure générale de cette nouvelle population peut être présentée comme suit: 
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0 0 0 
Xe: Xc2 XC3 Xc, 
1 1 di d 
Ze Ze Ze Ze 
Xk _ _ : 3 4 (2.43) 
Np=1 yNp=1  Np-1 Np-1 
p p p p 
Xe, Re, C3 C4 
Sachant que le £°"€ individu, fourni à l’itération k , est défini par: 
à x 4k  4k 
{individu} = {xer do + (2.44) 


: à £k :  _ 
Où chaque sous ensemble réels Xç;_ contenant les paramètres de la pondération irrationnelle 


W£ Kest défini par : 
Lk L,k Lk £,k mé k 
x = (MS € op : Pour W®* 


. ex 
T, ET,  WBro 10 ) : Pour W, 

êk êk (2.45) 
Un EU BC, Vo ) : Pour, 

LK\. £,k 
SGo E5Go SG Lo ) : Pour Wc 
Les pondérations 20) : A K WEEK et WE pour £ = 0,1,2,-+, £max = Ny — 1 Sont ensuite 
utilisés pour augmenter le modèle nominal Gy, on fait donc développé N,, systèmes augmentés et 
N, problémes de sensibilité pondérées. La résolution de chaque problème par la méthode H 


’ : : = — 2 ÈR ni 
fournie la représentation d’état du contrôleur H,, primaire correspondant, noté K,'"(s) ainsi que 


caque système bouclé correspondant permet de calculer les quatre sensibilités . 1 
£k c£k £k 
K5 Sy et GnSy, : 


A partir de ces dernières, on obtient N,, marges de robustesse sur les performances robustes 


Celles-ci sont obtenues, pour £ = 0,1,2, +++, N,, — 1, par l’expression suivante : 


BR = ne SA, + + RSA, + I GS, 2.46 


* Etape 4 : Sélection des meilleurs individus: 


Cette étape est s’inspirée du principe de Darwinien lors de la sélection naturelle des 
individus. D'où, les individus fournissant les meilleures améliorations de RP (ï.e., décroissance 


par rapport à RP,) sont maintenus, tendis que les individus restants sont totalement exclus dans 
les étapes de l’algorithme génétique. 
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* Etape 5 : Test d’arrêt 


Le critère d’arrêt est associé par la plus faible marge de robustesse, notée RP,, sachant qu’elle 


est donnée par l’expression ci-dessous : 
RPmin = Min{ RPÿ*}, pour £ = 0,1,2,+,N, — 1 (2.47) 


S1 la condition RPmin < RP, n’est vérifiée, le transfert du contrôle H,, robustifié est celui 
donné par K,(s). Sinon, on doit faire les changements ci-dessous avant de retourner à l’étape 1 
par la population qui contient seulement les meilleurs individus sélectionnés dans l’étape 4. Ces 


changements sont les suivants : 
e Incrémenter le nombre d’itération k par 1,ie k<k+1; 
e Mettre à jour la valeur de RP, par celle de RPmin , te. RPo <— RPin ; 


e Mettre à jour le transfert du contrôleur KÇ(s) par celui fournissant la marge de robustesse 


RP. min 
2.3.  Inconvénient de la robustification basée les pondérations irrationnelles ajustables: 


Le but de robustifier le contrôleur H, primaire cité en premier chapitre réside dans 
l’amélioration de sa marge de robustesse sur RP. A cet effet, on doit parcourir l’algorithme 
génétique pour exploiter l’espace de recherche, basé sur les règles de réglage précédentes, dans 
lequel les paramètres optimaux des pondérations irrationnelles peuvent être optimisés. 
Cependant, 1l existe souvent certaines difficultés lors de la sélection de ces pondérations. Celles- 
ci peuvent être causées par la saturation de certains paramètres dans ses limites inférieure ou 
supérieure, ce qui peut ralentir la convergence de l'algorithme et même enfreindre la condition 
RP. Pour remédier à cet inconvénient, nous devons donc élargir légèrement et avec précaution 
les limites là où elles sont saturées. Toutefois, cette alternative est valable pour certains 


paramètres à optimiser. 


Néanmoins, pour un ordre entier d’une telle pondération supposé saturé dans sa limite 
supérieure, sachant que l’accroissement de cette limite peut influer sur l’amélioration d’une autre 
marge de robustesse supposée également liée par une pondération contenant un ordre saturé 
auparavant dans sa limite supérieure correspondante. Dans ces conditions, la robustification 
proposée s’avère inefficace à cause de la recherche restreinte dans l’espace définissant les 
pondérations irrationnelles. Nous verrons dans le chapitre suivant que la recherche des solutions 


dans un espace comportant des pondérations fractionnaires ajustables peut améliorer la 
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convergence de l’algorithme d’optimisation, ce qui assure une bonne amélioration de la 


margeRP. 


2.4. Simulations et interprétations des résultats : 


Dans cette partie, le contrôleur H,, primaire est robustifié de telle sorte que sa marge initiale 
de robustesse RP, peut être améliorée via la robustification utilisant l’algorithme génétique. 
Sachant que cette amélioration est assurée par la sélection automatique des pondérations 
irrationnelles ajustables dont les limites inférieure et supérieure, contraignant ses paramètres, 
sont déterminées à partir des règles de réglage citées précédemment. A cet effet, nous devons 
reformuler les transferts des trois pondérations fixes citées en premier chapitres de telles sortes 


quelles devenues comme celles présentées par les équations (2.1), (2.5) et (2.10). On obtient 


donc : 
S 
—+0.06 
| We (s) = —25 
Ws (s) — 0.4(s+0.15) 5, ( ) Ss+0.06:10-6 
0 S+6-10—° s 4 1 
4.2482:(5+0.4218 — 09158 2.1713 
Wr (s) = 2276102418) Wr, (s) 7 02354 (1.48) 
0 s+3.891 0.9158 
15.152:-(5+0.5883) S 4 
Wu, S) = ———— 20528 
0 Ss+50.4 Wu (s) = DE 


2.05 


De l'équation (1.48), les paramètres constituant les transferts des trois pondérations 


irrationnelles fixes, à savoir Ws (s), Wr (s) et Wy, (s), sont obtenus comme suit : 


Ms =25 (Mr =21713 (My, = 4.82 

O8, = 0.06 OBT, = 0.9158 ” OBCo = 2.05 

es, = 0.01 *| er, = 0.2354 | ey, = 0.066 
ns, = 1 Nr, = 1 ny, = 1 


(1.49) 


Selon les règles de réglage citées précédemment, pour l’amélioration de la marge RP, 
nécessite donc l’amélioration simultanée des marges de robustesse NP,, RS, et le compromis 
NP:/RSo. À cet effet, on désire d’abord affaiblir le gain maximal du lieu de oux [Ws, (ju) | en 
haute fréquence et accroitre leur ordre entier pour assurer une atténuation plus descente du lieu 
en basse fréquence. De plus, le paramètre &ç, est maintenu égal à &s,, tendis que la pulsation 
w8.,est choisie libre dans l'intervalle] 0, &gr, [. Ce dernier choix permet de réduire la complexité 
du problème d’optimisation. L’allure souhaitée de S,, (s) peut être assurée par le choix des 


contraintes ci-dessous : 
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OMs " Ms, < Ms, < Ms, 2 < Ms, < 2.5 
0 < &p, < pr, , avec: 40 < wp, < 0.9158 (1.50) 
Mo < M < Mo + Pm 1<m <3 


Il est à noter que le cas prévu par cette sélection automatique présenté est de fournir les 
paramètres Ms, = 2.5 et m; = 1. Dans ce cas, on dit que la robustification a conservé la même 
norme H,, et la même atténuation des valeurs singulières maximales fournies par la sensibilité 
initiale du contrôleur H,, primaire. Par contre, l’amélioration de RP, peut être aboutie par 
l’amélioration de la marge de robustesse RS,. Ensuite, on désire également concentrer beaucoup 
plus sur l’amélioration de la marge de robustesse RS,, car la stabilité du système bouclé est 
indispensable dans la plupart des applications industrielles, en particulier le cas de présence des 
incertitudes de modélisation. Pour atteindre cet objectif, les différentes valeurs à savoir des deux 
gains maximal et minimal et la pulsation du lieu de Gux [Wr, (ja) | doivent être affaiblies autant 
que possible. Par ailleurs, l’ordre entier de la pondération doit être augmenter afin d’assurer une 
bonne minimisation de Oynax Fe (jw)| en haute fréquence. Dans ce cas, l’allure souhaitée de 


T., (s) peut être assurée par le choix ci-dessous : 


Pur * Mr, < Mr, < Mr, 1.75 < Mr, < 2.1713 
. < _ ee < 0. 
Per ET S< ËT, S< ET avec: 10 S< ET, S< 0 2354 (1.51) 
Por ‘ BTO < OBT, < OBT) 0.5 < OBT, < 0.9158 
No < M < No + Pn 1<nm <3 


La pondération irrationnelle ajustable Wy, (jw) est supposée fixe dont le transfert irrationnel 
est supposé égal à celui de la pondération Wy, (jw). Finalement, les paramètres de réglage de 


l’algorithme génétique sont choisis comme suit : 


e Type de codage du chromosome : Binaire; 

e Nombre de bits utilisé dans le codage binaire : 16; 
e Nombre de paramètres à optimiser : N =7; 

e Méthode de croisement : Croisement uniforme; 

e Taux de croisement : P. = 70%; 

e Taux de mutation: P, = 45%; 


e Nombre de génération: 30. 


L’algorithme génétique est compilé 20 fois avec différentes populations initiales générées 
dans l’espace de recherche limité par les contraintes citées précédemment. D’où, la Fig.2.9 
représente l’évolution de la fonction objective pendant les itérations de l’algorithme de 


robustification proposée. On obtient : 
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Fonction objective 


0 3 6 9 12 15. 18 21 24 
Générations 


27 30 


Fig.2.11. Evolution de la fonction objective pendant les itérations de l'algorithme de robustification proposée 


La Fig.2.11 montre clairement que la marge de robustesse RP n’atteinte 


pas une valeur 


inférieure à l’unité. Cela se traduit, plus tard, par le fait que la robustification proposée basée sur 


l’algorithme génétique devienne incapable de réaliser l’amélioration souhaitée de cette marge de 


robustesse. La principale question qui peut être posée maintenant est celle de savoir à quelle 


limite fréquentielle la robustification proposée peut améliorer cette marge de robustesse? 


ES 


Pour répondre à cette question, nous devons donner le transfert de chaque pondération 


irrationnelle fournie lorsque l'algorithme génétique atteint cette limite d’amélioration de cette 


marge de robustesse. De plus, nous devons donner le transfert du contrôleur H,, robustifié fourni 


lorsque la méthode H., résout le problème de sensibilité pondérée correspondante. On a donc : 


2 < Ms, = 1.72 < 25 
2 0 < wp, = 0.02 < 0.9158 
) 1<m=2<3 
42 avec: 4175 < Mr, = 2.07 < 2.1713 
Wr, (s) = =) 1074 < er, = 1073 < 0.2354 
0.5 < wgr, = 0.868 < 0.9158 
1<m=2<3 


S 
Fe - +002 
S1 5+0.02-V10 3 


1073 
s+1 


0.868 


(1.52) 


D'où, le Tab.2.1 récapitule le taux de réussite de l’algorithme génétique pour décroitre ou 


d’accroitre les paramètres cités précédemment. On obtient : 
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Tab.2.1. Taux de réussite de l'algorithme génétique 


Paramètre avant | Paramètre après | Taux de réussite 

la robustification | la robustification | de décroissement 
Ms 2.5 1.72 31.20% 
My 2.1713 2.07 4.66% 
Er 0.2354 107$ 99,57% 
Du 0.9158 0.868 5.22% 

Taux de réussite 

d’accroissement 
m 1 2 100% 
n 1 2 100% 


De plus, la méthode H,, basée sur l’approche LMIS fournie un contrôleur H, robustifié 
d’ordre 11 et un niveau de performance égal à y-0.9107. Le transfert de ce contrôleur est donné 


par : 


K, = 37.3-(5+50.46)(5+28.66)(5+26.37)(s+1.4)(5+0.99)(s+0.18)(s+0.016)(s+0.0102)(s2+1.582-5+1.195) 
17 (s+1017)(5+39.92)(s+18.44)(s+0.51)(s+0.019)(s+66xX10 5)(s+59xX10-5)(52+0.475+0.88)(52+3.765+5.13) 


(1.53) 


La condition de robustesse sur RS fournie après la robustification proposée est illustrée par la 


Fig.2.12. On obtient: 


Condition sur la Stabilité Robuste 
10 T T T T 


-20 F 


-30 Fr 


Valeurs Singulières Maximales (dB) 


-50 F 


-60 L L L 1 
107$ 10? 107! 10° 10° 
Pulsation (rad/s) 


Fig.2.12. Condition sur SR avant et après la robustification du contrôleur H,, primaire 


D’après la Fig.2.12, on peut remarquer que les valeurs singulières maximales de la nouvelle 
sensibilité complémentaire sont bien atténuées haute fréquence notamment dans la plage 
fréquentielle we[0.0532, wmax] radians/secondes. De plus, le lieu de cette sensibilité possède 
une pente plus descente que celle fournie par la sensibilité complémentaire initiale. Cette 
atténuation peut conduisant à une amélioration de la marge de robustesse sur les performances 


robustes du système bouclé. 
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Ensuite, la condition de robustesse sur NP fournie après la robustification proposée est 


illustrée par la Fig.2.13. On obtient: 


Condition sur les Performances Nominales 
10 I I T T 


Valeurs Singulières Maximales (dB) 


-60 L 1 L L 
107$ 102 107 10° 10! 
Pulsation (rad/s) 


Fig.2.13. Condition sur NP avant et après la robustification du contrôleur H, primaire 


D’après la Fig.2.13, il est claire de remarquer l’amélioration de la marge de robustesse NP par 
le contrôleur H, robustifié, en particulier dans la plage fréquentielle we[0.01, wux] 
radians/secondes. Celle-ci peut conduire aussi à une amélioration de la marge de robustesse sur 
les performances robustes du système bouclé. Maintenant, le compromis de robustesse fourni par 
le contrôleur H,, robustifié est comparé par celui issu du contrôleur H,, primaire. Cette 
comparaison est illustrée par la Fig.2.14. En revanche, la condition de robustesse sur les 
performances robustes avant et après la robustification du contrôleur H,, primaire est illustrée 


par la Fig.2.15. On obtient : 


Compromis de robustesse fourni avant at après la robustification proposée 
-0.45 T RES EE | T ML ET T LE A ES EE T CE EEE | T 


— Compromis fourni par le contrôleur primaire 
— Compromis fourni par le contréleur robustifié 


Valeurs Singulières Maximales (dB) 
o 
D: 
T 
1 


-0.95 i eh où di | i je À —$ ji dl i ii jo ji ji | i De à ÿ hf 
107 1072 10°! 10° 10! 
Pulsation (rad/s) 


Fig.2.14. Compromis de robustesse fourni avant et après la robustification du contrôleur H, primaire 
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Condition sur les Performances Robuste 
2 T T T T 


——— RP fournie par le contrôleur primaire 


——— RP fournie par le contrôleur robustifié 


Valeurs Singulières Maximales (dB) 
œ 
T 


-18 1 ji 1 1 it L 1 
107$ 1072 1071 10° 10! 
Pulsation (rad/s) 


Fig.2.15. condition de RP fournie avant et après la robustification du contrôleur H,, primaire 


D’après la Fig.2.14, on peut remarquer que la nouvelle marge du compromis de robustesse est 
assurée avec une marge de sécurité très élevée en comparaison avec celle fournie par le 
contrôleur H,, primaire. Il est à noter que la robustification proposée à la capacité d'augmenter 
les marges de sécurité des trois conditions de robustesse citées précédemment. Elle est capable 
aussi d’améliorer la marge initiale de RP dans toutes les fréquences limitées par 
DE[ min Wmax | (Voir la Fig.2.15), ce qui implique que les performances du système sont 
assurées par le contrôleur H,, robustifié pour un nombre des régimes perturbés plus large que 
celui du contrôleur H,, primaire. Néanmoins, le problème de garantir une marge RP inférieure à 
l’unité dans toute la plage de fréquence reste toujours persistant. Ce problème est dû à la 
violence de cette condition en basse fréquence (régime permanant), en particulier la plage de 
fréquence we[107%,0.009 ]radians/secondes. Il convient à noter que la principale cause de la 
violence de cette condition réside dans le type des pondérations irrationnelles utilisées lors de la 
robustification proposée. Nous verrons, dans la suite de ce travail, que l’introduction d’un 
nouveau type pondérations dans le processus de la sélection automatique assurée par 
l’algorithme évolutionnaire peut conduire à satisfaire la condition de RP dans toute la plage de 
fréquence. Ce nouveau type sera basé sur l’introduction de la puissance fractale dans les 
transferts des pondérations, on fait ressortir les pondérations fractionnaires ajustables. La 


sélection automatique de ces dernières faisant l’objet du troisième chapitre. 


Maintenant, le contrôleur robustifié par cet algorithme proposé est utilisé pour commander les 
192 régimes perturbés précédents en utilisant les mêmes entrées exogènes citées en premier 
chapitre. La Fig.2.16 représentes les sorties globales du système bouclé, tendis que la Fig.2.17 


représente la partie transitoire correspondante. On obtient donc : 
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Fig.2.16. Sorties fournies par le contrôleur H, robustifié pour 192 régimes perturbés 


D’après la Fig.2.16, on peut remarquer que le contrôleur A, robustifié à la capacité de fournir 
des meilleures sorties que celles fournies par le contrôleur H,, primaire. Ceci se traduit par 


l’atténuation des erreurs de poursuite £59%, de tous les régimes perturbés à partir de l’instant t = 


80 à , , Lx : D 
_ secondes. De plus, La réponse fournie par le pire des régimes perturbés est caractérisée par le 


e 


2 2 N 29.33 x z ; : 
temps de montée égal à tn = secondes. C’est une réponse néanmoins moins lente en 


e 


comparaison avec celle fournie par le contrôleur H,, primaire. De plus, elle est caractérisée par le 


temps de réponse égal à t 


80 : Rs : 
rs — 7, Secondes, ce qui assure un taux d’amélioration égal à 
e 


66.66%. Le dépassement maximal de cette réponse est égal à Dnux = 05%, ce qui assure un 


taux d’amélioration égal à 58.33% (voir la Fig.2.17). 


Fig. 2.17. Agrandissement de la partie transitoire de la dynamique de poursuite fournie par le 


contrôleur H,, primaire pour 192 régimes perturbés 
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La Fig.2.18 représente les commandes assurant les sorties globales citées précédemment. On 


obtient : 


Fig. 2.18. Commandes fournies par le contrôleur H., robustifié pour 192 régimes perturbés. 


D’après la Fig.2.18, on peut remarquer que toutes les commandes fournies par le 
contrôleur H,, robustifié sont moins couteuses que celles fournies le contrôleur H, primaire, 1.e., 


Umnax = 4.50 en régime transitoire. 
2.5. Conclusion : 


Dans ce chapitre, nous avons proposé la robustification du contrôleur H,, primaire dans 
laquelle l’algorithme génétique a assuré la sélection automatique des pondérations irrationnelles 
ajustables, ce qui peut résoudre le problème généré par le choix des paramètres de chaque 
pondération de façon manuelle. Le problème de la sélection des paramètres optimaux est 
reformulé, à partir de quelques règles de réglage, sous forme d’un problème d’optimisation sous 
contraintes dans lequel la fonction objectif à minimiser représente la condition RP ainsi que et 
les contraintes à satisfaire représentent les limites inférieures et supérieures imposées par les 
allures souhaitées de différentes sensibilités fournies par le système bouclé par le contrôleur 
robuste à synthétiser. Les étapes requises de cette robustification proposée ont été détaillées ainsi 
que la validée de cette proposition a été montrée sur le même système hydraulique cité en 
premier chapitre. Les résultats de simulation ont confirmé que l’inefficacité d’utiliser les 
pondérations irrationnelles ajustables dans l’étape de la sélection automatique des pondérations 
adéquates satisfaisant la condition RP avec une marge de sécurité acceptable. L'utilisation des 
pondérations fractionnaires ajustables ainsi que la détermination de ses paramètres optimaux 


faisant l’objet du troisième chapitre. 
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Chapitre 3 


Commande robuste utilisant les pondérations 
rationnelles ajustables 


3. Introduction : 


Ce chapitre à pour but d’introduire les pondérations fractionnaires ajustables dans le critère de 
sensibilité pondérée. Les paramètres de ces dernières peuvent être assurés par la résolution d’un 
problème d’optimisation sous contraintes ce qui permet d’aboutir systématiquement les 
pondérations souhaitées à travers de satisfaire les mêmes règles de réglage citées en deuxième en 
chapitre. L’objectif est donc d’assurer la condition sur RP avec une marge de sécurité contre les 


incertitudes non structurées de modélisation. 


Rappelons ici qu’on a confirmé par simulation que l’utilisation des pondérations irrationnelles 
ajustables dans le problème de robustification du contrôleur H,, primaire s’avère inefficace 
notamment pour atténuer simultanément de façon optimale les deux pentes des deux sensibilités 
directe et complémentaire du système corrigé. À cet effet, nous allons proposer une deuxième 
alternative permettant de résoudre le problème de la sélection convenable des pondérations 
fractionnaires adéquates assurant l’amélioration de la marge de robustesse RP non seulement 
pour le régime de fonctionnement nominal du système physique, mais aussi pour tous les 
régimes perturbés. Cet objectif sera atteint par l’utilisation de l’algorithme de PSO ce qui 


représente la contribution principale de ce travail. 


Sachant que l’algorithme génétique offre la possibilité d'optimiser les variables réelles ainsi 
que les ordres entiers survenant sur les puissances des pondérations irrationnelles ajustables 
citées en deuxième chapitre. De ce fait, l’amélioration de la marge de RP par l’algorithme de 
PSO n'enlève rien de l'efficacité de l'algorithme génétique. Car le type des pondérations qu’on a 
choisit dans l’algorithme de robustification cité dans le deuxième chapitre qui restreint 
l’efficacité de cet algorithme. Sachant que la méthode de synthèse H,, fonctionnant seulement 
par les pondérations d’ordre entier. C’est pour cette raison qu’on doit a priori appliquer une 
approche d’approximation assurant le rapprochement des transferts rationnelles, en module et en 
phase, par des filtres équivalents d’ordre entier. Dans notre travail, la méthode adoptée pour cette 


étape est celle basée sur l’approche d’identification fréquentielle. 


Finalement, ce chapitre sera achevé par la partie de simulation dans laquelle les 


performances fournies par le contrôleur Æ,, robustifié synthétisé par cette deuxième stratégie de 
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robustification seront comparées, dans les plans fréquentiel et temporel, avec celles fournies par 
le contrôleur H, robustifié synthétisé par la première stratégie de robustification citée en 


deuxième chapitre. 
3.1. Motivation : 


En général, les équations différentielles ordinaires est un outil mathématique très utilisé, dans 
les sciences de l’ingénieur, pour bien décrire l’évolution du comportement réel d’un système 
physique. Récemment, les scientifiques ont généralisé l'intégration et la dérivée d’ordre entier de 
ces équations par des valeurs réelles dont des comportements fractionnaires peuvent être prises 
en considérations lors de la modélisation des systèmes physiques. De ce fait, le concept 


fractionnaire peut être illustré par la Fig.3.1 comme suit : 


3 ü 5 + == à2 


NN NN RTS 
a © 


Dérivations 


Fig.3.1 : Concept de dérivations et intégrations dans le calcul fractionnaire 


La Fig. 3.1 présente les dérivations et les intégrations successives de la fonction f(t) dont des 
intervalles d’ordre entier. La suite de ce travail sera basée sur les ordres fractionnaires possibles 


qui existent entre deux dérivées ou intégrations successives d’ordre entier. 


3.2. Eléments de la théorie du calcul fractionnaire : 
3.2.1. Fonctions usuelles utilisées dans la dérivation d’ordre non entier : 


En général, 1l existe quatre fonctions usuelles fournissant la définition de la dérivée fractionnaire. 
Ces dernières peuvent être illustrées comme suit : 

3.2.1.1. Fonction Gamma d'Euler : 

En mathématiques, cette fonction gamma est souvent notée par le symbole l'(z). C’est une 
fonction monotone et strictement décroissante dans le plan complexe l'(z) € C. Elle est définie 
comme suit [67] : 


r(0) =4 


T(z) = j tG-De-tdt Avec: os 2: 


(3.1) 


Sachant que cette fonction est définie, dans l’espace des variables d’ordre entier, par : 
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T(n+1)= n!,avec Vn € N* (3.2) 
En tenant compte l’expression ci-dessous: 

TGS ER +z Lure e ‘dt = #1(&) (3.3) 
D'où, l’expression (3.2) est encore valable pour la variable z € C , avec : 

F(z+1)= 27T() (3.4) 


3.2.1.2. Fonction Beta : 


Cette fonction est définie par l’intégrale eulérienne de premier ordre, donné par [67]: 


BE,9) = PRG -x)P xt dx, p,q € R4 (3.5) 


En général, la fonction Beta B(x,*) peut être exprimée en fonction de l'(*) comme suit : 


r(p)r 
BG,9) = TE = B(a,p), p,qaeR, (3.6) 


3.2.1.3. Fonction Mittag — Leffler (ML) : 
La fonction Mittag-leffler de deux variables à et ff est, en général, définie dans le plan 
complexe par [67]: 


e Pour (a,B) > 0 


k 
Fa Z 
Ep) = 2-0 ra E) (37) 
De l’équation (3.7), on peut envisager les deux cas particuliers suivants : 


e Pour f = 1,on obtient: 


Zzk 


Ea1C@) = Ei(z) = Dk=0 rer D (3.8) 


e Pour «a = 1, on obtient : 


AG eNe Se 2. (3.9) 
LIN Ak=0 rt) ÆK=0g 


3.2.1.4. Fonction Miller — Ross (MR): 


En mathématiques, la fonction de Miller-Ross est définie comme suit [67] : 


akzk+v 


&,(v, a) = DE=0 Fur D (3.10) 
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Elle est encore exprimée, en fonction de la fonction Mittag — Lef fler, par l’expression (3.11) 


suivante : 
£&,(v, a) = Z°E;,v+1 (42) Gin 


3.2.2. Fonctions usuelles utilisées dans l’intégration d’ordre non entier : 

En général, il existe trois définitions de l’intégration fractionnaire. Ces définitions sont 
données comme suit: 

3.2.2.1. Définition au sens de Riemann -— Liouville (RL) : 

Cette définition est inspirée de la définition de la n°” intégrale de la fonction f(t), qui est 


donnée par le théorème de Cauchy comme suit [67, 68] : 


t(t-T)" 1 


OL (3.12) 


De l’équation (3.12), l’intégral du premier ordre 14 f(t) de la fonction f(x) sur l’intervalle 


[O +] est exprimé par la relation (3.13) comme suit : 
of = JFG dx (3.13) 


Sachant que l’intégral du second ordre ,12f(t) de la même fonction précédente sur 


l'intervalle [O +] est défini par : 
OC = [, Je FR) dr dx G.14) 
D'où, l’équation (3.14) peut être réécrite, suivant l’égalité de Dirichlet, comme suit : 
OS), aar (3.15) 


Il est claire que la fonction f (t) est constante par rapport à x. D’où, il est possible de réécrire 


l’équation (3.15) comme suit: 


O8fE) = SG -T)f( dr. (3.16) 
De l'équation (3.13) et l’équation (3.15), on peut généraliser la définition de la 
n° intégration de la fonction f(x) par l’expression donnée par (3.17) comme suit : 
t t(t-T) 1 
fs). LiOdr.-dr=) men (3.17) 


En 1847 Riemann a généralisé le théorème de Cauchy sur l’intégration d’ordre non entier 


comme suit : Pour un ordre fractionnaire « € R* compris entren—1letn (ien—-1<a<n 
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avec n € N°), l’intégration d’ordre non entier de la fonction f(t) sur l’intervalle t, et t peut être 


s’exprimée par : 
RIELFON = (LE - 91. far) (3.18) 


D'où, la dérivée fractionnaire correspondante peut être donnée, au sens de Riemann-Liouville, 


par : 
DEN = ee (JE - n-4"2. f(Hdr) (3.19) 


3.2.2.2. Définition au sens de Caputo : 


La dérivée d’ordre non entier à compris entre n — 1 et n de la fonction f(t) est exprimée, au 


sens de Caputo, comme suit [66, 67] : 


Eero= ref (O)]= 0 oppose pe dr (320) 


Supposons que D” est l’operateur dérivé d’ordre n avec n € N*. D’où, les propriétés ci- 


dessous sont toujours vérifiées : 
e D".D =] mais D ".D' + 
° D.D'F(E) = FD - DEF (0) 
Sachant que f“(t) représente la a*”Ædérivée de la fonction f(t). 
3.2.2.3. Définition au sens de Grunwald — Letnikov (GL) : 
En général, la dérivée première f (1) (#) de la fonction f(t) est exprimée par [66, 67] : 


Sachant que la dérivée seconde f 2)(t) est définie par l’expression (3.22) comme suit: 
@) 9 = Jim OT ER) 
FO = lim OT G22) 


Cette dernière expression peut être réécrite comme suit : 
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C2) 4 = im 1 @-fG-h) _ fŒ-h)-f(t-2h) 
F0 = time [LE] (3.23) 
D'où, la formule simplifiée de l’expression (3.23) et obtenue comme suit : 


, ()-2f(t-h)+f(t-2h 
fA(t) = di PI (3.24) 


De façon similaire, la dérivée troisième f ©? (£) de la fonction f(t) est définie par : 
C)(b-fC)(t- 
fa) = dim IE (3.25) 


En faisant quelques manipulations mathématiques sur l’expression (3.25), on obtient : 


FE) — lim PE EE (3.26) 


D'où, la n°"® dérivée f ®)(t) est donc exprimée par l’expression (3.27) comme suit : 
à n 
FE) = lim Ereo( 0)" (0) FC rh) (3.27) 


Sachant que la combinaison (”) est définie par l’expression (2.28) comme suit: 


ds _ nin-1)(n-2)-(n-r+1) 


. _ (3.28) 


A partir de ces dérivées d’ordre entier, la dérivée fractionnaire d’ordre a de la fonction f(t) 
continue et intégrable sur l’intervalle [t,, + est définie, au sens de Grunwald — Letnikov, 


par l’expression suivante : 


[ete 


GDELF ON = lim, frs Fe) (e-D C)-fœ-5 1)} (329) 


— 2 : . 2 t-—t 
Sachant que le terme [=] dénote la partie entière du nombre réel (=). T, est le pas 
e e 


LEA s a T(a+1) 2 A PA 2 
d’échantillonnage et ( n = rOtDr@ jh représente le binôme de Newton généralisé. 


Maintenant, pour le cas particulier de &« = 1, l’expression (3.29) peut être simplifiée comme 


suit : 


TO+1.T() 


KDE LE] = lim, POUR free fŒ-i. 1))} (3.30) 
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Où cette dernière expression représente l’intégral de Grunwald — Leitnikov. Cette 


intégration est définie, pour & = —2,—3,:., par l’expression ci-dessous : 


(CD (;) ft j. )}e20 


3.3. Transformée de Laplace des opérateurs d’ordres fractionnaires: 


É to 


le LE] = SDe [FE] = ut (n° Le 


3.3.1. Transformée de Laplace de l’intégration fractionnaire: 


e Théorème : 
La transformée de Laplace de l’intégration fractionnaire, donnée par les trois définitions citées 


précédemment pour n — 1<aæœ<n , est exprimée comme suit [69,701]: 


LT EVOD=SÈSE6) (3.32) 


Preuve : 


t_f@ 
LPCRIEL CE) = £ a JL, eo dr) F8 


Prenant en considération le théorème de convolution suivant : 


LC fl HO) PTO) 6.34) 

LEON) = LE 1 LE) (3:35) 
__ 1 2. 

+ | (3.36) 

= 574, £[f(t)] (3.37) 

= s7%,F(s) (3.38) 


3.3.2. Transformée de Laplace de la dérivée fractionnaire: 


Par l’utilisation de l’intégrale de Riemann — Liouville, il est possible de définir la dérivée 


fractionnaire de f(t), pour à < 0, par l’expression (3.39) suivante [69,70] : 


LGDEFOD = LG (fe dr) (3.39) 


En se basant sur le théorème de convolution, on obtient: 


L(RDETF(E)]) = Leo (3.40) 


TE a) 
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ee TO) (3.41) 
=$%.F(S) (3.42) 


Similairement, pour n ordre entier, la transformée de Laplace d’une dérivée 


fractionnaire quelconque peut être donnée comme : 
n—-k-1 
LCVONESIO Eee (3.43) 
0 k=0 


atn-k-1 


En général, il existe deux définitions de la transformée de Laplace d’une dérivée 


fractionnaire. Ces dernières sont exprimées comme suit : 
3.3.2.1. Transformée de Laplace de la dérivée fractionnaire selon RL : 


Selon la définition de Riemann — Liouville, la Transformée de Laplace de la a*"Ædérivée 


de la fonction f(t) est donnée par [69,70] : 

LODEL OD = LCPDPDE TE) (3.44) 
A partir de cette expression, on peut écrire : 

LCD) = ss EF (s)— 305 Dr EL EDE RE F(0) (3.45) 


Or, la transformée de Laplace de f(t) peut être définie, en fonction de la valeur de &, par: 


s%.F(s) ‘a <0 
LODEL (ED = 4 F(s) :æ=0 (3.46) 
SAPCI= TC DE PALECO)]) :a>0 


En général, la transformée de Laplace de la dérivée de Riemann — Liouville est 
souvent très limitée en pratique, parce que elle nécessite la connaissance a priori de 
l'information sur la condition initiale f(to) (ie, à l’instant ts = 0), ce qui est en fait 


dans la science de l’ingénieur une quantité physique qui est difficile à mesurer. 
3.3.2.2. Transformée de Laplace de la dérivée fractionnaire selon Caputo : 


Selon la définition de Caputo, la transformée de Laplace de la a*”®dérivée de la 
fonction f (t) est définie par l’expression (3.47) comme suit [69,70] : 
LOGDELF D = s°.F(s) - ZÉcos* **D*f(0) (3.47) 


Sachant que la valeur de & € R* est définieparn<a<n+letnenN 
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Preuve : 
LOGDELF (D = LCODE "D" f(E) (3.48) 
= EP) ne SE DEF (0) (3.49) 
Supposons que k = m—1—1 
LODEL OD = s%F(s) —- Xi st "DEF (0) (3.50) 


Cette dernière transformée peut être donnée en fonction de & comme suite : 


s%.F(s) ‘a <0 
L(GDETF (ED = F() ‘a =0 (3.51) 
S%, F(S) — XE2o(s "#1. 6DÉ[F(0)]) :a>0 


Dans les sciences de l'ingénieur, l’état de repos des systèmes physique correspond à & < 0 ou: 


L(GDELF D = s%.F(s) (3.52) 


Ce qui représente l’avantage majeur de la définition de Caputo sur les deux autres 
définitions précédentes. Par conséquent, la définition de Caputo est la plus adoptée par 


les chercheurs pour l’étude et l’analyse des systèmes fractionnaires. 


3.4. Approximations des transferts fractionnaires: 


Cette partie à pour but de présenter les différentes méthodes d’approximation qui 
assurent le rapprochement des pondérations rationnelles ajustables par des filtres 
équivalents d’ordre entier. Cette étape est donc indispensable lors de la robustification 
du contrôleur H, primaire. Pour atteindre cet objectif, il existe, dans la littérature, 
plusieurs méthodes d’approximation. Certaines méthodes de celles-ci ont été basées sur 
l’approximation de la dérivée ou l’intégration fractionnaire. Cependant les autres 
méthodes ont été basées sur l’identification paramétrique des filtres d’ordre entier dont 


ses structures doivent être pré-choisies par l’utilisateur. On aura donc : 


3.4.1. Méthodes d’approximation de la dérivée ou l’intégration fractionnaire: 


Ces méthodes d’approximation nécessitent d’abord à un changement de variable qui doit être 
effectué sur la pondération fractionnaire à rapprocher. Le principe adopté est s’énoncé comme 
suit : Etant donnée la structure générale de la pondération fractionnaire ajustable comme suit : 


MO (Er ER (3.53) 


No'S+n1 
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Mo'S+M: 


Supposons la variable complexe, notée p € C avec (p = ). D'où, la pondération 


No'S+n1 


donnée par l’équation (3.53) peut être réécrite comme suit: 
W() = p° (3.54) 


Il convient à noter qu’il existe plusieurs méthodes assurant l’approximation du nouveau 
transfert fractionnaire donné par l’équation (3.54) en un filtre d’ordre entier dont la structure 


générale est souvent définie par un gain k,, et N nombre des pôles et des zéros. On obtient donc : 


7 a 
W(p) = k,: RE (3.55) 


A partir de l’équation (3.55), il est facile de déduire la pondération irrationnelle ajustable 


souhaitée en substituant la variable complexe p, dans l’équation (3.55), par celle donnée par s = 
Mi -N1.D 


Mo—No-P 
Les méthodes d’approximation de p® peuvent être résumées comme suit: 


a. Méthode de Carlson : 


La méthode de Carlson [71] est basée sur le processus itératif de Newton — 
Raphson lors de l’approximation de p%. Récemment, plusieurs chercheurs [72-75] ont 
adoptés cette méthode pour rationaliser le terme dep%. Ceci nécessitant les étapes 
suivantes : 

Supposons, tout d’abord, le transfert d’ordre entier W(p) qui approxime, en module 


et en phase, le transfert fractionnaire W(p) tell que : 
= q 
W(p) = [W)] (3.56) 


[W(p)] 4 = W(p) = 0 (3.57) 


q : représente l’ordre fractionnaire de la fonction de transfert. Supposons que : 


q=— (3.58) 
En tenant compte la condition initiale suivante: 
WoG) = 1 (3.59) 


D'où, il est possible d’approximer W(p) d’une façon itérative en utilisant formule 


récursive donnée par (3.60) comme suit : 


@=r)(Wz-1(p)) +(v+r) (pb) 


3.60 
+7) (Wz-1(p)) +7) W(p) Qi 


We) = We-1(p). 
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Sachant que le nombre d’itération, noté kXmax, doit être a priori fixée par 


l’utilisateur. Si l’on suppose maintenant que : 


_ 174 ; 
is TT avec : (3.61) 


= (3.62) 


7 1+a 
W(p) = p{ (3.63) 


Alors, sig > 0, ce qui implique que & < 1. Dans ce cas, le transfert se comporte comme 


dérivée W,(p) à l’itération finale kmax: 
Wo (D) = Wimax D) = WE) =p (3.64) 
Par ailleurs, sig < 0, ce qui implique que a > 1. Dans ce cas, le transfert se comporte 
comme intégrale W,(p) à l’itération finale kmux : 


MCD) = WimxP) © WE)=> (3.65) 


b. Méthode de Charef : 


ES 


Cette méthode à pour but de rapprocher l’intégrateur fractionnaire W,(p) = p * par 
un transfert d’ordre entier W,(p) possédant N zéros et N +1 pôles. Ces deux derniers 
doivent être répartis uniformément sur la plage fréquentielle[w,,w,]. Dans ce cas, le 
terme p * peut être considéré comme étant un Pôle à Puissance Fractionnaire (Qppp), 


donné par : 


QpPpr(s) = 


K 
: 7 , AVEC Op < We < Wy (3.66) 


Sachant que les deux paramètres K°,, et w, sont définis comme suit : 


Ko = (Yo) (3.67) 


&, = 4]—1 + 10Ga) (3.68) 
Avec : 
e K, : Le gain du PPF. 


e &, : La pulsation de coupure du PPF. 
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e € : L’erreur maximale commise entre la pente de la réponse fréquentielle de l'opérateur p * et 


celle du PPF. 


Le transfert d’ordre non entier QPpr(p) peut être donc rapproché par l’utilisation de la formule 


récursive ci-dessous : 


nca(1+ d ) 


Les pôles w,,, et les zéros w,7 peuvent êtres déterminés d’une façon récursive comme suit : 


W,) = G. (3.69) 


Dpk = wo. (p.p)*. Jo , k=0,1,-,N _—. 
or = Yao.(b.p}.p  , k=0,1,-,N —1 
Sachant que Ÿ, @ et N sont donnés par : 
Ÿ — 10e) G 7) 
p = 1062) 
l 100.0ph 
N = partie entier Lao) +1 (3.72) 
log(y.p) 


c. Méthode générale du CFE [44]: 


Le principe de cette méthode est basé sur la représentation de la dérivée 
fractionnaire W,(p) =p* par un transfert d’ordre entier en hautes fréquences, noté 
Hyign(p) et un autre transfert d’ordre entier en basses fréquences, noté H,,%(p). Ces 


deux transferts sont définis comme suit: 
p (24 
Haign@P) = (1 + 2) (3.73) 
@e\ À 
Hiow(P) = (1+%€) G.74 


L’application du CFE sur les deux transferts précédents peut fournir les deux expressions 


présentées par l’équation (3.75) comme suit : 


. = &T. CFE{Hyign(p)} si © » _—. 


Wow D) = O6. CFE{Hiowp)} Si & € w 
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Dans ce cas, le transfert d’ordre entier équivalent W,, (p) est obtenu comme suit : 


(2 K(K+a) [D\ kk-a) [D N 
Po = 08.[0:3, SE, (EEE) BERG) 676 
k=1 
a.(2°) K(k-@) (2°) K(k+a@) (2°) N 
,.,(p) = us. [04 TZ jee D ue ne À (3.77) 
k=1 


d. Méthode d’Oustaloup : 


Cette méthode est largement utilisée dans le calcul fractionnaire, notamment dans le 
domaine  d’automatique. Elle est connue aussi sous l’acronyme  CRONE 
(Commande Robuste d'Ordre Non Entier)|41,761. L’idée de méthode est basée sur 


l’approximation du terme fractionnaire W,(p) = p* 


par un transfert d’ordre entier 
équivalent possédant N pôles et N Zéros attribués uniformément sur la plage 
fréquentielle [w,,w,]. Le transfert d’ordre entier résultant W,(p) est s’exprimé 


comme suit: 


Wits)= Ci” (Etes) (3.78) 


P+Op,k 


Où le gain Ci, les N pôles w,,4 etles N zéros w,, peuvent être déterminés comme suit : 


Ca = Op? 
2k-1+a 
Opk = Op. (@y) 2N (3.79) 
2k-1-a 
DK — Op. (CM) 2N 
Avec : 
&, = (2) (3.80) 


En général, le gain C; doit être ajusté de telle sorte que la condition [(jw)*{,,-1 = 
0 dB, Vae R est toujours remplie. Pour assurer un placement alterné équidistant entre 
les N paires w,7 et w,, dans une échelle logarithmique de fréquences, on doit vérifier 


le rapport suivant : 


H 


Ozkt1 _ Opkti ()" (3.81) 


z,k Dp,k Oh 
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3.4.2. Méthodes utilisant la réponse fréquentielle du transfert fractionnaire: 

Cette classe des méthodes d’approximation est basée sur l’évaluation de la réponse 
fréquentielle du transfert fractionnaire dans une plage de fréquence bien précise. Cette 
réponse est ensuite comparée avec celle prédite par le modèle irrationnel on fait 
ressortir la norme H, de l’erreur d’approximation. La minimisation de ce dernier par un 
algorithme adéquat fournie les paramètres optimaux du modèle d’ordre entier souhaité. 


Parmi les méthodes existant dans la littérature, on peut citer : 


a. Méthode de Hartley : 


L’algorithme de cette méthode a été développé par Hartley et Lorenzo en 2003 [77] dont 
ses détailles sont disponibles dans la bibliothèque ninteger. Il est basé sur l’évaluation de la 
réponse fréquentielle du transfert fractionnaire dans une plage de fréquence présélectionnée par 
l’utilisateur. Cet algorithme fournit soit un modèle rationnel ou un modèle irrationnel présenté 
par la structure générale ci-dessous : 


1 
CnSQ+cn-15 M DQ+..+c 520 +050 +00 


W(s) = (3.82) 
De l’équation (3.82), si la valeur de Q est choisie entière, on obtient, par conséquent, une 
approximation par un modèle d’ordre entier. Sinon, l'algorithme fournit un modèle fractionnaire 


de type commensurable. 
b. Méthode de Levy 


La fonction de Levy à pour but d’identifier les paramètres d’un transfert rationnel 
ou 1irrationnel présenté, qui peut être donné par la structure générale (3.83). Cette 
identification  paramétrique est assurée à partir de l’extraction de la réponse 
fréquentielle du transfert fractionnaire W(s). La réponse fréquentielle du modèle 
prédite d’ordre entier W(s)est ensuite évaluée dans une plage de fréquence. La 
comparaison de ces deux réponses fournie une erreur d’approximation et la 
minimisation de la norme H, de ce dernier par un algorithme d’optimisation fournie les 
paramètres optimaux de ce modèle. La structure générale de ce dernier est donnée 
comme suit [78]: 


W(s) = bn" +bn 15-00 4p,, 50-20 4.4 ps +b, (3.83) 
ans"@+an_15 8-00 +an 25-20 +...+a1s0+1 : 
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La norme H, de l’erreur de modélisation est donc s’exprimée par : 
E = W(jw)[a, Go)" ++ a (jo) +1] - [b,Go)"e + + b Go) + bo] (3.84) 


L'inconvénient majeur de cette méthode est qu'elle n'est pas appliquée que dans une plage de 


fréquences très restreinte. 


c. Méthode de Venagre 


Le principe de cette méthode est basé la minimisation de la norme H, de l’erreur 


d’approximation qui est préalablement pondérée par le vecteur v . On aura donc [78] : 
E' = y. {W(jow)[a, Go)" ++ a, (Go) +1] - [bn Go)" ++ b (jo) + bol} (3.85) 


Les composants constituant le vecteur v dépendent par les pulsations w;, sachant que i = 


1,2; …; f. Ces composants sont exprimés par: 


"21 si i=1 
2wf 
Di+1—Wi-1 : : 
AS NE 1<i<f (3.86) 
LS A D : _ 
ST 1 — 
2wf F 


Pratiquement, cette méthode a été proposée afin d’améliorer la qualité d'approximation du 


transfert fractionnaire notamment dans les basses fréquences. 
d. Identification fréquentielle par la commande fitfrd du Matlab : 


En général, la fonction fitfrd est disponible dans la bibliothèque Robust/Control/ 
Toolbox du logiciel Matlab. Elle est utilisée pour approximer un bloc de transfert d’ordre non 
entier dont sa réponse fréquentielle doit être rapproché, en module et en phase, à celle prédite par 
le modèle d’ordre entier. La syntaxe sous Matlab de cette fonction est très simple, premièrement 
il faut évaluer le transfert fractionnaire W(jw) dans une plage fréquentielle [w, wn] réparties de 
façon logarithmique par la commande frd. Ensuite, on exploite les données d’entrée/sortie 
fournies par la première commande, en utilisant la commande fitfrd et en ressortissant une 
représentation d’état d’ordre entier qui modélise le transfert d’ordre entier W(jw). L’ordre de ce 
dernier est ajusté par l’utilisateur afin d’aboutir une erreur de modélisation &,, qui doit être assez 


réduite et satisfaisant la condition suivante : 


D Go) — WGw)|, < &n (3.87) 
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3.5. Formulation du problème d’optimisation sous contraintes : 
Nous rappelons 1ci que, le problème de la robustification du contrôleur H, primaire basée sur les 
pondérations fractionnaires ajustables, peut être reformulé comme suit: 


Omax [ws,C Xe): S,, (jæ))| + Omax [Wr,C Xe &). 7, Uw)| 
minx, [Max + (3.88) 


Omax [Wu,C Xe ©). Ki * S,y(jw)| + Omax [Wse, € oh Sy, (j&) *GN Go)| 
En tenant compte-les contraintes de saturation ci-dessous : 


PT UR (3.89) 


XCmin 
Avec We f° Wy P° Wy Fr © Wsc F Sont les pondérations fractionnaires ajustables limitant les 


sensibilités S,,, T,,, (Ki 5,) et CS . Gn) respectivement. La résolution de ce problème doit 
permettre non seulement d’assurer la sélection automatique de ces pondérations, mais également 
elle doit assurer l’amélioration de la marge initiale de robustesse RP fournie par le contrôleur 


H> primaire. 
3.5.1. Solution du problème de robustification par l’algorithme d’essaim particulier: 


L’optimisation par essaim particulaire est une méthode développée en 1995 aux Etats-Unis 
sous le nom anglais "Particle Swarm Optimization PSO ". Elle a été proposée par les deux 
psychologues Russel Eberhart et James Kennedy [80, 81]. Historiquement, l'idée de base 
cette méthode est inspirée de la nature. Elle consiste à chercher de modéliser des interactions 
sociales entre des « agents » devant atteindre un objectif donné dans un espace de recherche 
commun. Chaque agent ayant une certaine capacité de mémorisation et de traitement de 
l'information. La règle de base est qu’elle ne devrait y avoir aucun chef d’orchestre, ni même 
aucune connaissance par les agents de l’ensemble des informations, seulement des connaissances 
locales. Un modèle simple fait alors élaborer. Dès les premières simulations, le comportement 
collectif de ces agents évoqué celui d'un essaim d’être vivant convergeant parfois en plusieurs 
sous essaims vers des sites intéressants. Ce comportement se retrouve dans bien d'autres 
modèles, explicitement inspirés des systèmes naturels. Ici, la métaphore la plus pertinente est 
probablement celle de l’essaim d’abeilles, particulièrement du fait qu’une abeille ayant trouvée 
un site prometteur sait en informer certaines de ses consœurs et que celles-ci vont tenir compte 


de cette information pour leur prochain déplacement. Finalement, le modèle s’est révélé être trop 
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simple pour vraiment simuler un comportement social, mais par contre très efficace en tant 
qu’outil d'optimisation. 

La version historique peut facilement être décrite en se plaçant du point de vue d’une 
particule. Au départ de l’algorithme, un essaim est réparti au hasard dans l’espace de recherche 


de dimension N, chaque particule j ayant également une vitesse aléatoire de déplacement. On 
T . . CO . 
note 6; = Ge 7 a) le vecteur qui caractérise la position d’une particule et V; = 


T ss 2 £ 2 , : 
(via Vi, vin) le vecteur de changement de position (appelé vélocité) d’une particule. 


Ensuite, à chaque itération on a : 


e Chaque particule j est capable d'évaluer la qualité de sa position et de garder en mémoire sa 


T 
meilleure performance que l’on note M = cm) , c’est-à-dire la 


meilleure position qu’elle a atteinte jusqu'ici (qui peut en fait être parfois la position 


courante) et sa qualité (la valeur en cette position de la fonction à optimiser). 


e Chaque particule est capable d'interroger un certain nombre de ses congénères (ses 
informatrices, dont elle même) et d'obtenir de chacune d'entre-elles sa propre meilleure 
performance (et la qualité afférente). La meilleure position atteinte par toutes les particules 

D 


de l’essaim est indiquée par un vecteur xb65t.,, = (1h69... x best es qe 


e À chaque itération, chaque particule choisie la meilleure des meilleures performances dont 
elle à connaissance, modifie sa vitesse en fonction de cette information et de ses propres 


données et se déplace en conséquence. 


Une fois la meilleure informatrice détectée, la modification de la vitesse est une simple 
combinaison linéaire de trois tendances, à l’aide de coefficients de confiance qui sont les 
suivants: coefficient d'inertie de la particule (noté ici parc, ), coefficient d'attraction de 
la particule vers sa meilleure position (noté ici parc.) et finalement, coefficient d'attraction 
de la particule vers la meilleure de ses voisines (noté 1c1 parc,). Encore une fois on détermine 


cette vitesse, la nouvelle position de la particule est obtenue en ajoutant cette valeur de vitesse à 
la valeur de sa position obtenue à l'itération précédente. Ce processus se répète jusqu'à 


l'obtention d'une précision fixée à l'avance par l'utilisateur. 


Les équations de mouvement d’une particule, pour chaque itération "£", sont définies par [39]: 


VÉtT = co Vf + curé: Gi — xé) +o' ri: .…. — x£) (3.90) 
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A = XV (3.91) 


Le mauvais choix des paramètres d’initialisation peut dégrader les performances de 
l’algorithme PSO mais n’empêche pas d’obtenir des solutions néanmoins acceptables de point de 
vue pratique. C’est pour cette raison que cet algorithme est considéré comme un algorithme 


d'optimisation robuste. 


Le paragraphe suivant récapitule les étapes à suivre de l’algorithme de la robustification du 
contrôleur HX primaire dans laquelle on introduit l’algorithme PSO pour assurer la sélection 
automatique des pondérations fractionnaires adéquates, puis on applique l’approche 
d'identification fréquentielle pour approximer ces pondérations en des filtres d’ordre entiers. Ces 
derniers sont utilisés pour formuler le problème de sensibilité pondérée ce qui est résolu par la 
méthode H,, utilisant l’approche LMIs. Ce processus se réitère jusqu’à l’obtention d’une marge 
de RP, meilleure que celle fournie par le contrôleur H,, primaire. Ces étapes sont résumées 


comme suit : 
3.5.2. Algorithme de la robustification basé sur les pondérations fractionnaires ajustable : 


%  Initialisation de l’algorithme de robustification : 


Dans cette étape, les transferts du modèle nominal G,(s) et les pondérations irrationnelles 
fixes Ws,(s), Wr, (s) , Wu, (s) et Wsc,(s) sont supposés connus. Sachant que, le transfert du 
contrôleur H, primaire K,(s) est assuré par la méthode H,, utilisant l’un des deux algorithmes 
cités en premier chapitre. Ce dernier permet de fournir la marge de robustesse initiale RP. 
Finalement, les limites inférieure et supérieure de l’espace de recherche sont données à partir des 
règles de réglages précédente ainsi que le nombre des paramètres à optimiser est supposé égal 
à N. 

* Initialisation de l’algorithme PSO 
A l'itération ? = 0, la position de n, € N particules Ge) est initialisée de façon aléatoire 
Ï=1,..Np 


dans l'intervalle | L on obtient : 


X .…. XX 
Cmin, Cmax, 


£ — ? 4 € £ £ L L £ LA L A £ £ £ £ £ 
Xc; nn M$, Es; &B,, M) » MT, ET @Br » M5; Eu OR Vj » MSG ES y WSG » l' (3.92) 
© 


£ ê £ £ 
W W W, W. 
sf} Ti Us sf 
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e Etape C.1 : pour chaque particule Xp faire les sept étapes suivantes: 


> Etape C.I.I : Calculer les pondérations fractionnaires ajustables suivantes : 


(Ws, o ) Wr,, (s, ) Wüs, (s, x), Wscy, (s, xé,) 


> Etape C.1.2 : Calculer les réponses fréquentielles suivantes : 


(We, («, 2} Wr,, (&, 2) Wos, (o, “ee Wscy, («, 2e) 


> Etape C.1.3: Appliquer l'algorithme d'’'approximation utilisant l'approche d'identification 


fréquentielle. On obtient donc les filtres d'ordre entier ci-dessous : 
LRO ROL MOLAO)E 


> Etape C.1.4: Déterminer le système augmenté à partir de ces filtres d'ordre entier et le modèle 
nominal G\(s) ; 


> Etape C.1.5 : Formuler le problème de sensibilité pondérée, ensuite le résoudre par la méthode H, . 


On obtient le contrôleur H,, suivant: {ké, (s)}; 
> Etape C.1.6: Déterminer les nouvelles sensibilités 15 (s), 4.6), Ki, (s) "5e (6) Gn(S) 
j j j 
€ ; 
So (5): 


> Etape C.1.7 : Evaluer la fonction objective RP Ge ) avec : 


F = € ct ere L pl .ctl nr... ct 
RPÉ, = IL Sÿo,|L. + IL2 Th]. + IL Ko | + IL2 Gn "So | (3.93) 
e Etape C2: sachant que l’étape précédente est achevée par l’obtention n, marges de 


best, A d best,l 
swarm 


robustesses telles que : (RP, RPS, ++ RPÉ . D'où, les solutions x 


sont obtenues comme suit: 


xpestl = min {RPS ,0<i<£} 
J xt. J 


best,t : € ï (3.94) 
Xswarm — ru {rpg vi} 
6j 


Ensuite, La position et la vitesse de chaque particule sont mises à jour en fonction des deux 


expressions (3.90) et (3.91) respectivement. Aller à l’étape C.3. 


e Etape C.3 : test d’arrêt : 
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Le critère d'arrêt est associé par la plus faible marge de robustesse, notée RP;, sachant qu'elle est 


donnée par: RPnin = J CG) Si la condition RPmin < RPo n'est vérifiée, le transfert du contrôle Hs 


robustifié est celui donné par K{s). Sinon, on doit faire les changements ci-dessous : 
Ÿ__ Incrémenter le nombre d'itération € par 1, ie. £ <£+1; 


Ÿ”_ Mettre à jour la valeur de RP, par celle de RPin , te. RP9 — RPin : 


Ÿ” Mettre à jour le transfert du contrôleur KQ(s) par celui fournissant la marge de robustesse 


RPmin 


Retourner à l’Etape C.1. 


Remarque : 


Dans l’étape C.1.4, si la solution de la méthode H, est vide (ï.e., {K4,(s)} 


CL 1). 


hefiny) — 
Dans ce cas, la particule * fournissant cette anomalie doit être réinitialisée aléatoirement dans 
ses limites correspondantes x, < va < Xe, avant de retourner à l’étape C.1.1. 

Il convient à noter ici que, le problème de sensibilité pondérée ne peut pas être résolu par la 
fonction Hinflmi du logiciel Matlab avec des pondérations fractionnaires ajustables. Il faut 
tout d’abord les approximer par des filtres irrationnels dont ses transferts doivent être fournis 


avec des ordres très réduits. Cet objectif est assuré par l’identification fréquentielle basée sur la 


commande fitfrd du logiciel Matlab. Les étapes de cet algorithme sont [40] : 
3.5.3. Algorithme d’identification fréquentielle: 
Dans cette partie, on suppose qu’on a les transferts des pondérations fractionnaires ajustables 
suivantes W£ (s, x wi (s, a WÉ (s, x) et W{ (s, } Sachant que l’approximation 
fj Cj fj Cj fj j fj j 
souhaitée de ces pondérations est effectuée dans une plage de fréquence définies a priori par w € 


[min » Wmaxl. Les filtres irrationnels correspondants d’ordre réduit peuvent être obtenus 


comme suit : 


# Etape 1: Déterminer les réponses fréquentielles des pondérations suivantes W£ s, xt.) 
Ip: P q P Sf, c 


€ £ € € £ ( £ ) 
We, (s, 2) ; Wu, (s, xe) et Wscs, S,Xc;) 


+ Etape 2: Mettre un choix plus faible des ordres entiers : ls.» NT,» Nu, el Nsc,. Pour les filtres 
j j j j 


TE 
souhaités Wé (s),Wf, (s), Wé, (s), Wie, (s) approximant les  pondérations  fractionnaires 
j j j j 
(4 € A (A (A LA € LA ; ie 
Ws,. S,Xc; }» Wr,. S,Xc; }» Wu. s,Xc;) et Wic r, \S*c;) respectivement. De plus, choisir les 
j j j j 


précisions d’approximation souhaitée : AWs , AWr, AWy et AW, puis aller à l’étape suivante. 
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+ Etape 3 : Appliquer séparément la commande fitfrd sur chaque réponse fréquentielle WS.. («, ae) 
j 
WË («, A W£ " (o, x.) et W£ ” («, 8). On obtient, par conséquent, les filtres d’ordre 
j j j 
entier: Wé Es) Wé, CS); Wé, {s), Wie, _(s). Aller à l’étape suivante. 
j j j j 


+ Etape 4 : Vérifier, pour chaque approximation, les critères d'arrêts suivants: 


| Ws,, (w,xe,) _ Ws,, || < AWS; IL22 (&, ke) ” Wr, (@, | < AWr 


LA : (&, 2) — Wé, (ox) Î, < AWy . LE (o,x£) — We, (&,%) Î, < AWsc 


e Si:les conditions sur AWe, AWry, AW, et AWG sont satisfaites, arrêter l'algorithme. 
e Sinon: mettre, la ou la précision souhaitée n’est pas atteinte, les incréments ci-dessous: 
— 
Reese +1 
— 
Nr ; + 1 
= 
Ny f + 1 


nsc; an nsc; +1 


Puis retourner à l’étape 3 


3.6. Simulations et interprétations des résultats : 


Dans cette partie, la robustification proposée dans ce chapitre est validée sur le même système 
hydraulique cité précédemment. Nous rappelons que l’objectif de cette robustification réside 
dans la satisfaction de la condition RP et de l’assurer avec une marge de sécurité acceptable du 
point de vue pratique. Sachant que cet objectif est assuré par l’optimisation adéquate des 
paramètres des pondérations dont l’allure souhaitée de chacune est imposée par un gabarit 
correspondant prédéfini par les mêmes règles de réglage. La sélection automatique lors de cette 
robustification est effectuée dans un espace de recherche plus large que celui mentionné en 
deuxième chapitre. Autrement dit, l’algorithme de PSO à la capacité de fournir non seulement 
des pondérations irrationnelles ajustables. Il peut également fournir des transferts fractionnaires 
définissant les pondérations souhaitées. Ces transferts doivent être préalablement rapprochés par 
des filtres d’ordre entier en utilisant l’approche basée sur l’identification fréquentielle. 

En comparaison avec l’ancien espace de recherche cité précédemment, le nouvel espace est 


donc imposé comme suit 
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| Ms, <2 
S° 11.3 < m, < 15 


4< My, < 5.66 
Ui Ê < OBC; < 3.33 (3.95) 


1.75 < Mr, < 2.07 
Wr.: 0.5 < OBT, < 0.868 
15 <n <2 


Sachant que tous les paramètres des pondérations fractionnaires désirées Wé,. (s, se) et 
J 


Wôs, (s, 4) et Wf f; (s,xt,) doivent être optimisés pour j = 1,2,--,n, particules pendant 


toutes les itérations £ = 0,1,:::,#m4x. À cet effet, les paramètres de réglage de l’algorithme 


PSO ont été choisis comme suit : 


=" Nombre des paramètres à optimiser : N = 7 

= Nombre maximum des itérations : lnax = 160 
" Taille de l'essaim : n, = 16 

CoOmax — C0 


"Facteur d'inertie Co = Coux — ( mn) 4 OÙ Comux — 09, Coin = 0-4 


Ümax 
“ Taux d'apprentissage cognitif: C1 = 1.65 
“ Taux d'apprentissage social :c3 = 1.75 


Dans la suite de ce travail, nous allons détailler l’algorithme de robustification proposé en 


utilisant les données citées ci-dessus. On obtient donc : 


A l’itération initiale (1.e., { = 0), on obtient les 16 positions initiales des particules ci-dessous: 


M, M, M$ 
ns, ns, ns, 
M5, M5, Môe 
x = | we, |, = | wie, |... 2e = | @c, (3.96) 
MY, M7, MF, 
CB, WBr, pre 
mr, mr, Me 


Où chaque particule est aléatoirement choisie entre ses bornes inférieure et supérieure comme 


suit : 


M M, n8, ac, 
: £<2,4< : £5.66,13<| : |<15,2< : < 3.33, 


1.5 < : 
0 0 0 
Ms. My, Nic 


0 
OBCi6 
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MY, Br, nr, 
1.75 < : < 2.07,0.5 < : < 0.868 et 1.5 < : < 2. 
Me OBTie Ne 


Ÿ” L'étape C.I.I fournie les 32 pondérations fractionnelles ajustables suivantes: 


0 0 
s 1 es s 1 Tic 
w*9 0 w*9 n°. 
BT: T1f 0 BT16 T16 [0 
t, (522) _ be (5,x8) L 
Wr. (s,xc ) = | ——— W. s, À = 
Tr, \SXc 0 pes WT Xci6 0 
fi 1 le, f16 rs 
T1 T16 
mn EL 5 s+1 
OBT: OBTi16 
et 16 pondérations irrationnelles ajustables suivante : 
S 1 S 1 
M Um M M 
0 0 \ — 1 1 0 0 _ 16 16 
Wo, (s, xe,) ou Ù ue Wie (s, 7) e) 
's+i de s+l 
BC] WBC16 


Dans l’étape C.1.2, la réponse fréquentielle de chaque pondération fractionnaire est évaluée 
pour 1000 fréquences réparties uniformément dans la plage de fréquence w € [107#, 20] 


radians/secondes. Cette étape fournie les réponses fréquentielles ci-dessous : 


| n° . n° 
jJ'or de +0 S1 J'or L w*° S16 
nÿ Bi nÿ Bic 
1 2 16 M2 
We (ox )= .  (e)= . 
S TAPE on 0 gran} S Ta NC 7. 0 
f 1 Lo f16 16 S16 
*0 0 *0 0 
J (0/2 SL Op Es, J (0/2 LD OB,, S16 
0 
J'or 1 Mrs J'or 1 Ti16 
*0 0 +0 0 


6 . 
+6 —Ù ——J'ort+1 
BT1 SBT16 


Avecr = 1,2,..,1000, &w, = 107% et 900 = 20. De plus, l’algorithme d’identification 
fréquentielle, cité dans le paragraphe 3.4.3, est ensuite appliqué avec le choix des deux ordres 


entiers suivants: ns L 2etnr = 2 et les deux précisions souhaitées suivantes: AW, = 107%et 
J p 
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AW% = 107$. Finalement, cet algorithme est achevé par l’obtention de 32 filtres d’ordre entier 


équivalents approximant les 32 pondérations fractionnaires citées ci-dessus. Ces filtres sont 
ensuite utilisés avec les 16 pondérations irrationnelles Wy, (s, x0.) pour augmenter le modèle 


nominal Gy(s). On obtient donc les 16 problèmes de sensibilités pondérées suivants : 


W5,, (s).Sy, () WS,, (S).Sy, (s) WS,,().55,, CS) 
W, (s). K6, (s). Sy, Wa, (s). Ko, (s). Sy, ee, WC). KQ, (s). Syn 
UMROBAON | RUROSAOE ER Wri S) Try, CS) 


Dans l’étape C.1.5, on obtient 16 contrôleurs H,, , qui sont: Kç.(s), Kç,(s), ..…, Kç,,(s) qui 
avec lesquels on peut déterminer les sensibilités (So, : K, "Sy, Ty); (So, | KG, "Sr, Ty.) 


0 : 0 ,ç0 
6 SR 


Vos mn) Ces dernières permettant de fournir les 16 marges de 


robustesses RP, qui sont obtenues comme suit : 


0 0 , <0 0 0 0 .p0 . «| 
RP9 = ILE ve + [w?, me + lus, KP, -S, 
0 0 0 0 0 0 0 , € 
rpg, = |W£.-s + [we er + [we ee 
16 116  Yo46 . 116 Yo6 | 16 1 ul, 
A . 2 : best,0 best,0 
Ces marges sont comparées entre elles afin de déterminer les deux vecteurs x; SX Sparm 


: : : : best,£ : . N 
sachant que ce dernier fourni la meilleure marge, .e., RPyin = J a Finalement, la mise à 


jour de la position et celle de la vitesse de chaque particule sont assurées par les deux équations 
(3.90) et (3.91) et toutes les étapes précédentes sont réitérées jusqu’à £max = 160 itérations. La 
Fig.3.2 représente l’évolution de la fonction objective pendant les 160 itérations pour 15 


exécutions de l’algorithme de robustification proposé. On obtient donc : 


(Ce) 20 40 60 80 100 120 140 160 
Générations 


Fig.3.2. Evolution de la fonction objective pendant les itérations de l'algorithme PSO 
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La Fig.3.2 montre clairement que parmi les 15 exécutions de l’algorithme de robustification 
proposé, il existe certains cas où la condition RP est satisfaite à partir de l’itération { = 30. Ceci 
confirme aisément que l’utilisation des pondérations fractionnaires est indispensable pour 


améliorer cette marge de robustesse. La meilleure amélioration de RP est assurée par le vecteur 


best,£ : : & : : Joe ; 
Xswarm — » Ce qui fourni les paramètres optimaux des trois pondérations ci-dessous : 


{ Ms, = 1.82 < 2 si 


à ; 15 —+0.02 
_ = W. re ) = —"#2 
1.3 < m1 = 1.37 < 15 Sf Coptimal s+0.02-V1073 


S 1 
Pmax as _ zos "287 
Wu (SAcoptimar) io (3.97) 


My, = 482 < 5.66. 
OBC; = 2.05 < 3.33 


Mr, — 1.92 < 2.07 — 

0.56 1. 
©pr, = 0.56 < 0.868 = We (s, ee.) - E =) 
1.5 <n = 1.75 <2 


Sachant que les deux filtres irrationnels W;, (s) et Wr., (s) approximant les deux pondérations 


| 4 £ £ £ ee 
fractionnaires M," (s,x ss ) et MT (s,x te ) sont obtenus par les deux précisions 


Coptimal Coptimal 


2.6952 x 107“ et 2.6418 x 107% respectivement. Les transferts correspondant sont donc donnés 


par : 
___ 0.454(5+0.033)(5+0.019)(s+0.006)(5+0.0015) 
Ws, (s) = (s+0.013)(s+0.0036)(5+0.00081)(5+0.0000037) (3.98) 
Wy, or 1012.3(5+6.817)(s+1.742)(s+0.5057)(5+0.3484) (3.99) 


(s+35.55)(s+18.92)(s+4.86)(5+1.2444) 

La Fig.3.3 représente le rapprochement, en module et en phase, des deux réponses fréquentielles 

fournies par la pondération fractionnelle ajustable Wimax (s, x imax ) et le filtre d’ordre entier 
f optimal 


équivalent W. (s). De plus, La Fig.3.4 représente le rapprochement, en module et en phase, des 


deux réponses  fréquentielles fournies par la pondération fractionnelle ajustable 


ê 4 ; _——— 7. 
ME (s, nu) et le filtre d’ordre entier équivalent W7. (s). On obtient : 
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Approximation de la pondération rationnelle Wsf1 


Magnitude (dB) 


—— W1_f_sys 


— — Wi_response 


-45 | 


Phase (deg) 


-90 


———— WA1_f_sys 


—_— Wi_response 


-135 L 
107$ 


102 


10°! 


10° 


Frequency (rad/s) 


101 


Fig.3.3. approximation de la pondération fractionnaire Ws, parle filtre d'ordre entier Ws, 


Approximation de la pondération fractionnaire WTf 


————— W3_f_sys 


— — W3_ response 


Magnitude (dB) 
Ë 
T 


-20 
135F TT 
S o0+- — W3_f_sys 
S 
La — — W3_response 
2 45 
D 


102 10°! 


10° 


Frequency (rad/s) 


10! 


Fig.3.4. approximation de la pondération fractionnaire Wr, parle filtre d'ordre entier Wr, 


Finalement, la méthode H,, basée sur l’approche LMIS fournie un contrôleur H,, robustifié 


d’ordre 14 et un niveau de performance égal à y—0.924. D'où, le tableau.3.1 ci-dessous 


récapitule le gain, les pôles et les zéros de ce contrôleur. On obtient : 


Tableau.3.1 récapitule le gain, les pôles et les zéros de ce contrôleur A, robustifié 


Gain 0.038 

Pôles | -0.0008123 | -2.83e-05 | -1.3305+0.088146j | -1.3305-0.088146.j | -0.791+0.75321,j 

zéros | -0.00529 -0.001467 | -0.30745+0.83491,j | -0.30745-0.83491,j | -1.7565- 1.0202,j 
Pôles | -0.791- 0.753217 | -34.94 | -17.62 | -5.233 | -1.252 | -0.503 | -0.02058 | -0.01322 | -0.003681 
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zéros | -1.7565+ 1.0202/ | -38.44 | -31.16 | -18.58 | -4.86 | -9819 | -0.181 | -0.0165 | -0.0145 


Les conditions de robustesses sur RS fournies par la robustification utilisant les pondérations 
fractionnaires ajustable et celle utilisant les pondérations irrationnelles ajustable sont illustrées 
par la Fig.3.5, sachant que T,,.(s) et T,.(s) représentent respectivement les sensibilités 


complémentaire fournies par la première et la deuxième robustification proposées. On obtient : 


Condition sur la Stabilité Robuste 


10 T TT] T LE T LR T A D 1] 
0 
a 
Z 
() -10 + 
© 
É 
x -20 + 
Ê T_y_01 
® 30} Ty 
£ 1W_T 1 
E 
E 4) 
(2) 
2 
3 -50F 
2 
(Q 
> 
-60 F 
-70 L FRERES LA PS EN 1 1 CARRIERE CS EN 1 TASSE A A 1 AR PEN PRE Pl il 
107 10 10° 100 10! 


Pulsation (rad/s) 


Fig.3.5. Comparaison des deux conditions RS fournies par la 1°"° et la 2°"€ robustification proposée. 


De plus, les conditions de robustesses sur NP fournies par la robustification utilisant les 
pondérations fractionnaires ajustable et celle utilisant les pondérations irrationnelles ajustable 
sont illustrées par la Fig.3.6, sachant que S,,,(s) et S,, (s) représentent respectivement les 


sensibilités directes fournies par la 1°"® et la 2°" robustification proposée. On obtient : 
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Condition sur les Performances Nominales 


Valeurs Singulières Maximales (dB) 


-60 1 1 ii ii ii i 1 FO D A PO L ñ puni vil ñ 1 fe ppp il 
10% 102 10° 10° 10! 
Pulsations (rad/s) 


Fig.3.6. Comparaison des deux robustesses NP fournies par la 1°"° et la 2°"® robustification proposée. 
Maintenant, les conditions sur la limitation de l’énergie de commande sont comparées par 


Fig.3.7, sachant que K5: * S,,,(s) et K::S,,(s) représentent respectivement les sensibilités 


fournies par la 1°7€ et la 2°"€ robustification proposée. On obtient : 


Condition sur le signal de commande 
15 2 D ——— "577 D 


10 


-5F ——— K_01xS_y_01 
—— K_1xS_y_1 


— 1\W_U_1 


Valeurs Singulières Maximales (dB) 


25 L i ju ie hi pis i ju à ji il i ii ir iii il i ji ii pp 
1075 102 10° 10° 10! 
Pulsations (rad/s) 


Fig.3.7. Comparaison des deux limitations de commande fournies par la 1°"® et la 2"€ robustification proposée. 
Le compromis de robustesse trouvé dans le chapitre précédent est comparé à celui fournie par 


cette robustification. Cette comparaison est illustrée dans la F1g.3.8. On obtient : 
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Compromis de robustesse 


T A EE EE ER ON | T T RIT] T A A ER ES M | T LL NE | T 


Valeurs Singulières Maximales (dB) 


Compromis assuré par la 1ère ronbutification proposée 
Compromis assuré par la 2ème ronbutification proposée 


24 ji ji ji iiiil ji ii ji iiil ji ii ji iiiil ji ji ji iiiil 
10% 102 10° 100 10! 
Pulsations (rad/s) 


Fig.3.8. Comparaison des deux compromis de robustesses fournis par la 1°"° et la 2°"€ robustification proposée. 


Les deux robustesses de RP fournies par la 1°7® et la 2°"€ robustification proposée sont 


finalement dans toutes les fréquences. Cette comparaison illustrée dans la Fig.3.9. On obtient : 


Condition sur RP 

0 à] 
a ?| 
«A 
nu A4F 
ao) 
œ 
£ 6r 
x 
LL] 
= 8H 
8 
& 10 ——— RP_01 
2 

— RP_1 

Ba - 
o 
® 14 
3 
2 
© - = | 
S 16 

18 - 
-20 PE PA CE ON 2 gi + rit ii à iii di dd pr) 

10% 102 10° 102 10! 

Pulsations (rad/s) 


Fig.3.9. Comparaison des deux conditions RP fournies par la 1°" et la 2°"® robustification proposée 


À partir de toutes ces figures précédentes, on peut confirmer l’efficacité de cette nouvelle 
robustification proposée. La clé de succès de cette efficacité est s’articulé notamment sur 
l’exploitation d’un espace de recherche très vaste que celui utilisé dans les méthodes de synthèse 
existantes. La réussite de cette robustification proposée se traduite par la satisfaction de la 
condition RP dans toute la plage de fréquence (i.e.,RP < 0,p, VWE[@ymin » @maxl), Chose qui 


n’est par atteinte par l’autre stratégie de robustification. 
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Maintenant, le contrôleur robustifié par cet algorithme proposé est utilisé pour commander les 
192 régimes perturbés cités dans les deux chapitres précédents. Pour cela, le système bouclé est 
excité par les mêmes trois entrées exogènes citées précédemment. La Fig.2.10 représentes les 
sorties globales du système bouclé, tendis que la Fig.2.11 représente l’agrandissement de la 


partie transitoire de ces réponses. On obtient donc : 


Fig.3.10. Sorties fournies par le contrôleur H. robustifié pour 192 régimes perturbés 


D’après la Fig.3.10, on peut remarquer que les sorties fournies par ce contrôleur sont 
complètement meilleures que celles produites par les deux autres contrôleurs H,, cités dans les 
deux chapitres précédents. Le tableau ci-dessous récapitule les taux d’améliorations du système 
bouclé par ce contrôleur par rapport à ceux bouclé par les contrôleurs H,, primaire et celui 


robustifié par la robustification utilisant les pondérations irrationnelles ajustables. On obtient : 


Tableau.3.2 récapitule les taux d’améliorations assurés par le contrôleur H,, robustifié basé sur les pondérations 


fractionnaires ajustables 


contrôleur contrôleur contrôleur Taux d'amélioration % 
H> H, H> Ko2(s) par Ko2(s) par 
primaire robustifiè 1 | robustifiè 2 rapport à rapport à 
Ko(s) Ko1(s) Koz(s) Ko(s) Ko1(S) 
Temps nécessaire pour annuler 
22.5 12.00 05.25 76.66 % 56.25% 
l'erreur de poursuite à 5% (£5%) 
Temps de montée t,, 03.92 04.49 03.45 11.98% 23.16% 
Temps de réponse à 5% (£,...,) 22.50 12.00 05.25 76.66 % 56.25% 
Dépassement maximal Dax 12.00 05.00 03.00 75.00% 40.00% 
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La Fig.3.11 représente un agrandissement de la dynamique de poursuite dans l’intervalle du 
temps normalisé t € [o, =] secondes. On obtient 
e 


Fig. 3.11. Agrandissement de la partie transitoire de la dynamique de poursuite fournie par le 
contrôleur H,, robustifié pour 192 régimes perturbés 


La Fig.3.12 représente les commandes nécessaires pour fournir les sorties globales citées 


précédemment. On obtient 


Fig. 3.12. Commandes fournies par le contrôleur H. robustifié pour 192 régimes perturbés. 


D’après la Fig.3.12, on peut voir que toutes les commandes fournies par le contrôleur 4, 


robustifié sont moins couteuses que celles fournies le contrôleur H,, primaire, i.e., Umax = 4.50 


en régime transitoire. 
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3.7. Conclusion 


Dans ce chapitre, nous avons proposé une deuxième stratégie de robustification du contrôleur 
H> primaire et cela dans le but de satisfaire la condition RP en présence les incertitudes 
multiplicative de modélisation. Cette proposition nécessite a priori la résolution du problème de 
la sélection automatique des pondérations fractionnaires ajustables existant dans un espace de 
recherche présélectionné a priori par l’utilisateur. Cette dernière a été assurée par des gabarits 
fréquentiels imposés par certaines règles de réglage. Ces règles peuvent fournir les allures 
souhaitées de différentes sensibilités du système bouclé. Une fois l’algorithme d’optimisation 
assure la sélection adéquate de ces pondérations, on fait un appel à une méthode d’approximation 
basée sur l’approche d’identification fréquentielle avec laquelle l’approximation a été assurée à 


partir de la minimisation de l’écart en module de ces deux transferts. 


La validation de cette nouvelle stratégie a été montrée sur le même système physique cité 
dans les chapitres précédents et les résultats fournis par le contrôleur H,, robustifié par la 
robustification basée sur les pondérations fractionnaires ajustables sont comparés dans les plans 
fréquentiel avec ceux fournis par les deux contrôleurs H,, à savoir, H,, robustifié synthétisé à 
partir des pondérations irrationnelles ajustables et H, primaire synthétisé à partir des 
pondérations irrationnelles fixes. A partir de ces résultats, on a pu conclure que l’introduction 
des pondérations fractionnaires ajustables dans le problème de synthèse par la méthode A, offre 
la possibilité de satisfaire plusieurs condition, en particulier la condition RP.s En conséquence, 


cette amélioration nous ont conduit à atteindre l’objectif de ce travail. 
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Chapitre 4 


Application du contrôleur robustifié sur un 
moteur à aimants permanents 


4. Introduction : 


En général, le moteur synchrone à aimants permanents (PMSM) est un dispositif très utilisé 
dans la plupart des applications industrielles tels que: la traction, l’aéronautique, la 
robotique,..….etc. Ceci est dû à plusieurs spécificités caractérisant les propriétés électrique et 
mécanique de cette machine on l’occurrence: l’absence des pertes rotoriques, l’obtention des 
réponses rapides en courant et en couple qui sont souvent liés à des inductances faibles, la 
production d’un couple massique élevé, une inertie très faible et une grande capacité de 
surcharge, ...etc. De ce fait, l’utilisation de cette machine est indispensable dans l’industrie, ce 
qui nécessite la synthèse d’un contrôleur robuste garantissant des performances dynamiques et 


statiques avec des marges de sécurité très élevées. 


Désormais, la synthèse du contrôleur souhaité est basé sur le développement d’un modèle 
mathématique adéquat décrivant le comportement réel du PMSM de façon les plus précises que 
possible. En général, Ce comportement est souvent régit par un modèle MIMO dont ses 
grandeurs d’entrées-sorties sont fortement couplées, ce qui peut compliquer sa phase de synthèse 
notamment par les méthodes conventionnelles de commande. A cet effet, il existe de nombreuses 
techniques dédiées à la commande de cette machine. Dans ce contexte, la commande par flux 
orienté (ou en anglais Field Oriented Control FOC) permet d’éliminer le couplage entre les 
grandeurs électriques du stator et celles du rotor du PMSM. Ce découplage permet de rendre 


cette phase de synthèse similaire que celle utilisée dans la machine à courant continu. 


En revanche, la commande de la vitesse rotorique ainsi que celle du courant statorique via des 
contrôleurs conventionnels tels que PI et PID dans la commande FOC deviennent insuffisantes 
pour répondre aux performances exigées du fait que les paramètres du modèle simplifié du 


PMSM sont souvent entachés par des incertitudes de modélisation. 


Récemment, diverses stratégies avancées de synthèse telles que, la synthèse par les réseaux 
de neurones, la synthèse à base de la logique floue, la synthèse neuro-flou adaptatif et d'autres 
approches ont été proposées dans la littérature pour pallier les problèmes susnommés et assurer 
la robustesse de RP avec une marge de sécurité élevée lorsque les paramètres du modèle 


changent dans une large gamme. Malheureusement, l’analyse de robustesse du système bouclé 


par ces contrôleurs n’est pas possible dans le plan fréquentiel en raison de l’absence 
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d'expressions analytiques décrivant le transfert de ces contrôleurs. Par conséquent, les valeurs 
singulières maximales de différentes sensibilités fournies par le système bouclé correspondant ne 
peuvent pas être tracées dans le domaine fréquentiel. Ce constat pousse à l’utilisation des 
méthodes très efficaces existant dans la commande robuste dans laquelle le contrôleur à 
synthétiser permet de quantifier à quelle marge de robustesse RP le système bouclé peut 
conserver les performances souhaitées en présence des incertitudes de modélisation. Sachant que 
l’amélioration de cette marge faisant l’objet de la commande du PMSM, la robustification du 
contrôleur H, primaire basée sur les pondérations fractionnaires ajustables sera adopter pour 
fournir un contrôleur H,, robustifié assurant le bon compromis « poursuite-régulation » de la 
partie mécanique du PMSM et cela malgré l’existence des anomalies produites via la commande 
de son compartiment électrique par des contrôleurs conventionnels. Par conséquent, le contrôleur 
H> robustifié à la capacité de fournir des meilleures robustesses que ceux synthétisés par les 


méthodes de synthèse basées sur les pondérations irrationnelles fixes ou ajustables. 


Dans cette optique, l’objectif principal est de comment assurer la sélection automatiques de 
ces pondérations fractionnaires et de comment les introduire dans le critère de sensibilité 
pondérée pour que la méthode H,, puisse d’élaborer un contrôleur capable de fournir une bonne 


amélioration de la marge RP du contrôleur H,, primaire. 


4.1. Description du système PMSM : 


Le moteur synchrone présenté par la Fig.4.1 est composé par d’une partie mobile, appelée 
rotor ou l’inducteur, et une partie fixe, appelée stator ou l’induit, portant des enroulements. La 
zone localisée entre ces deux parties est appelée entrefer [82-86]. La Fig.4.1. représente une 


coupe du PMSM. On obtient donc : 


Aiïmant permanent 
surfacique 


Fig.4.1. Vue en coupe du PMSM. 
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Dans une machine synchrone triphasée, le stator est constitué de trois enroulements identiques 
décalés entre eux par un angle de 120°. Ces trois enroulements sont installés dans les encoches 
du circuit magnétique fixe ou ce dernier est feuilleté afin, d’un côté, de réduire les courants de 


Foucault et d’un autre côté, de limiter les pertes dans le fer. 


Dans la partie du rotor, les enroulements est généralement parcourus par un courant continu 
(dans le cas de la machine à rotor bobiné). Par contre, dans le cas du PMSM, le rotor n’a pas 
besoin d'un approvisionnement extérieur pour être excité. En revanche, ses aimants permanents 
génèrent le champ magnétique ce qui ne nécessite aucune alimentation en tension ou en courant. 
Ceci présente l’avantage d’éliminer, d’une part, les pertes rotoriques liées aux enroulements et 
d’une autre part, de réduire le coût de réalisation et d’implémentation dans les chaines de 


production. 

Dans l’industrie, on peut distinguer deux types de PMSM, à savoir : PMSM à rotor avec des 
aimants montés en surface, qui se caractérisent par un très faibles taux de saillance ( & 1), et 
PMSM dont le rotor est constitué avec des aimants enterrés offrant un entrefer variable donc un 
coefficient de saillance supérieur à l’unité ( > 1} Dans ce présent travail, nous verrons que la 
synthèse du contrôleur H.. sera effectuée pour le premier type dans lequel la vitesse mécanique 


Nn dépende seulement de la tension statorique appliquée aux enroulements ainsi que du nombre 


de paires de pôles n,,. 


n, => (4.1) 


Avec : 
e {fh: représente la vitesse de rotation de la machine 
e _] : représente le nombre de paires de pôles 


e _w, : représente la vitesse électrique du rotor 


4.2. Modélisation de la machine synchrone à aimant permanent 


En général, le système PMSM est souvent composé d’un stator dont les enroulements sont 
connectés en étoile et un rotor avec n, paires de pôles dont l’excitation est créée par des aimants 
permanents surfacique. Afin de développer un modèle mathématique assurant la bonne 
description de la dynamique réel du PMSM on doit admettre les hypothèses simplificatrices 


suivantes : 


e Le système de tension est symétrique est équilibré 


e La distribution de la force magnétomotrice créée par le stator est sinusoïdale. 
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e La densité de courant est uniforme dans la section des conducteurs. L’effet de peau est 
négligé. 

e Le circuit magnétique n’est pas saturé. 

e L'effet de l’enroulement amortisseur au rotor doit être négligé. 

e Les pertes causés par les courants de Foucault et ceux causés par l’hystérésis doivent 
être négligés. 

D'où, la Fig.4.2 représente le schéma équivalent du système PMSM dans le repère 


triphasé (a — b — c). On obtient : 


Axe de référence 


Fig.4.2. Schéma équivalent du PMSM dans le repère (a — b — c) 


Les équations des tensions statoriques dans le repère (a — b—c) liées au stator sont 


exprimées comme suit : 


; dha 
Va = Roi + 
: da 
Vy = Rs-iy + (42) 
-: 1 dhc 
V = R..i + de 


Avec : 
e _R, : représente la résistance de l’enroulement statorique ; 
e inc : représente les courants de trois phases statoriques à, b, c ; 
e Parc représente les flux produits par les trois phases statoriques à, b, c ; 


e nc représente les tensions des phases statoriques. 
Les équations du flux produit dans chaque enroulement statorique sont données par les 
équations suivantes : 


> Le flux dans la phase « a », noté D, est exprimé par: 


Pa = Paa + Pba + Pca + fa (43) 
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Avec : 
e  Daa : représente le flux propre de la phase « a » 
e _ Ppa : représente le flux mutuel de la phase « b » 
e _ Dca : représente le flux mutuel de la phase « c » 


e__ Dra: représente le flux mutuel de l’aimant sur la phase «a » 


En remplaçant chaque flux par sa valeur correspondant, il en découle immédiatement les 


expressions des flux statorique de trois phases « a », « b » et «c ». On obtient : 


Pa = Ls-ia + Mss-ip + Miss. ic + Dfa 
Pp = Mss-ia + Ls-ip + Mssie + Dfa (4.4) 
ce — Ms: la + Ms: Üb + Este + fa 
Où : 
e _L, : représente l’inductance propre d’une phase statorique. 


e M, : représente l’inductance mutuelle entre les trois phases du stator. 


4.2.1. Transformation de PARK 


Les deux systèmes d'équation (4.2) et (4.4) sont fortement non linéaires et couplés entre eux. 
En conséquence, la résolution d'un tel système d'équation est délicate. Pour contourner ce 
problème, il existe, dans la littérature, une transformation appelée PARK P(8) qui consiste de 
projeter toutes les combinaisons présentées dans le repère (a — b — c), on l’occurrence, les 
tensions triphasées, les flux et les courants triphasés, en vue d’obtenir d’autres combinaisons 
fictives représentées dans le repère (d — q). Entre autres, cette transformation est interprétée 
physiquement, comme étant une substitution des enroulements immobiles « a », « b» et «c » par 
des enroulements en directe « d » et en quadrature « q » tournant avec le rotor. En fait, cette 
transformation permettant de faciliter la modélisation et la commande du système PMSM par les 
méthodes existant dans la littérature. La Fig.4.3 représente le schéma équivalent du PMSM dans 


le repère orthogonale de Park : 


47 
| 7 
St" 
FN 
dd 
= —_—__ 
en = 


Fig. 4.3. Schéma équivalent du PMSM dans le repère (d — q) 
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La transformation de PARK, notée P(8), est définie comme suit : 


Xà 
[Kago] = (Xal ; [Xa Xa Xol = P(B)EXe X5 Xel' 


X est un vecteur définissant peut être, le courant, le flux, la tension, ...etc. 


cos(6) cos (e — Æ) cos (o + æ) 
P(8) = E _sin(8) —sin (e - 2) —sin (e + =) 


V7 Vz Vz 


(4.5) 


(4.6) 


0 représente la position du rotor. La transformation inverse de PARK, notée P71(8) , est 


exprimée par la matrice ci-dessous : 


cos(8) —sin(0) _ 
1 .- Dan _ 35) 
P @)= cos (9 =) sin (0 =) 
2T . 27 1 
COS (e + Æ) —sin (o + æ) F 
Il vient alors : 
Xa 
L t 
LXavel = [al ; La Xo Xclf = P74(0)[Xa Xa Xo] 
Xo 


4.2.2. Modélisation de système PMSM dans le repère de PARK 


Dans le repère de PARK (d,q) tournant à la vitesse angulaire w, = 


dt 


d’équation donné par (4.2) peut être réécrit comme suit [82-86]: 


> Tensions statoriques, directes et en quadrature : 


d 

U; — Rs. ia Fees Or. Da 
$ dp 

U, = Rs.iq + + @r. Pa 


(4.7) 


(4.8) 


, le système 


(4.9) 


Les équations des flux statoriques s’écrivent, dans le repère de PARK (d, q), comme suit : 


> Les flux directs et en quadrature : 


fs = La.ia + 
a = La: iq 


° _ Da, PA représentent les flux statoriques dans le repère (d, q) ; 


(4.10) 
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e La, L, représentent les inductances cycliques statoriques d’axe direct d et en quadrature 


q respectivement. 


En remplaçant les flux par leurs valeurs dans les équations (4.9), on obtient donc : 


Ua = Rs-ia + La SE — @e La ig 
(4.11) 


: L Fe 
U, —= Rs. lq + LQ 2e + ©, (La. la + Pr) 


Avec : 

e U;, U, représentent les tensions statoriques dans le repère (d, q) ; 

° ia, 4 représentent les courants statoriques dans le repère (d, q) ; 

e_ Pa, DA représentent les flux statoriques dans le repère (d, q) ; 

+ _ représente le flux induit par les aimants permanents ; 

e La,L, représentent les inductances cycliques statoriques d'axe direct d et en quadrature q 


respectivement ; 


> Couple électromagnétique du moteur 
3 : E 
Com = >np-(Dr-ia + (La — La)ia-ia} (4.12) 


> La dynamique motrice peut être simplement décrite par l'équation : 


dAm 
JE + fem = Cem — Ct (4.13) 


Les trois systèmes d’équation (4.11), (4.12) et (4.13) montrent que le comportement réel du 
système PMSM est généralement modélisé par un modèle non linéaire. Comme ce dernier est 
caractérisé par des aimants montés en surface (1.e., L; = L,). D’où, le modèle décrivant ce 


système peut être illustré par la Fig.4.4. On obtient : 


Fig. 4.4. Schéma bloc décrivant le modèle du PMSM. 
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Dans ce travail, l’étape de synthèse du contrôleur stabilisant la boucle de commande nécessite 
a priori de définir les paramètres du modèle nominal du PMSM. Le Tab.4.1 récapitule la 
signification et la valeur correspondante de chaque composant intervenant dans la définition de 


ce modèle de synthèse. 


Tab.4.1 récapitule la signification et la valeur correspondante des différentes composantes du PMSM 


Paramètres significations Unités-Valeurs 

U; Tension statorique dans l’axe d 14 

U, Tension statorique dans l’axe q V 

la Courant statorique dans l’axe d A 

lg Courant statorique dans l’axe q A 

LA Inductance statorique de l’axe d 8.5 x 10° H 
Lo Inductance statorique de l’axe q 8.5 x 10H 
Pa Flux statorique dans l’axe d V.s 

Ÿa Flux statorique dans l’axe q V.s 

®r Flux des aimants 0.175 V.s 
R; Résistance statorique 0.2Q 

Or Vitesse électrique du rotor tr/min 

Np Nombre de paires de pôles 4 
D Vitesse mécanique du rotor tr/min 

] Moment d’inertie 0.089 kg.m° 
fe Coefficient de viscosité 0.005 N.m.s 
Ce Couple électromagnétique N.m 

Ce Couple de charge N.m 


4.2.3. Linéarisation du modèle de PMSM : 


D’après les équations (4.9), (4.13), la dynamique d’interaction produite via le couplage des 
combinaisons électriques présentées dans les deux d et q doit être compensée afin de linéairiser 
le modèle non linéaire du systèmePMSM. Cet objectif peut être abouti en utilisant le proncipe de 
la commande FOC dans laquelle le couplage de ces combinaisons peut être résolu par la synthèse 
de deux contrôleurs de courant, notés PID, et PID,. Dans ce cas, le rôle du contrôleur PID, est 
de pouvoir atténuer, autant que possible, l’erreur générée par le courant statorique directe de 


référence iq,., # © celui mesuré à partir du systèmePMSM. De ce fait, un forçage de i4,, — 0 


permet d’assurer le découplage désiré qu’à partir duquel la vitesse mécanique 9, est asservie 
seulement dans l’axe q par un réglage optimal du couple électromagnétique C: ou le courant 


d’axe en quadrature i,. L’idée de cette stratégie e commande est illustrée dans le schéma de la 


Fig.4.5 comme suit : 
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Three — phase 
; inverter 
PID; current 
controller _ aq € ù 
PWM 
Generator 


| 


| PID, current 
controller 


Fig.4.5. Commande de vitesse du PMSM basée sur la stratégie de commande à flux orienté 


D’après la Fig.4.5, le modèle linéaire du PMSM est régi par le système d’équation simplifié 
suivant : 
_ : diq 
Uy = Rs-ig + Lo as Ppnp.Q 
3 : 
Cem = 7° Np- Pr. la (4.14) 


dOm 
+ fe Qm = Cem — Ce 


J. 


En outre, l'avantage principal de la stratégie du FOC réside dans la possibilité de synthétiser 
un pré-compensateur SISO, ce qui rend cette stratégie de commande très attrayante pour les 
applications industriels. L'objectif de notre synthèse est de développer un contrôleur H,, robuste 
assurant le bon compromis « poursuite-régulation » du système bouclé malgré les incertitudes de 
modélisation et/ou l’existence d’une dynamique d’interaction dans le compartiment électrique du 
PMSM qui est probablement causé par un mauvais réglage des deux contrôleurs de courant cités 


précédemment. 


Tout d’abord, on suppose, dans l’étape de synthèse du contrôleur H.,, , que le couple de charge 
C+ est nul dans l’équation (4.14), puis, la fonction Linmod du Matlab® on est appliquée sur le 
schéma bloc présenté dans la Fig.4.6 pour déterminer la représentation d’état du transfert 


nominal Gy(s). On obtient: 


Fig.4.6.Schéma block utilisé pour obtenir la fonction de transfert du régime nominal G\(s) 
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D’après la Fig.4.6, le contrôleur PID, doit minimiser l’écart entre le courant d’axe en 
quadrature mesuré i, et sa référence à, > Ce dernier présente plus tard le signal de commande 


fourni par le contrôleur H., à synthétiser. 
Notons ici que la synthèse du contrôleur PID, est assurée par la fonction Rltooldu Matlab®. 
Cette dernière fournie le transfert ci-dessous: 


4557(s+45.74) 


PIDQ(s) T s(s+101 ) 


(4.15) 


De plus, le modèle nominal est assuré par la fonction Linmod du Matlab® . on obtient 


donc : 


6.325X10%(s+45.74) 
(s+49.87)(5+0.054)(52+992.2 xs+5.116X10°) 


Gn(s) = (4.16) 


Il est évident que les paramètres de ce modèle dépendent des paramètres du contrôleur de 
courant PID, et ceux présentés dans le Tab.4.1. En conséquence, deux inconvénients majeurs 
peuvent être contribués à la dégradation des performances de la stratégie du FOC. Le premier 
inconvénient réside dans la grande sensibilité de cette stratégie à l'effet indésirable provoqué par 
la fatigue du capteur de vitesse mécanique et/ou celui du courant d’axe en quadrature, tandis que 
le deuxième inconvénient réside dans la dégradation de la robustesse du système bouclé à cause 
de la variation des paramètres du modèle dans une large bande. De ce fait, pour contourner les 
deux problèmes précédents, la boucle externe de la stratégie FOC doit contenir un contrôleur 


plus robuste offrant une robustesse avec une marge de sécurité très élevée. 


4.3. Synthèse du contrôleur H, primaire : 


Le problème initial de sensibilité pondérée est formulé en utilisant les deux pondérations 


irrationnelles fixes Ws, (s) and W7, (s) suivantes : 


ns, — 1 
S Jos Op = 25 
Ws,(s) = ——, where Ms, 419 (4.17) 
é 195 *X10 * 
nr — 1 
is 40 hr = 400 
WC) = RE wherc À y 2 125 (4.18) 
400 0 ' 


en = 166 107 
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Ensuite, la fonction Hinflmi du Matlab® est utilisée pour résoudre ce dernier problème, ou 
lPindice de performance du contrôleur H, était fourni par Yopr = 0.9216. En effet, la 


représentation d’état du contrôleur H,, primaire K, (s) est obtenue comme suit : 


—119.14 15.27 —4638.9 —11886  —1.2038: 10° —4.6427 : 10° 
9553.3 —1775.4 3.795-10°9.72:10° 9.8482:106  3.7981 : 107 
4, =111622 —15.569 —211.96 —464.26 —4291.7 —18991 
id 40358 —0.61284 —-5.948  —417.1 —8516.3 —2450.6 
21.309 —0.95047 —855.66 —62034 —1.5671 : 10° 68799 
7.0595 —7.5889  2316.3 38011  1.2289:10°  —1.4153 : 10° 
5.2409 x 104 —39.652 F 
1.5057 x 1074 2.9869 
—2 _ 
ne 0 land) 1664 e19 
— 39.289 —38472 
57.504 — — 14838 x 10° 


Il est clair que le contrôleur H,, primaire K, (s) est d’ordre 6 et sa matrice de transfert peut 


être déterminée par la relation Ko (S) = Co (s X 1 — A5) 1.Bo 
4.4.Synthèse du contrôleur H, robustifié: 


D’après la Fig.4.5, il est clairement perceptible que la condition de robustesse sur les 
performances initiale RP, n’est pas satisfaite par le contrôleur H, primaire. Pour améliorer cette 
marge de robustesse, nous allons introduire les pondérations fractionnaires ajustables dans la 
formulation du problème de la robustification proposée dont la sélection automatique de chacune 
est assurée par l’algorithme PSO et l’approximation en filtre d’ordre entier est assurée par 
l’identification fréquentiel. Pour atteindre ces objectifs, l’espace de recherche limitant le vecteur 


x. est choisi comme suit : 


1.15 Ms 2.50 
1.00 ns 3.00 
25.00 OB 35.00 
1076 alé" 1074 


* 


1.20 |" | Mr] | 2.00 (28) 
1.00 nr 3.00 
200 | [ol | 400 
105) Le, Li10-4 
Xcmin Xc Xmax 


Les paramètres de réglage de l’algorithme PSO ont été choisis comme: 


=" Nombre maximum des itérations : Unax = 20 
" Taille de l’essaim : n, = 20 


»+ + Co CO mi x 
"Facteur d'inertie Co = Cour — (ne nn) OÙ Comax = 09 Comm = 04etl=1,…, Enax- 


Cmax 
= Taux d'apprentissage cognitif: c, = 1.8 


=" Taux d'apprentissage social :Cc, = 1.8 
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La fonction objective a été évaluée à 1000 fréquences réparties uniformément sur une plage 


de fréquence valant [107*, 10*%] radian/seconde. 


Par conséquent, la solution optimale fournissant la bonne minimisation de cette fonction 
objective est ensuite utilisée pour déterminer les paramètres des deux pondérations fractionnaires 
désirées Ws - (w,x)etWr : (w,x). Entre autres, l’algorithme PSO est compilé 30 fois avec 
différentes populations initiales de x, € Ê? À ] . Par conséquent, la solution optimale est 

min max 


obtenue à l’itération £,y4x = 20, aboutissant les deux pondérations fractionnaires ajustables 


suivantes : 
ns = 1.229 
+ 1386 "7 w* = 281386 
#) = (NES — 48. 
MEGA : =) Por Ms = 2.251 (421) 
es = 1.9 x 1074 
s à 2.3784 a nn 
30623784 Dr = 
Wr,(s,x”) = res > Ou M, = 1.88 (4.22) 
— ———,5s+1 
Éd é5=29%10 


Ces deux transferts fractionnaires sont ensuite approximés en deux filtres d’ordre entier en 


exigeant les deux précisions suivantes : 


= Précision souhaitée pour approximer Ws, (w, x ): AW, = 1072 


= Précision souhaitée pour approximer Wr : (w,x ): AWr = 1072 


Finalement, l’approximation des deux pondérations fractionnaires précédentes fournie les 


deux filtres d’ordre entier ci-dessous: 


__ 0.444(s+57.01)(5+29.46)(s+6.794)(s+0.9529)(s+0.122) 

LACS (s+22.28)(5+4.476)(s+0.5879)(s+0.0795)(s+0.024) (4.23) 
___ 34482(5+763 )(s+217 )(s+562.9)(s+3.024)(52+406.35+42  x104) 

Wr, (s) 7 (s+1.167x104)(5+3561)(5+936.3)(5+3.025)(s2+5.1xX1045+663 X104) (4.24) 


La Fig.4.7 (côté gauche), compare les réponses  fréquentielles des deux 


pondérations Ws : (s,x*) et Ws,(s) dans le plan de Nyquist. Cependant La Fig.4.7 (côté droit), 


compare les réponses fréquentielles des deux pondérations Wr - (s,x*) et Wr., (s). On obtient : 
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— Optimal W. x) 


= 
D 


— Optimal MAS) …_ Appronimated W, {s) 


Imaginary Axis 


— - Approximated W, nl 


Imaginary Axis 


4 05 0 05 1 15 2 25 3 3 
Real Axis xt0! 


Fig.4.7. Comparaison entre ( Ws.Co,x°); Ws, (w)) présenter dans le côté gauche, et entre 


( Wr,(@,x*) ; Wr, (w)) présenter dans le côté droite. 

D’après la Fig.4.7, il est facile d'observer que les deux réponses fréquentielles des 
pondérations fractionnaires ajustables Ws,(s, x) et Wr ACTES) sont étroitement rapprochées 
celles fournies par Ws (s) etWr.(s) respectivement. En outre, le problème de sensibilité 
pondérée est résolu par la fonction Hinflmi du Matlab® fournissant donc un niveau de 
performance Yope = 0.6946 et un contrôleur H,, robustifié d’ordre 15. Par conséquent, les deux 
conditions de robustesse, à savoir RS et NP sont examinées avant et après la phase de 


robustification du contrôleur H,, primaire. 


A cet effet, la Fig.4.8 présente le tracé des lieux des valeurs maximales singulières maximales 
des deux sensibilités directes S4,(s) et Sa, (s). Ces dernières sont également comparées, dans 
les basses fréquences avec les deux pondérations inversées Ws*(s) et Ws,*(s) (voir la Fig.4.8 
côté gauche). De même, la Fig.4.8 présente le tracé des lieux des valeurs singulières maximales 
des deux sensibilités complémentaires 5, (s) et S- (s) qui sont également comparées dans les 
hautes fréquences avec les pondérations inversées Wir" (s) et Wr"(s) (voir la Fig.4.8 côté droit). 


On obtient : 


Page| 117 


CH.4: APPLICATION DU CONTRÔLEUR ROBUSTIFIÉ SUR UN MOTEUR À AIMANTS PERMANENTS 


Dee DENTS TENTE Ter TENTE ATT] 
(A fi +, 
407 LL 
A 
# 20r | 
ï u) 20! 
! 
H 0" 14 
5 $ 0! 
| | e | 
s 40 
! d 3 40 ——— 
: 3 
La LA 13 IS 1 
Le ; m so au JW, 
ü 40 LH E ; 
| sh LE] 
j " LE UMTS y TN, 
PATES = MM = tI 1 
#0! AU pl Fer | 
40 


10 40 402 407 40 two! #02 10 


Frequency (rad/sec. Frequency (rad/sec.} 


Fig.4.8. Comparaison entre (Sas (w); Sa, (w); Ws'(w); Ws;' (&)) présenter à gauche, et entre 
(se (); S,(w);: Win (w): Wr* (w)) présenter à droite 
Il est évident remarquer que la meilleure marge NP est obtenue quand les valeurs singulières 
maximales de la sensibilité directe sont faibles autant que possible notamment en basses 
fréquences. Par ailleurs, la meilleure marge de robustesse RS est obtenue quand les valeurs 
singulières maximales de la sensibilité complémentaire sont faibles autant que possible 


notamment en hautes fréquences. 


A partir de la Fig.4.8 (côté gauche), il est clair que le contrôleur H,, robustifié à la capacité 
d’assurer une meilleure marge de robustesse NP en comparaison avec celle assurée par le 
contrôleur H, primaire, en particulier la plage de fréquence w € [107*#, 20] radian/seconde. 
Ainsi, le lieu de Gruxl Sa .(w)| se situe en dessous de —80 dB, cela signifie que la perturbation 


est atténuée plus de 10000 fois à la sortie du système. 


De plus, à partir de la Fig.4.8 (côté droit), on peut remarquer que le contrôleur H, robustifié 
assure la meilleure marge RS en comparaison avec celle fournie par le contrôleur H, primaire. 
Encore, le lieu de Oéaxl Sc, (w)| est inférieure à—60 dB dans les fréquences au-dessus de w = 
4220 radian/seconde, cela signifie que l’effet indésirable provoqué par la dynamique négligée 


des modes inconnus est atténué plus de 1000 fois à la sortie du système. 


Finalement, on peut conclure que le contrôleur AH, robustifié a assuré non seulement un 
meilleur compromis de robustesse NP/RS , mais aussi une meilleure marge de RP, ce qui est 


illustré dans la Fig.4.9. 
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RP condition 
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Fig4.9. Satisfaction de la condition de RP par le contrôleur H, robustifié proposé 


Pour confirmer les résultats obtenus précédemment dans le domaine temporel, 
l’environnement Power Sim® est ensuite utilisé pour comparer les réponses temporelles 
fournies par les deux contrôleurs H, primaire et robustifié. Dans cette optique, ces deux 
contrôleurs sont connectés simultanément avec le modèle nominal du PMSM et celui existant 
dans la boîte à outils Power Sim®. Il convient à noter que la simulation est assurée par le choix 
du pas de discrétisation T, = 20 microsecondes. D’où, les deux systèmes bouclés, à savoir : le 
système PMSM commandé par les deux contrôleurs H, primaire et H,, robustifié et le modèle 


linéaire nominal du PMSM contrôlé par les même contrôleurs précédents, sont montrés par la 


Fig.4.10. On obtient : 


- - Hinfinity_contrloller 


PMSM linear model 
VECT 


Gaint 
Model and 


process outputs 


Reference 


Control signal Gain2 


powergui Hinfinity_contrloller. 


Process output 
andreference 


iq_ref 


DC Voltage 
Source <Rotor angle thetam (rad)> 


c = Load torque 
<Rotor speed wm (rad/s)> 


<Electromagnetic torque Te (N*m)> 


LL 1 
Measures Actual PMSM EE 


3-phase inverter <Stator current is_a (A}> 


Fig.4.10. Schéma block du PMSM et leur modèle linéaire avec les contrôleurs H,, primaire et robustifié 
A partir du schéma bloc de la Fig.4.10, nous allons examiner la conformité de l’étape de 
modélisation du comportement réel du système PMSM autour de son point de fonctionnement en 


utilisant le modèle nominal simplifié cité précédemment. A cet effet, l’entrée de référence est 


Page| 119 


CH.4: APPLICATION DU CONTRÔLEUR ROBUSTIFIÉ SUR UN MOTEUR À AIMANTS PERMANENTS 


excitée par un signal bruité pseudo-aléatoire variant au voisinage de la plage de fonctionnement 
du PMSM. Cependant, l’entrée du couple de charge est supposée nulle (démarrage à vide). La 
Fig.4.11 représente une comparaison entre la sortie fournie par le système bouclé par le 
contrôleur H., robustifié, qui est connecté avec le modèle nominal simplifié du PMSM, ainsi que 
la sortie du système bouclé par le même contrôleur précédent, qui est connecté avec le bloc 


Simulink du PMSM existant de l’enivrement Power Sim®. On obtient : 


410 1— 
360 F 
310F 
260 F 
210 Ë 
160 F 
110 E 
60 

10 F 

-40 F 

-90 F 

-140 F 
-190 F 
-240 F 
-290 F 
-340 F 
390 L——1 je 5 pe à + 
O0 04 08 12 16 2 


Mechanical speeds (rad/sec.) 


1 1 1 L 1 1 


n L 1 1 n L n L 
2.4 2.8 3.2 3.6 
Time (sec.) 


4 4.4 4.8 5.2 5.6 


Fig.4.11. vérification de l'exactitude de l'étape de modélisation du comportement réel du PMSM par le modèle 
linéaire simplifié 

Maintenant, pour examiner le compromis de robustesse « poursuite-régulation », le schéma de 

la Fig.4.10 est utilisé pour les deux contrôleurs : H, primaire et H, robustifié. Dans ce cas, les 

deux entrées exogènes précédentes sont modifiées de telle sorte que la première entrée, qui 

représente la vitesse mécanique de référence 97, fs exprimée en tours/minutes (r.p.m), est 


imposée par le signal suivant: 


50 
25 


:0<t<0,5 
:05<t<1 


0 
0.) = 4-25 
—50 
—40 


:1<t<12 

:12<t<1.7 
17 Et< 22 
ADSL LT] 


(4.25) 


—30:27<t<4 
mn — 7. 
Mref 30 


0 d’où, la référence citée ci-dessus peut être 


En tenant compte la relation 9 Mref? 


transformée à celle donnée par (4.26), qui est exprimée en radians/secondes (rad/sec.) comme 


suit: 


47746: 0<t<0.5 
238.73 : 05 <t<1 
0:1<t<12 
—238:73:12<t< 17 
—477.46 : 1.7 < t < 2.2 
39179: 22<t<27 
—286.48 : 2.7 < t < 4 


a (4.26) 


Mref — 
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En outre, le filtre du premier ordre F(s) = est inséré dans la boucle pour rendre la 


0.01.5+1 


référence ci-dessus plus proche à celle utilisée dans la réalité. De plus, la deuxième entrée 
exogène représente le couple résistant qui est imposé par un échelon unitaire d’amplitude C; = 
15 N.m appliqué à partir de l’instant t = 3 secondes. La Fig.4.12 représente les vitesses 
mécaniques données par les deux contrôleurs H., primaire et robustifié. Cependant, la Fig.4.13 
représente les deux courants en quadrature fournis pour les deux contrôleurs suscités. On 


obtient : 


Mechanical speeds 


— — Reference speed 
—— Robustified H-infinity controller 
— — Primary H-infinity controller 


0 05 1 1.5 2 25 3 3.5 4 
Time (sec.) 


Fig.4.12. les vitesses mécaniques données par les deux contrôleurs H, primaire et robustifié 


275 1 ———— —— 


Robustified H-infinity controller 


770 K 


— — Primary H-infinity controller 


570 195 


370 
115 


35 


Quadrature current 


H 
1 
-30 di 


-230 F 


Quadrature current 


-430 - 


e] 
e] 
() 
TS 
Es 
(6) 
Lo 
N 
on 
œ 
(20) 
() 
EN 


770 
620 
470 
320 
170 

20 | 
-130 | 
-280 | 
-430 


- — Primary H-infinity controller a 
—— Robustified H-infinity controller | 1 


Fee 


Quadrature currents 


2 
Time (sec.) 


Fig.4.13. les courants en quadrature fournis par les deux contrôleurs H,, primaire et robustifié. 


Page| 121 


CH.4: APPLICATION DU CONTRÔLEUR ROBUSTIFIÉ SUR UN MOTEUR À AIMANTS PERMANENTS 


Selon les Fig.4.12 et Fig.4.13, il est facile de voir qu’une bonne dynamique de suivi de la 


vitesse de référence est bien assurée par les deux contrôleurs proposés. Cependant, la meilleure 


dynamique d’atténuation du couple résistant est celle fournie par le contrôleur H,, 


dernier fournit également un courant en quadrature i, très raisonnable, délimité 


robustifié. Ce 


par les deux 


bornes —125A et 270 À. En contre partie, le signal de commande fournit par le contrôleur H,, 


primaire est caractérisé par un dépassement très élevé, ce qui est souvent intolérable 


éventuellement dans le cas des applications industrielles. De ce fait, on constate 
présence des pics du courant en quadrature i, très élevés, délimités par les deux b 


et 770 À. Ceci est causé par les changements rapides de l’entrée de référence. 


également la 


ornes —4304A 


La Fig.4.14 représente les courant statoriques triphasés du PMSM, tendis que la Fig.4.15 


représente les tensions statoriques triphasées qui alimentent l’onduleur triphasé. On obtient 


donc : 


7 " L 
200 É 06 Stea 
<Stator current is, (A)> 2 04 
= 
<Stator current is, (A)}> Oo 
150 6 (A) 5 02 
<Stator current is, (A}> ü 
à 0 
2 8 
ô $ 02 
+ rs 
3 o 041 
£ E 06 ! 
ü 1 L L L L 1 L 
2 3925 393 3935 394 3945 395 3955 396 3.965 
a 
© - 100 T T 
2 
= 50 Transient state L 


_200 ñ ñ ñ 1 1 1 1 1 i 1 1 1 400 


Three phase stator currents 


1 
0 03 06 09 12 1.5 1.8 21 24 27 3 33 36 394 


Time (sec.) 


0.02 0.024 0.028 0032 0.036 0.04 


Time (sec.) 


0.044 


0.048 


Fig.4.14. Les courants statoriques triphasés du PMSM 
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Fig.4.14. Les tensions statoriques triphasées qui alimentent l’onduleur triphasé 
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4.5. Conclusion : 

Dans ce chapitre, nous avons validé la robustification proposée pour le contrôleur H,, 
primaire dans laquelle la formulation du problème de sensibilité pondérée est assurée par des 
pondérations fractionnaires ajustables. Cette robustification permettant d’améliorer les marges de 
robustesse NP, RS et RP. Ces objectifs ont été atteints via la sélection automatique de ces 
pondérations tout en satisfaisant certaines règles de réglage existant préalablement dans la 
littérature. Cette robustification proposée est basée sur les étapes suivantes : Tout d’abord, 
l’algorithme PSO est appliqué pour la recherche des bonnes pondérations fractionnaires 
ajustables dans un espace prédéfinit par l’utilisateur. Ensuite, l’approche d’identification 
fréquentielle est appliquée pour fournir des filtres d’ordre entier rapprochant ces pondérations 
fractionnaires dans une plage de fréquence précise. Ces derniers permettant d’augmenter le 
modèle nominal de synthèse ce qui permet de mettre à jour l’ancien problème de sensibilité 
pondéré. Ce dernier est résolu par la méthode A., et la solution obtenue est utilisée pour la mise à 
jour du transfert de contrôleur H,, primaire. Les sensibilités fournies par le système bouclé 
correspondant donnent la possibilité de mettre à jour la marge de robustesse précédente RP avec 
laquelle le critère d’arrêt de l’algorithme de robustification proposée dépend de cette nouvelle 


marge de robustesse. Ce processus se répète jusqu'à l’obtention du contrôleur H., robustifié. 


En fin de compte, les résultats de la simulation obtenus montrent l’amélioration évidente de la 
marge RP par le contrôleur H,, robustifié dans lequel deux gains significatifs ont été assurés. Le 
premier gain réside dans l’amélioration de la dynamique d’atténuation du couple de charge 
tandis que le second gain réside dans l’obtention d’une commande économique et moins 
couteuse par rapport à celle fournie par le contrôleur H,, primaire. D’autre part, il est également 
clair que d’autres améliorations du schéma de contrôle proposé nécessiteront une structure de 
contrôleur fixe qui devrait être donnée avec un ordre très réduit et qui ne doit pas dépendre par 


l’ordre du système augmenté. 
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Conclusion Générale 


Dans ce travail, la robustification d’une classe des contrôleurs H,, primaires a été proposée 
afin de sécuriser les systèmes bouclés en présence des effets non souhaités causés par les 
incertitudes de modélisation, l’influence des dynamiques négligées de quelques modes 
inconnus décrivant le comportement non linéaire du système physique, l’effet des bruits de 
mesurés qui souvent produit pat la fatigue ou la défaillance des instruments de mesure 


installés dans la boucle de commande...….etc. 


La robustification proposée a visé l’amélioration de la marge de robustesse sur les 
performances robustes, qui est fournie préalablement par les contrôleurs H,, primaires. Cette 
amélioration offre la possibilité d’assurer parallèlement et avec succès l’amélioration des 
plusieurs d’autres marges de robustesse à savoir: la robustesse sur les performances 
nominales du système corrigé, la robustesse de la stabilité robuste en présence des incertitudes 
multiplicatives, le compromis entres ces deux dernières robustesses. En effet, la 
robustification proposée à la capacité de produire des contrôleurs A, robustifiés assurant 
tous les objectifs susmentionnés et ce non seulement pour le régime de fonctionnement 
nominal du système physique, mais aussi pour tous les régimes perturbés potentiels qui 


peuvent apparaitre à cause des erreurs commises lors de la phase de modélisation. 


En général, la synthèse d’un contrôleur H, primaire nécessite a priori le choix adéquat 
assurant la sélection des pondérations appropriées. Ces dernières ont été exploitées pour 
augmenter le modèle nominal de synthèse et le résultat de cette augmentation fourni un 
modèle augmenté. Sachant que le transfert de ce modèle est interconnecté avec celui du 
contrôleur à synthétiser. D’où la transformation linéaire basse appliquée sur ce système 
d’interconnexion permet de formuler le problème de sensibilité pondérée. La solution 
optimale de ce problème faisant la représentation d’état du contrôleur H, primaire. Le 
système bouclé par ce contrôleur à la possibilité de fournir des sensibilités avec lesquelles on 
a pu faire l’analyse de robustesse dans le plan fréquentiel. Celle-ci est basée sur le tracé des 
lieux des valeurs singulières maximales de différentes sensibilités. Ces lieux sont ensuite 
comparés par ceux fournis par les transferts inversés des pondérations. Cette comparaison 
permettant de calculer les différentes marges de robustesse, en particulier, la marge de 
robustesse sur les performances robustes. Les étapes requises, pour développer le 


contrôleur H, primaire ont été détaillées et validées dans le premier chapitre de ce manuscrit. 


La validation a été effectuée à partir de la commande d’un système hydraulique dont le 
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comportement réel était donné par un modèle incertain entaché par des incertitudes 


multiplicatives amenant en sortie du modèle de synthèse. 


A présent, la sélection des pondérations est assurée manuellement. Celle-ci n’est pas 
souvent facile, notamment si le problème de synthèse est de satisfaire la condition sur RP. La 
satisfaction de cette condition avec une bonne marge de sécurité faisant donc l’objet de cette 
thèse. En effet, le problème majeur mentionné dans plusieurs travaux antérieurs, est la 
sélection automatique des pondérations satisfaisant toutes les spécifications imposées par le 


cahier de charges. En conséquence, l’objectif de notre travail devient plus précis. 


Nous avons proposé à cet effet une première solution basée sur la proposition d’une 
robustification du H,, primaire dans laquelle le choix des paramètres de chaque pondération 
irrationnelle ajustable était garanti par l’algorithme génétique. Les étapes requises, pour 
développer le contrôleur A, robustifié ont été détaillées et validées dans le deuxième chapitre 
de ce manuscrit. La sélection automatique était donc assurée par la recherche à des paramètres 
optimaux constituant chaque pondération. Cette recherche était réalisée dans un espace 
assurant des allures souhaitées de différentes sensibilités. Ces allures ont été imposées par des 
gabarits fréquentiels fixées a priori par quelques règles de réglage existant dans la littérature. 
Néanmoins, cette première robustification proposée à la capacité d’assurer une légère 
amélioration de RP sur une bande fréquentielle très étroite. Celle-ci est due à la saturation de 
quelques contraintes de bornes, en particulier celles fournies par certains ordres entiers portés 


sur quelques transferts des pondérations irrationnelles. 


Nous avons proposé à cet effet une deuxième solution fondée sur la robustification basée 
sur les pondérations rationnelles ajustables dont ses paramètres ont été optimisés par 
l’algorithme d’essaim particulaire. Les étapes requises, pour développer le contrôleur A,, 
robustifié ont été détaillées et validées dans le troisième chapitre de ce manuscrit. Dans cette 
proposition, les pondérations rationnelles adéquates doivent être approximées en des 
transferts d’ordre entier afin d’étendre l’application de la méthode H,, en cas de présence ce 
genre des pondérations. Pour approximer les transferts fractionnaires, on a fait un appel à 


l’identification fréquentielle basée sur la fonction du Matlab fitfrd. 


Finalement, le manuscrit est achevé par la validation de l’approche proposée sur un 
prototype très connu dans des revues de renommées internationales. C’est un moteur 
synchrone à aimants présenté dans plusieurs travaux antérieurs. Les résultats de simulation 


fournis par la robustification proposée sont très encourageants. Les améliorations de la marge 
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RP requises pour perfectionner cette idée peuvent ouvrir une fenêtre à d’autres éventuels 


candidats dans le futur proche. 
Perspectives : 


Comme perspective, nous conseillons les futurs candidats intéressés par cet axe de 
recherche de penser à introduire des structures fixes pour les deux contrôleurs : H, primaire 
et H, robustifié, car le degré de chacun représente la somme des degrés de tous les transferts 
utilisés dans la synthèse par la méthode H,, . En outre, d’autres perspectives peuvent être 


résumées dans les points ci-dessous : 


e Tenir compte des erreurs de modélisation structurées lors la robustification du 
contrôleur H, et amélioration de la marge RP face à ce type d’erreurs de 


modélisation. 


e Etendre la robustification sur les systèmes multivariables modélisés par des 


modèles linéaires à des paramètres variables (LPV model). 


e Proposer une méthode efficace et simple permettant d’approximer les transferts 
fractionnaires des pondérations rationnelles ajustables en des filtres irrationnels 


d’ordre très réduit 


e Proposer d’autres structures fractales dans les pondérations rationnelles ajustables, 
en particulier, imposer des ordres non entiers sur les opérateurs de Laplace et 


d’autres sur le transfert global. 
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ANNEXE A 


Annexe À : Transformation fractionnaire 
linéaire 


A.1. Définition d’une Transformation fractionnaire linéaire LFT : 


Soit la matrice complexe M donnée par la matrice de transfert ci dessous : 


(A.1) 


Soit la matrice complexe définissant les incertitude de modélisation A. L’interconnexion 


entre M et À, peut être représentée par la Fig. A.1 comme suit : 


Fig. A.1. LFT inferieure 


A partir de la Fig. A.1, il est possible de définir une LFT inferieure entre les deux transferts 


complexes M et À, telle que la relation ci-dessous : 
FM, À) = Mis + Mi2Ai( — M2) M1 (A.2) 


Où la matrice complexe (1 — M;,,A;) doit être inversible. De plus, le transfert reliant l’entre avec 


la sortie de cette interconnexion, représentée par Fig.1 peut être exprimé comme suit : 


G)=m(u)= [mel (us) a 
Sachant que : 
U = AY: (A.4) 


Supposons maintenant qu’on a une autre matrice complexe, notée À,,. Alors, l’interconnexion 


entre M et À, est présentée par la Fig. A.2 comme suit : 


Page | 135 


ANNEXE A 


Fig. A.2. LFT supérieure 


A partir de la Fig. A.2, il est possible de définir une LFT supérieure entre M et À, telle que la 


relation ci-dessous : 
ECM, A5) = M2 + M1 — Miidu) Mi (A.S5) 


Avec, la matrice (1 — M;,A,,) doit être inversible. Dans ce cas, le transfert d’entre/ sortie de 


cette interconnexion peut être exprimé comme : 


(on) = M (ue) = Fa ul … (A6) 


En Tenant en compte la relation ci dessous : 
U> — y: V2 (A.7) 


* Remarque : 
e La transformée linéaire fractionnaire inferieure F,(M,A,) fournie un transfert en boucle 
fermée dont l’entrée est w, et la sortie est z: . 
e La transformée linéaire fractionnaire supérieure, (M, À.) fournie également un transfert en 
boucle fermée dont l’entrée est w, et la sortie est z, . 
e La matrice M s’appelle « matrice du coefficient ». Elle représente le système à commander. 


e La matrice À , avec (A, ou A) représente le contrôleur ou la perturbation. 
A.2. Produit star de Redhef fer : 


Etant données les deux matrices de transfert P et telles que les relations suivantes: 


P P K K 
P = | 11 “| : K = | 11 . AS 
en ki Re 4) 


De l’équation (A.8), il existe un produit, appelé produit star de Redhef fer appliqué sur les 


deux transferts complexes P et K. D’où le résultat de ce produit peut être défini comme suit: 
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Fe(P,K11) Pi2Q — KiaP52) Ki 


. (A.9) 
K21( — Pa2Ka1) Por EiCK, P22) 


S(P,K)=P+0=| 


Ce dernier résultat peut être illustré par la Fig. A.3 comme suit: 


Si el 


PxK 
S2 €2 
Fig. A.3. Produit star de Redhef fer 
L’interconnexion du produit star de Redhef fer de P et K peut être représentée par la Fig. 


AÀ.4 comme suit : 


Fig. A.4. Diagramme fonctionnel de l’interconnexion entre P et K 


Si les matrices P et K n’ayant qu’un seul bloc, il est possible d’adopter les 


conventions suivantes: 
P'x Kis: = F(P,Ki1) et  P22 * K: = E,(K,P2) 


La fonction starp du Matlab peut être utilisée pour calculée le produit star de 
Redheffer comme suit : 
P*K = starp(P;K; Dy; Du) 
Avec : 
e  Dy: la dimension de y; 
e Du: la dimension de u 
e Propriété : 


Le produit star de Redhef fer est une opération associative vérifiant l’expression ci-dessous: 


(PxK)xL = Px(K+xL)= PxK+L (A.10) 
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A.2.1. Espace d’état du produit star de Redhef fer : 


Soient les deux réalisations minimales de P et X données respectivement par : 
A BB, B Ax  Bri  Bk2 
P=|C; Di D: ; K=|Cx1 Dxi1 Dri2 (A.11) 
C2 D1 Di Cx2 Dx21 Dka2 


La réalisation minimale du produit star de Redhef fer de P et K est donc donnée par : 


A B, B;  — 
P * K —= Ci Di Dis — é 2 (A.12) 
2 Dy1 D» 
AVE : 
— (£ B; ) . ee an _ (A + B2.R°Dyr1 Co BR Cri (A3) 
C> D; BK AK PT LE AK + Ban Die 
Re Fe B ) . Fe. _ [B1+ Bo. À Da Don B3.R 7 Dy12 (A.14) 
D; D;2 Bx1 Bx2 Bah "Ds Bx2 + Bah DD 
= _ (C1 Ds Dai Cri _ [Ci + Di2Dg11R "C2 Di. R Cyan 
— , = cr - (A.IS) 
C2 Di Dx21 Ck2 Dx21R7 C2 Cr2 + Dx21R7 D2Ck1 
= _ {D11 Di Dxa1  Dri2) _ [Dis + Di2Dk13R 7 Dos Di2. R°Dyns 
D = * = (A.16) 
D1 D22 Dx21  Dr22 Dx21R Dai Dx22 + Drk21R7*D22Dk12 


Sachant que les deux matrices R et À sont définies comme suit : 
R=1-—D;2Dyun et R = 1 — Dxs1D22 


A.3. Composition de LFTS : 


Considérons le système d’interconnexion entre les trois blocs des transferts présentés ci- 


dessous. Sachant que chaque bloc possédant un model LFT incertain (voir la Fig. A.5) : 
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Fig. AS. Interconnexion entre trois systèmes incertains 


En faisant quelques simplifications sur le schéma précédent et rassemblant tous les transferts 
connus dans un seul bloc P et les autres transferts inconnus (perturbations) dans le bloc A;. On 


obtient donc le système LFT présenté dans la Fig. A.6 comme suit : 


A, O0 oO 
0 A, 0 
0 O0 A, 


Fig.6. Système LFT équivalent 
A partir de la Fig. A.6, on peut envisager les opérations ci-dessous : 
A.3.1. Operations sur LFT : 
Supposons qu’on a deux systèmes incertains donnés par les réalisations minimales ci-dessous : 


A1 B: A2 | (AT) 


G1(A1) = É sa ;  G2() = e D, 
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D'où, la matrice des incertitudes A peut être donc définie comme suit : 


(A.18) 


a. Produit des LFTSs : 
Le système résultant P.(A) fourni par l’association en cascade des deux systèmes incertains 
G:(A;) et G>(A>) représente aussi un système LFT dont la réalisation minimale correspondante 


est donnée par: 


À B 
PQ) = (GG) = [ee D] (A9) 
PC Pc 
Avec : 
À B.C B.D 
Ape = ( < À 2); Bpe = ( À Le Di C2) ; Dpe = DiD 
0 À B; 


b. Somme des LFTSs : 
De plus, l’addition des deux transferts G,(A;) et G(A) représente également un système 


LFT défini par la réalisation minimale suivante: 


A B 
Pa (A) = (Gi + G2)(A) = [7° | (A.20) 
Pa Pa 
Avec: 
A 0 B 
Apa = e à) ; Bpa = Fa : Cpa = (C1 C2) ; Dpe = Di+D 
c. Inversion des LFTSs : 
a) Supposons qu’on a la matrice de transfert G telle que: 
Gi | 
G = A.21 
ot 


Sachant que le transfert G..doit être inversible (non singulière). Supposons, maintenant, que 
toutes les incertitudes de modélisation peuvent être englobées dans une seule matrice, notée A . 
Dans ce cas, la LFT inferieure du système d’interconnexion  (P — G — A) peut être écrite 


comme : 
P = F(G,A) (A.22) 
L’inverse de LFT inferieure de l’expression donnée par (A.22) est donc donnée par : 


Pt = (FE,(G,A)) ! = F;(G, A) (A.23) 
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Avec : 

â _ Mi —Mi5 Mi2 (A 24) 

MM Mo — M1 Mii Mia 
b) Supposons maintenant que la matrice G est définie par : 

G = F(P,K) (A.25) 

Avec P , P,2et P;. sont des matrices carrées et inversibles. Alors il est possible d’écrire : 
Z W 

4 = P[?| (A.26) 
Or: 

u=K.y (A.27) 
Sachant que : 

z = F(P,K).w = G.w (A.28) 
Il est possible de résoudre cette équation, pour K = E,(P71,G), comme suite : 

W] _ 112 

Li=P", (A.29) 
Avec : 

z=G.w (A.30) 
Tenant compte la relation : 

u=E,(P"1,G).y (A.31) 
Avec : 

K =EF,(P"!,G) (A.32) 


A.3.2. Fonctions des sensibilités sous formes LFT : 


Supposons le système bouclé représenté sous la forme LFT présentée par la Fig. A.7 : 


Fig. A. 7. LFT inferieure d'un système bouclé 
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Les fonctions des sensibilités du système bouclé de la Fig. A.7 peuvent être exprimées sous 


formes des LEFT comme suite : 


% Fonction de sensibilité directe : 


_ =: — I —(G 
S,=(+GK)I=R (|, El) (A.33) 
* Fonction de sensibilité complémentaire : 


T,=GK(+GK)Ÿ=F, ([ | K) (A.34) 


Page | 142 


43 pui O5 5Y Au 33 Aus 2633 À 3 Cie 

yaital) 

ati js Qle 36 y Linil Ji Yl Aile, À À Lite AG Bas 5 3 985 al Cp Bas Lil Jul > ji, La 5 ENT oûa (à 
Ci AN ui a piball Aid A8 JL, Qi VI 580 UE Y5 las Guatis all sil Aile , jus O5 SY Gags 
GAsi juisi Jleriul Gi jeule Off dlariuls À 5 jadll Aulas Aie Ja ne (Gisi Cie ai aps Yi: All 
4 pa Of SN ss jt ile sladll 71 Sul de E, Gus O5 de Y Glaall Gi pull ge 5 3 dhariuls (gi 
dlaxiuls ae Às 5 jee Lulus Joie de a Lasls Aa y Aus GE Jletul Lauls G153i Gil 5531 aie GE 
AA Ci LAN 8e) Gps Vs JsŸ je Alto JiSf at masi ae sil sil si cuis Of Ge sanall 55 
Ai pe Si O5 51 Ge y oo de chasidll sil Asia Le Ce lee chasiall lil 45 lie ei sil 3 5 al 
lili lai ç À lie ua dus ils (puglolikes Cyel ia dl se ie Au 551 Gb ve 4e dhasidl sil 


dainlida Cibals 


iaŸ Ci pull ue 5 38 all Le) Ja € aïla Dabelites Gyal Jia dl js, Ho Al ssidll ; dauluY aKaïll las 3 


INTITULE DU MANUSCRIT : 


ROBUSTIFICATION D'UN CONTROLEUR CONVENTIONNEL PAR UN TRANSFERT 
D'ORDRE NON ENTIER 
Résumé: 


Dans ce travail, une nouvelle stratégie de robustification d’un contrôleur H, primaire basée sur une 
sélection systématique des pondérations irrationnelles adéquate est proposée afin d'améliorer plusieurs 
marges de robustesse, en particulier marge des performances robustes RP. Cette marge doit être garantie 
face aux incertitudes multiplicatives non structurées affectant le modèle nominal de synthèse. Le 
contrôleur H, robustifié est obtenu à partir de la robustification du contrôleur À, primaire tout en suivant 
les étapes suivantes : Premièrement, le contrôleur H, primaire est conçu à partir de la résolution du 
problème de sensibilités pondérée utilisant les pondérations irrationnelles fixes. Ces pondérations sont 
ensuite utilisées afin d’initialiser d’un côté, l’algorithme de robustification proposée. Elles sont également 
utilisées pour définir l’espace de recherches à d’autres pondérations fractionnaires ajustables. La sélection 
optimale de ces dernières est assurée par l’algorithme d’optimisation d’essaim de particules PSO de telles 
sortes que certaines règles de réglage existant dans la littérature doivent être respectées. Comme la 
méthode H,, est adaptée uniquement pour les pondérations d’ordre non entier, ce qui peut poser des 
problèmes sérieux avec celles d’ordre non entier, qui sont sélectionnées préalablement par 
l’algorithme PSO. A cet effet on doit recourir une approche basée sur l’identification fréquentielle avec 
laquelle toutes les pondérations fractionnaires ajustables peuvent être approximées par des filtres d’ordre 
entiers. Ces derniers permettant de reformuler un nouveau problème de sensibilité pondérée qu’avec 
lequel une représentation d’état d’un nouveau contrôleur est obtenue via sa solution optimale. Ces étapes 
sont réitérées en tant que la condition sur RP n’est satisfaite. L’algorithme de robustification est donc 
achevé par l’obtention du contrôleur H,, robustifié assurant une bonne marge de robustesse de RP non 
seulement pour le régime de fonctionnement nominal du système physique, mais aussi pour un ensemble 
englobant un nombre très élevé des régimes de fonctionnement perturbé. La validation de cette nouvelle 
proposition est montrée à partir de la commande d’un moteur synchrone à aimants permanents où son 
comportement réel est modélisé par un modèle linéaire incertain. Les résultats de simulation obtenus sont 
comparés dans le plan fréquentiel en utilisant les lieux des valeurs singulières de différentes sensibilités du 
système bouclé. Ces résultats sont aussi comparés dans le plan temporel en utilisant l’environnement 


Powersim. 


Mots-clés : 


x 


Contrôleurs H,X Primaire et Robustifié; Moteur Synchrone à Aimants Permanant ; Algorithme 


Génétique ; Algorithme d’optimisation par l’Essaim Particulaire. 


MANUSCRIPT TITLED 
ROBUSTIFICATION OF A CONVENTIONAL CONTROLLER BY FRACTIONAL 
ORDER TRANSFER FUNCTIONS 
Abstract: 


In this work, a new robustification strategy of a primary H,, controller based on a systematic 
selection of adequate fractional weights is proposed. Its robust performance RP margin is well 
enhanced with respect to unstructured multiplicative uncertainties presented in linear feedback 
systems. The proposed robustification strategy provides à robustified H, controller when 
following proposed hierarchical control is well respected. First, the primary A, controller is 
synthesized from solving a weighted-mixed sensitivity H, problem using initial integer weights. 
Thus, an initial RP margin is obtained. Second, an automatic selection of adjustable fractional 
weights is performed by the particle swarm optimization algorithm, in which some proposed 
tuning rules are accordingly well satisfied. Third, frequency response data of these weights are 
computed and then fitted by corresponding approximated integer weights using a frequency 
identification technique. Finally, these weights reformulate a new weighted-mixed sensitivity 
problem. The optimal solution to this problem updates the previous initial RP margin. These last 
three steps are repeated as the updated RP margin is diminished. Otherwise, the proposed 
hierarchical control is achieved by selecting the best adjustable fractional weights, providing, 
therefore, the best approximated integer weights and leading, therefore, to the robustified H,, 
controller. In order to confirm the effectiveness of our proposed hierarchical control, primary and 
robustified H  controllers are applied on a permanent magnet synchronous motor where its 
actual behavior is modeled by an unstructured multiplicative uncertain model. The results 
obtained are compared in frequency domains using the singular value plots of their sensitivity 
functions. Otherwise, the same results are compared in time domains using the Powersim 


software. 


Key-words 
Primary and Robustified H , controllers; Permanent Magnet Synchronous Motor; Genetic 


Algorithm; Particle Swarm Optimization algorithm. 


